NASA Technical Memorandum 109029 



NONUNEAR RANDOM RESPONSE PREDICTION USING 


J. H. ROBINSON 

NASA LANGLEY RESEARCH CENTER 
HAMPTON, VIRGINIA 


C. K. CHIANG S 

OLD DOMINION UNIVERSITY * « S 

NORFOLK, VIRGINIA 7-5 

S* u o 

O C rH 

S. A. RIZZI z => o 

NASA LANGLEY RESEARCH CENTER 

HAMPTON, VIRGINIA ^ 

m 

O 


X 

a 

o 

z 

< 


OCTOBER 1993 


NASA 

National Aeronautics and 
Space Administration 


Langley Research Center 

Hampton, Virginia 23681-0001 


0 c 

< 

OJ o (X 

z z 

t— • *— i 

J </) nO 
Z 3 

o 

z 

a < 

CO 

K < 
OUZ 
<\J w 

0 a 

o UJ 

OQCZ 
a <r 

1 ac 
x uj h- 
►— oO 00 
I Z < 

< a z 

a s 

< OO o 
Z UJ 00 
W* CC X 




Abstract 


n ) A ? eqUlVal ^ nt hneanzatl ™ technique has been incorporated into MSC/NASTRAN to 
predict the nonlinear random response of structures by means of Direct Matrix Abstract 
Programming (DMAP) modifications and inclusion of the nonlinear differential stiffness 
dule inside the iteration loop. An iterative process was used to determine the rms 
displacements. Numerical results obtained for validation on simple plates and beams are 
in good agreement with existing solutions in both the linear and linearized regions The 
versatility of the tntplementation wiH enable the analyst to determine the nonlinear random 
responses tor complex structures under combined loads. The thermo-acoustic response of 

th e h program hema pr ° teCtl0n system P anel is used t0 highlight some of the features of 
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Section 1 Introduction 


The current trends in advanced vehicle development show a need for lighter, more 
economical structural components. This trend, coupled with the increasing propulsion 
and environmental loads associated with these vehicles, has renewed interest in nonlinear 
structural response. This is most evident in, but not necessarily limited to, the aerospace 
industry with such proposed vehicles as the National Aero-Space Plane (NASP) and the 
High Speed Civil Transport (HSCT). Surface panels, particularly those exposed to the engine 
noise and jet exhaust and those in the region of shock boundary layer interactions, are 
anticipated to respond nonlinearily in at least part of the ilight regime. Figure 1 depicts the 
thermo-acoustic loads on a single-stage-to-orbit vehicle. Other intense random loads may he 
transmitted through the structure from engine mounts or other hard points. To effectively 
and economically evaluate these structural components, a practical method of predicting their 
large deflection random response is required. 

There are several methods currently in use to predict the large deflection random response 
of structures. A perturbation method [1] based on classical perturbation theory for nonlinear 
deterministic motion can be used to obtain approximate solutions to weakly nonlinear systems. 
A stochastic averaging method [2] yields approximate solutions when the damping is light 
and the excitation is broadband. This method has been applied principally to single-degree- 
of-freedom systems. The Fokker-Plank-Kolmogorov (FPK) approach [3] is the only method 
that yields an exact solution, but solutions are only available for a few restricted classes of 
problems. The numerical simulation technique, also referred to as the Monte Carlo method 
[4], is the most general method and yields the best results of all the approximate methods. A 
substantial drawback to the Monte Carlo method is the computational time required to solve 
realistic structural problems. The most widely used method is the equivalent linearization 
method [5], It yields good approximate solutions for the statistics of the random response of 
simple and complex structures and lends itself to an incremental solution procedure similar 
to the methods employed in static nonlinear problems. 

The equivalent linearization method for obtaining nonlinear random responses was an 
obvious choice for implementation in a commercial package. The technique has been used, 
refined, and validated by many authors [6 — 10], The validation of the method is well 
documented by many authors for beams, plates, and other nonlinear dynamic structures. 
The refinements include methods for solving structural problems with thermal and acoustic 
loads, initial stresses, and imperfections. Techniques have been developed, for example, for 
the random response of pre- and post-thermally and mechanically buckled plates, linear and 
nonlinear statically deflected panels, and various combinations of concentrated and distributed 
random loads. The equivalent linearization procedure has been applied primarily in research 
or special puipose codes, so a general purpose finite element code incorporating this procedure 
was unavailable. 

The MacNeal-Schwendler Corporation version of NASTRAN (MSC/NASTRAN) [11] 
was selected for this work due to its extensive use in the aerospace and automotive indus- 
tries, where nonlinear random phenomena are most prevalent. The equivalent linearization 



procedure was programmed as a "stand alone" solution sequence for version 67 using the Di- 
rect Matrix Abstraction Program (DMAP) [12] language. It was found that all the necessary 
components of the equivalent linearization procedure already existed as DMAP modules. 
The essence of the new solution sequence therefore consisted of incorporating the necessary 
modules and iterative procedures into an existing MSC/NASTRAN solution sequence for 
linear random analysis. Two solution sequences were available to serve as starting points: 
the Super Element Modal Frequency Response (SEMFREQ) and Super Element Direct Fre- 
quency Response (SEDFREQ) solution sequences. The SEMFREQ was chosen for reasons 
described in this report. 

The large deflection finite element formulation is first reviewed to establish the general 
nonlinear equations of motion. The theory of equivalent linearization is then presented and 
the expression for the equivalent linear stiffness is derived. An overview of the iterative 
implementation of the equivalent linearization procedure is presented in flow chart form with 
consideration to the various methods of solving dynamic systems. The ease with which the 
expression tor the equivalent linear stillness is evaluated in multi-degree-ol-lreedom systems 
is somewhat dependent on the method used to form and solve the equations of motion. 
The evaluation of the equivalent linear stillness and the particulars of programming the 
new solution sequence are presented for broad-band Gaussian loads and modal equations 
of motion. In the validation section ol this report, textbook examples are used to compare 
the MSC/NASTRAN equivalent linearization solution sequence with published results. A 
series of simple plate problems are presented to show potential users how to use the solution 
sequences to solve a variety of problems. A final example problem is shown to demonstrate 
the ability of the solution sequence to efficiently solve complex structural problems. 



Figure 1: Design environment for a generic hypersonic vehicle 
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Section 2 Theory 

The large deflection nonlinear finite element formulation is first reviewed tor the determi- 
nation of the system matrices. The equivalent linearization technique is then introduced for 
the solution of the nonlinear equations of motion. Several special cases are then considered 
for the determination of the equivalent linear stiffnesses. 


2.1 Large Deflection Finite Element Formulation 

The large deflection nonlinear strain-displacement relationships as taken from elasticity 
[13] are: 
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where u, v, and w are the three displacements, e x , fy. and r z are the normal strains, and 7 X y, 
7 XZ , and 7 yz are the shear strains. All are functions of x, y, and z. 

For an arbitrary finite element, assume nondimensional shape functions, N„, such that 
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where {u t u 2 ••• u n }, {v, v 2 ••• v„}, and {w, w 2 ••• w„} are the vectors of nodal 
displacements for the n nodes of the element. 

The matrix form of equation (1) is [14] 
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The subscripts L and N denote the linear and nonlinear part of the total strain, respectively, 
and superscript T denotes transpose of a quantity. 

The variation of the strain, {<r}, is expressed as 
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The matrix [I] in equations (10) and (11) is a (3x3) identity matrix and the vector {q} is 
the vector ol nodal displacements. Note that the shape lunction and displacement vector are 
dependent on the particular element chosen. 

The internal force is computed from the equation of static equilibrium. 


I [B] T M 


and the variation of the internal force is 


SF = J [B] T {«W}dV + J [<SB ] 1 {rr} dV 

V V 

Substituting equation (5) and the stress -strain relation, 

{M = [D]{fc} 

into equation (13) and using the identity, 

= [G] T [dH] T {<7> 

= |C:] T [r](M} 

= [G] T M[GH«q) 

yields a simple expression for the variation of internal force. 

(SF = k + k 1 { q } + k2{q}{q} T l {<Sq} 


The quantities k, k 1 { q } , and k2{q){q) T in equation (16) are given as 

M = J [Hi.] t [d][h l ] dv + [kg) 
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In equations (17 — 19), [D] is the material property matrix, [k] is the linear stiffness matrix, 
kl and k2 are the nonlinear stillnesses, and [k g ] is the geometric stiffness (which depends 
on the initial stresses). 

The element internal force vector {7} is defined as 


{7} 


[k]+[kl{q}]+ [k2{q}{q} T ]]{q} 


and the system internal force {T} is 


(22) 


{r} = [[K] + [Kl {Q}] + [k2{Q}{Q} t ]]{Q} 


(23) 


The system mass and damping matrices are obtained using the standard finite element 
formulation [15]. 

The equation of motion based on the nonlinear strain-displacement relations is 


[M]{Q} + [C]{Q] + [[K] + [K 1 { Q}] + [l<2{Q}{Q} T ]]{Q} = { P } (24 ) 


or, in more general form, as 


[M]{q| + [C]{q] + {r(Q, Q 2 ,Q’)} = {P} (25) 

where the matrices [M], [C], and [K] are the system linear mass, damping, and stiffness 
matrices. The vector {P} is the time dependent load and Kl and K2 are the system first- 
and second-order nonlinear stiffnesses. 

Equation (25) has no general solution when the excitation is random. An approximate 
solution to these equations is obtained by seeking an equivalent linear system [6], of the form 

[M]{Q} + [C]{Q} + [K e ]{Q} = {P} (26) 

where [K^] is an equivalent linear stiffness matrix. 
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2.2 Equivalent Linear Stiffness Matrix [K e ] 

The equivalent linear stiffness matrix [K c .] is to be determined such that the difference 
between the actual nonlinear system and the approximate linear system is minimized. The 
approach may be thought ot as a statistical version of a classical least square minimization. 
The error in obtaining the approximate system is defined as 


{A} = {r}-[I< e ]{Q} 


(27) 


Since the error is a random function of time, the required condition is that the ensemble 
average or expectation ot the mean square error be a minimum. This is expressed as 


e[{A}(A) t ' 


minimum 


(28) 


where E[ ] denotes the expectation operator. As in the cases of classical least square 
minimization, the necessary condition for satisfying equation (28) is 


dE 


{A}{A} 


d[K t 
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Substituting equation (27) into equation (29), and interchanging the expectation and differ- 
entiation operators yields 


e[{I’HQ) T ] = e[{Q}{Q} t ][KJ t (30) 

Using the fact that the matrix E[{Q}{Q} T ] is non-singular, the equivalent linear stiffness 
matrix [K e ] can be determined from the equation 


[!<«]= E[{QHQ} T ] _l E[{r}(Q) T ] (31) 

The equivalent linear stiffness [K e ] defined in equation (31) can be directly programmed in 
a finite element code if the stiffnesses K1 and K2 are available and the expectation operator 
can be evaluated. 

Two assumptions regarding the distribution and dependance of the displacements are 
necessary in order to evaluate equation (31). The most commonly assumed distribution of the 
displacements is a Gaussian distribution, since the most commonly encountered random loads 
are typically modelled by Gaussian distributions. The most commonly assumed dependence 
between displacement responses is that they are independent. This is simply because, in a 
linear modal analysis, the modal responses are solved for independently; their modes are 
uncoupled. These assumptions are not the only possible assumptions; other assumptions can 
easily be substituted, but would yield more complicated results. 

It is generally assumed that the response is Gaussian if the load is Gaussian. By using 
the formula for the expected value of a Gaussian vector {//} 

E[f(?/) '/] = E jv7 // T } E{V f(/?)} (32) 
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where V is the gradient operator, E[ { r } { Q } T ] on the right hand side of equation (31) 
becomes 
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. The equivalent linear stiffness matrix [KJ can then be determined from the equation 
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where [K e ] is an equivalent linear function of the displacement vector { Q } , which is one 
order less than the nonlinear system stiffness matrix {T}. 

The nonlinear stiffnesses are generally formed in tangential or differential form and 
the expectation operator in equation (35) requires knowledge of the joint probability density 
function ot the response vector, which is the unknown. Therefore, the equivalent linearization 
solution procedure is programmed in an iterative fashion and some additional assumptions 
regarding the expectations of the response vector are required. It should be noted that, if K1 
and K2 are available, the mean square response can be obtained directly [7] with appropriate 
assumptions for the expectation operator. 

In all instances cited above, assumptions regarding the expectations of the response 
vector are required. These assumptions are usually based on a knowledge of the excitation 
and the solution method used. A discussion of the general iterative equivalent linear solution 
procedures is next presented. 


2.3 Iterative Equivalent Linearization Solution Methods 

There are two basic means to solve linear dynamic equations of motion: one uses the 
physical degrees of freedom and the other uses the modal degrees of freedom. The first 
method is generally referred to as the direct frequency response method and requires solving 
a complex coupled system ot equations in the nodal degrees of freedom at each frequency 
of interest. The second method is generally referred to as the modal frequency response 
method. It involves solving for the linear eigenvectors first and transforming the equations 
of motion into modal coordinates. The resulting system of equations is uncoupled and can 
be easily solved at each frequency of interest. 

The primary consideration as to which method to use for a particular linear system 
is based on the computational time required. This decision in an equivalent linearization 
solution procedure is further complicated by the iterative nature of the problem and the 

evaluation of the expectations. The choice of method can either greatly simplify or complicate 
the process. 
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The direct method would seem to be the easiest and most straightforward to implement, 
and the computational time required would be simple to compute. The difficulty in the direct 
method arises in the assumptions regarding the expectation operator in the expression for the 
equivalent linear stiffness and the implementation of these assumptions in a general sense. 
Accurate approximations of the expectation operator require assumptions regarding the full 
set of four moments (mean, standard deviation, skewness, and kurtosis) of the response vector 
in nodal degrees of freedom. It should be noted that in physical coordinates, the correlations 
between all the degrees of freedom are necessary and must be determined. 

As a simple example of the direct method, consider a beam of length L with ten nodes 
and three degrees of freedom, u, w, and 0, at each node. The evaluation of equation (31) for 
the equivalent linear stiffness requires the evaluation of the complete set of expectations of 
all the nodal degrees of freedom to the fourth moments. The equivalent linearization solution 
relies on determining expressions for the third and fourth moments in terms of the first and 
second moments. These may be obtained by assuming appropriate probability distributions 
for the nodal displacements. In the beam example, if the excitation is broadband, Gaussian 
distributed, and spatially correlated over the beam, it can be assumed that the responses w 
and 0 are Gaussian and u is Chi-square distributed. From these assumptions, an expression 
for the equivalent linear stiffness in terms of the first and second moments of the response 
can be found. However each entry in the 30 x 30 equivalent linear stiffness matrix could 
have a different coefficient representative of the degrees of freedom, correlation coefficients 
between the degrees of freedom, and the order of the expectations involved. The complexity 
in using physical degrees ot freedom can be deduced from this simple problem when it is 
noted that it is terms such as the square ot the slope and the in-plane displacement that 
are strongly correlated. This entire process is programmable, but it is not easily done in a 
general sense. The selection of modal coordinates will be seen to make the evaluation of 
equation (31) simpler. 

The modal solution method of the equivalent linearization procedure is simpler to 
implement than the direct method because reasonable assumptions regarding the correlation 
of the modal degrees-of-freedom as well as their joint distribution are possible. This is 
not to say that the modal approach is without deficiencies or difficulties. To illustrate the 
advantages and difficulties with the iterative modal solution procedure, the simple beam 
problem discussed in the direct method is used. The first difficulty arises immediately from 
the linear eigenvalue problem. The extracted eigenvectors are functions of either the out-of- 
plane nodal displacement (bending modes) or the in-plane nodal displacement (membrane 
modes), but not both. This is because the bending motion of the beam is coupled to the 
membrane motion through the nonlinear terms. 

There are three ways to handle the decoupling of the membrane and bending motion 
induced by the use of the linear eigenvectors. The first way is to simply exclude the 
membrane modes from the modal response. This is easy, but not particularly accurate. 
A popular corollary to this solution is used for one- and two-dimensional structures [7]. This 
procedure assumes the in-plane inertia and damping to be negligible. It is then possible 
to solve for the membrane modes in terms of the bending modes and thus account for the 
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in-plane stiffness. This procedure is efficient, but highly specialized and difficult to include 
in a general finite element code. 


The second method involves selecting particular bending modes and membrane modes 
to include in the formulation. The difficulties that arise from this solution are similar to those 
encountered in the direct method when trying to evaluate the expectations and solving the 
system of equations. In the beam problem, it is again assumed that the bending is Gaussian 
- and the membrane is Chi-square distributed when the excitation is Gaussian. The bending 
modes can be assumed uncorrelated with respect to each other, as can be the membrane 
modes, but the membrane modes are strongly correlated to the square of the bending modes. 
The resulting system of equations is coupled and the expression for the equivalent linear 
stiffness matrix is only marginally simplified with respect to the direct method. Another 
difficulty with the linear modal solution procedure is that the type of modes, bending, 
membrane, or otherwise, are not always readily identifiable or available. Many current finite 
element programs use Lanczos-type eigenvalue solvers in which only the lowest modes or 
modes within a certain range are computed. It is difficult to construct a general program 
using this method that will extract the particular eigenvectors needed for an accurate solution. 


The third modal solution method for the equivalent linearization procedure uses updated 
or “equivalent linear” modes. The obvious drawback to this method is that it requires the 
eigenvalue problem to be solved at each iteration. The advantages of this method are that 
the system of equations that are solved at each frequency are uncoupled and that simple 
assumptions regarding the moments and correlation of the modal responses are adequate for 
accurate solutions. The simple beam problem discussed in the previous solution methods 
could be solved with only a small number of updated modes. If the load were Gaussian, 
these modes could be assumed Gaussian-distributed and uncorrelated. Although the means of 
the equivalent linear modal amplitudes are also assumed to be zero, this does not require that 
all the nodal displacements comprising the mode shape have zero means. The relationship 
between the mean of the in-plane displacement, u, and the mean square of the slopes, 0, in 
the simple beam problem, is implicitly maintained in the equivalent linear modal approach. 


The relationship between the steps involved in the direct, linear modal, and equivalent 
linear modal approaches to implementing the equivalent linearization solution procedure are 
outlined in the flowchart in figure 2 for a general finite element program. The solution 
procedure is iterative as previously discussed, since the nonlinear stiffness is only available 
in a differential form. The convergence of the iterative procedure is based on the Euclidian 
norm of the response variance vector. It was determined that the equivalent linear modal 
method of solving the iterative equivalent linearization procedure would be the simplest and 
most versatile of the three methods to implement in MSC/NASTRAN. 
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Figure 2: Flowchart for equivalent linearization solution procedure 


2.4 Implementation 

The equivalent linear stiffness matrix in equation (31) must first be expressed in equiv- 
alent linear modal coordinates in order to evaluate the expectation operator. The stiffness 
vector {T(Q, Q 2 , Q 3 )} in equivalent linear modal coordinates has the form {F(A, A 3 )}, 
where the bar indicates a quantity transformed in modal coordinates. The expression for the 
equivalent linear stiffness, equation (35), with the Gaussian, zero mean, and uncorrelated 
modal response assumptions reduces to 


[k7| = e 


a{r} 

5{A} 


(36) 


The partial derivatives are easily performed and yield a linear modal stiffness matrix and 
a differential modal stiffness matrix that is based on the mean square of the modal response 
[5]. The modal representation of the equivalent linear stiffness is then 

[KT] = [K] + 3[K2(E[A 2 ])] (37) 
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This model expression is not directly programmable in MSC/NASTRAN. It must instead be 
expressed in physical coordinates, since the eigenvalue problem in MSC/NASTRAN is solved 
in physical coordinates. In addition, the differential stiffness matrix in MSC/NASTRAN 
is formed using the physical displacements. The linear stillness matrix in equation (35) 
in physical coordinates is simply the linear stillness matrix as assembled and computed 
in the MSC/NASTRAN program. The dilterential stiffness matrix expression in physical 
coordinate is the MSC/NASTRAN differential stiffness matrix formed using an equivalent 
linear displacement vector. This equivalent linear displacement vector is given by 

{Q> = fflW} + [5] {//•«} (38) 

wherej <r A } is a vector of the standard deviations of the equivalent linear modal amplitudes 
and [$] is the matrix of normalized eigenvectors. The standard deviation of the modal 
amplitudes is always positive. The sign convention of the physical displacement is determined 
by the eigenvectors. The vector {//.„} is the mean displacement obtained from a static 
solution sequence. The matrix of eigenvectors is normalized such that the magnitude of each 
eigenvector in the matrix is unity. The final expression for the equivalent linear stiffness 
is then 


[K*] = [K] + 3[K R ] ( 39) 

where [ Kr ] is the standard MSC/NASTRAN differential stiffness matrix. 
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Section 3 Programmer’s Notes 

The MSC/NASTRAN version 67 solution sequences are written using a common set 
of “subroutines” or SUBDMAPs. It is the MAIN SUBDMAPs, “main programs,” that 
vary significantly and usually contain the essence of the solution procedure. The authors 
attempted to follow this structure in the development of the new Super Element Modal 
Equivalent Linear Random Response (SEMELRR) solution sequence, but some alterations 
to the common SUBDMAPs were also necessary. These alterations to the SUBDMAPs, as 
• well as a description of the MAIN SUBDMAP of the SEMELRR solution sequence, are 
outlined. 

All solution sequences are broken down into three general sections. These sections are 
simply expressed as Phase 1, Phase 2, and Phase 3. The Phase 1 portion of the program is 
dedicated primarily to the setting-up of the problem and the assembly of the linear matrices. 
Key portions of these procedures are the reading of the NASTRAN data deck, the restart 
capability, the creation of the element summary tables, the partitioning of the global degrees 
of freedom into the various analysis set tables (USET, etc.), and the formation and assembly 
of the linear elements and their reduction to the analysis set. The Phase 2 procedures are 
primarily associated with the actual solution of the problem. These solution procedures 
are, for example, the eigenvalue and eigenvector extraction routine of SOL 103, the linear 
matrix equation solvers in SOL 101, and the modal matrix formation and complex frequency 
response solver routines in SOL 111. SOL 106, the nonlinear static solution sequence, has 
a complicated Phase 2. This Phase 2 involves an iterative solution procedure similar to the 
solution sequence that was written into the Phase 2 of the SEMELRR solution sequence. 
Phase 3 procedures are primarily associated with post-processing routines such as data 
recovery, plotting and printing, and stress/strain calculations. Phase 3 also includes the 
dynamic sensitivity analysis. The calculation of power spectral densities and root mean 
square responses for random analysis using SOL 111 and SOL 108 are also included in 
Phase 3 procedures. The scattered placement of these procedures caused difficulty in the 
implementation of the equivalent linearization solution procedure. 

3.1 SEMELRR Main SUBDMAP 

As a starting point from which to write the SEMELRR DMAP, the authors selected the 
MSC/NASTRAN-delivered SOL 1 1 1 main SUBDMAP. This solution sequence is capable 
of performing linear random analysis. The primary additions to this solution sequence were 
envisioned to be the incoiporation of Phase 2 procedures, similar to those found in SOL 106, 
for the formation of the nonlinear stiffness matrices and the iterative solution method. It was 
immediately apparent that the logical flow of the set of MSC/NASTRAN SUBDMAPs and 
modules did not readily permit simultaneous geometric nonlinearities and dynamics. The 
SEMELRR Solution sequence would have to be a hybrid-type solution sequence comprised 
of linear and nonlinear Phases. The calculation of the rms quantities, which usually occurs 
in Phase 3, and the necessity of having that information available in the iterative procedures 
required the new solution sequence to have no clear distinction between Phase 2 and Phase 3. 
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Figure 3: SUBDMAP call tree of SEMELRR solution sequence 
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The logical How ot MSC/NASTRAN solution sequences is partly controlled by parame- 
ters and flags that are set in the main SUBDMAP of the solution sequence. These flags are 
passed to the standard SUBDMAPs and appear in IF, THEN, ELSE type logical construc- 
tion. Typical character parameters are the solution type (SOLTYP= “DIRECT,” “MODAL,” 
etc.), solution approach (APP= “STATIC,” “FREQRESP,” “TRANRESP,” etc.), and logical 
flags are (NONLNR, AERO, FS, etc.). The logical flow for a nonlinear dynamic solution 
sequence does not exist, but by changing these parameters during the solution procedure the 
necessary logical flow can be generated. 

The SEMELRR solution sequence requires that linear and nonlinear element tables be 
generated in Phase 1 procedures and that linear dynamic data recovery be performed in Phase 
3 procedures. In order to generate the necessary matrices and tables for both geometric 
nonlinear and linear dynamic procedures, the pre-processing sequence Phase 1 was initiated 
with the NONLNR flag set to TRUE in the call to SUPER1, Figure 3, and the APP (approach) 
was set to FREQRESP. The reader is referred to the MSC/NASTRAN user’s manual [11] 
for a description of these parameters. 

The logical flag NONLNR was set TRUE for Phase 1 only, in the call to SUPER 1, and 
not lor Phase 2 or 3, because linear data recovery is required in SUPER3. In addition to 
this modification, it was required that the element summary table, ESTL, needed for linear 
dynamic analysis (not generated when NONLNR is TRUE) be equivalenced to the element 
summary table, EST, for the linear portion of nonlinear analysis (generated when NONLNR 
is TRUE). This equivalence was programmed as an ALTER to the SEMG SUBDMAP. 

The programming of the SEMELRR main SUBDMAP consisted of writing an iterative 
procedure around the frequency response solution procedures, the geometric nonlinear matrix 
generation procedures, and the data recovery SUBDMAP, SUPER3, which includes the 
updated displacement calculations. To implement this iterative procedure, some of the files 
needed tor the next iteration have to be saved. The module FILE to save or overwrite 
files was used for this purpose. Phase 3 procedures were included in the iteration loop 
because the updated displacements, necessary as input to the differential stiffness modules, 
are obtained from SUBDMAP SEDRCVR in Phase 3. SUBDMAP SEDRCVR had to 
be substantially rewritten to generate the correct updated displacements for the equivalent 
linearization procedure, equation (38). The calculation of the updated displacements will 
be discussed in depth in the following subsection. A full listing of the SEMELRR main 
SUBDMAP is provided in Appendix A. 

The formation of the geometric nonlinear stiffness matrix in Phase 2 follows closely with 
the procedure in Nonlinear Transient solution sequence (NLTRAN, SOL 129). The linear 
dynamic equations of motion are solved first and the linear rms displacement vector {A} 
is obtained. II the parameter LGDISP is greater than —1, the geometric nonlinear stiffness 
matrix KDJJ is formed from module EMA on the next iteration by applying this linear 
displacement vector. This geometric nonlinear stiffness matrix then reduces to KDLL, [Kr] 
in equation (39). (If the parameter LGDISP equals -1, only the linear frequency response 
is calculated.) The equivalent linear stiffness matrix [KJ now consists of two matrices: the 
linear stiffness [K] and the differential stiffness matrix [Kr]. The frequency response is then 
obtained using both geometric nonlinear and linear dynamic matrices. 
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This iteration method can be used to determine the rms displacements; however, it 
is slow to converge. An improved method for speeding up the convergence is to use an 
underrelaxation approach where displacements are not updated to their full values, but instead 
to the scale of the full values after each iteration. This method can be expressed as 


(Qlc„„„. = (i - 0HQ), < 4 °) 

A user-defined convergence enhancement parameter, ft (BETA in DMAP programming), is 
introduced to scale the updated displacements. If the nonlinearity is mild to moderate, the 
convergence of the iteration procedure is faster for 0.5 < BETA < 1.0. It the nonlinearity 
is severe, the convergence of the iteration procedure is faster for 0.0 < BETA < 0.5. The 
parameter BETA is set by the user in the Bulk Data Deck. 

Two user-defined parameters were introduced to control the termination of the iterative 
loop. The user-defined parameter MAXITER defines the maximum allowable number of 
iterations and the user-defined parameter MAXNORM sets the convergence criteria, i.e. 


{Q} 


current 


{Q} 


previous 


error < MAXNORM 


(41) 


If the iteration count exceeds MAXITER or if the error norm, equation (41), is less than 
MAXNORM, the solution sequence will terminate normally. There is a warning message if 
the solution is not converged after the MAXITER iterations. There are two ways to handle 
convergence errors; the first is by increasing the number ot allowable iterations, MAXITER, 
and the second is by choosing a different convergence enhancement parameter BETA, which 
is less than the previous BETA. A summary of the user-defined parameters and defaults is 
given in Appendix B. 


3.2 Updated Displacement Calculation 

The updated displacement vector is formed by multiplying the maximum rms displace- 
ment by the updated mode shapes. In order to do so, one deflection point number has 
to be obtained first by asking for XYPRINT (or XYPLOT) in the Control Deck ot the 
MSC/NASTRAN data cards. In the SUBDMAP SEDRCVR, individual modes of the actual 
displacement vector are extracted. For each mode, the modal rms responses are calculated 
from module RANDOM. Each mode is normalized to unity for the largest component ot the 
eigenvector. The actual rms response of each mode is then obtained by multiplying the rms 
response by the normalized eigenvector. The updated response of the structure can be calcu- 
lated by using superposition of the modes and storing the updated rms displacement vector. 
This procedure entails the assumption that the modes and modal responses are independent. 
The modified SUBDMAP SEDRCVR is in Appendix C. 

Although some minor modifications on SUBDMAP SUPER3 are made, no functional 
procedure was carried out. The modification passes parameters needed for communication 
between the main SUBDMAP and the SEDRCVR SUBDMAP. The modified SUBDMAP 
SUPER3 is included in Appendix D. 
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3.3 Output Requests 

One new feature from the output request is for plotting the overall rms displacement 
output. In module RANDOM, only the rms values for a single degree of freedom are 
calculated. The actual overall rms displacement is formed by the updated mode shape 
at each iteration. Therefore, at the converged stage, the overall rms displacement can be 
extracted by using a DISP=ALL card in the Control Deck. 

There is no rms strain response obtained lrom frequency random analysis of SOL 111. 
If rms element strain is required, the user-defined parameter RMSTRAIN has to be set to 
1. For this case, the STRAIN=ALL card is needed in the Control Deck and the strains will 
be calculated. 
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Section 4 Validation of The SEMELRR Solution Sequence 

The linearized random vibration capability developed for use with MSC/NASTRAN 
is validated by solving lour problems and comparing the results with known solutions. 
The frequency response ot tree vibration and rms displacement response of forced random 

vibrations of a plate and a beam are considered. The results show that reasonable accuracy 
is achieved. 

Problem 1: Random response of beams 

The rms displacements ot a 12-in x 2-in x 0.064-in aluminum beam with either 
end clamped or simply supported and subjected to uniformly distributed random loads is 
investigated. To demonstrate the accuracy of the SEMELRR results, approximate rms 
maximum deflections were obtained by using a separate Equivalent Linearization (EL) 
analysis [16] and finite element (FE) solution [7], Results are shown in Figure 4. Since 
all three results (EL, FE, and SEMELRR) based on small deflection linear theory lie directly 
on top ot one another, it is shown as one straight dotted and dashed line. The EL and FE 
results are identical tor the linearized case, so one curve is plotted for these two methods. For 
acoustic excitations less than 90 dB for a simply supported case and 1 10 dB for a clamped 
case, the small deflection assumption yields good results. At high SSL, however, the small 
deflection theory overestimates the rms deflection, while it underestimates the frequency of 
vibration. It is cleaily demonstrated that the SEMELRR results give reasonable predictions 
as computed to the EL and FE solutions in both linear and linearized cases 



Figure 4: Eltect ot acoustic excitation level on maximum deflection for beams. 
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Problem 2: Random response of a clamped plate 

The comparison is made for rms displacements as a function of SSL of an aluminum 
plate [7]. For acoustic excitations less than 1(H) dB, the small dellection assumption yields 
good results as shown in Figure 5. Above 100 dB, the large deflection formulation must 
be used. At the 130-dB level, the results between the SEMELRR and Locke’s analysis [7] 
show a 6-percent difference. The discrepancy is attributable to the approximation of the 
nonlinear stiffness matrix in equation (39) and the assumption of curvatures and midsurface 
strains in the von Karman sense in Locke’s formulation. 



Figure 5: Effect on acoustic excitation level on maximum dellection for clamped plate. 


Problem 3: Free vibration of rectangular plate 

The free vibration of a 15-in x 12-in x 0.04-in aluminum rectangular plate reported in 
Chiang’s paper [17] is used. The variation in SEMELRR free vibration results of a plate with 
all edges clamped for the frequency ratio uj / uj 0 at different amplitudes is shown in Figure 6. 
u> 0 is the fundamental frequency of the clamped plate. There is a maximum of 10-percent 
difference between the SEMELRR and Chiang’s results. The frequency ratio for Chiang’s 
results are lower. The differences are caused by two factors. First, Chiang’s formulation used 
von Karman strain-displacement relations, which use thin plate assumptions, and therefore do 
not have all the terms in equation (1). The second is due to the approximation in equation (39). 
Because of this approximation, in which the first-order stiffness matrix in the SEMELRR is 
calculated one more time than the equivalent linearization approach, the linearized frequency 
is expected to be higher. These results show the SEMELRR procedure gives reasonable 
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Modal Amplitude 


predictions 
Figures, 4, 


in comparison to finite element 1 17, 7] and equivalent linearization [16] solutions 
5, and 6. 



Figure 6: Amplitude versus frequency ratio for clamped plate 
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Section 5 Example Problems 

This section is intended to provide a series of simple analyses that demonstrate the 
capabilities and use of the equivalent linearization solution procedure as implemented in 
MSC/NASTRAN. The types ot analyses used in this section were selected Irom a review of 
previously publisher papers [7, 16]. For simplicity these analyses share a common structural 
configuration, that of a simple rectangular panel. The thermo-acoustic response of a large 
hexagonal thermal acoustic protection panel is also presented to further demonstrate some 
of the features of the program. The format of this section follows closely that of the 
MSC/NASTRAN demonstration problems manual [18], It is assumed in this section that 
the reader has a basic understanding of the basic NASTRAN CARDS and DECKS. 

5.1 Problem Execution 

The equivalent linearization solution sequence was written by incorporating portions 
of the NLSTATIC (SOL 106) solution sequence into the SEMFREQ (SOL 111) solution 
sequence. It is assumed in this manual that the reader has a basic understanding of the 
application, options, and limitations of both of those analyses. 

MSC/NASTRAN performs random response analysis as post-processing to the frequency 
response. The Equivalent Linearization solution sequence is performed by including this post- 
processing in the iteration loop because the rms displacements, which are necessary as input 
to the differential stiffness modules, are obtained in Phase 3 procedure. 

The SEMELRR solution sequence is not included in the MSC/NASTRAN-delivered 
data base, but is available as a separate DMAP program. The program must be read into 
the Executive Deck of the NASTRAN data file, and the individual SUBDMAPs and main 
SUBDMAP must be compiled and linked as part of each execution. The solution sequence 
can also be incorporated into the NASTRAN data base of solution sequences by creating a 
permanent USER. OBJ and USER.EXE as discussed in Chapter 7 of the “DMAP and DATA 
BASE APPLICATIONS” seminar notes [19], The necessary commands to include, compile, 
link, and execute the solution sequence are provided in the example problems. 

Model Description 

The basic Bulk Data cards for the rectangular panel will be included in each example but 
will appear only here. The demonstration Bulk Data Deck includes the CQUAD4, GRID, 
SPC1, PSHELL, MAT2, and MAT8 cards. The rectangular aluminum plate is 12-in x 15-in 
x 0.04— in and is modeled using 64 QUAD4 elements with inplane and bending material 
property entries on the element PSHELL cards. The boundaries are assumed completely 
clamped. The zero displacements and rotations are enforced using SPC1 cards. 
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5.2 Linear Random Analysis 

Problem Description 

The random response of the plate subjected to broad-band acoustic excitation is first 
demonstrated. The spectral density functions of the selected displacements and element 
stresses are computed. 

Executive Control Deck 


The Executive Control Deck specifies that Structured Solution Sequence 1 1 1 (Modal 
Frequency Responses) is to be used to analyze the plate response under random loads. 


ID PLATE.DEMD $ 
SOL IDS 
TIME 10 J 
CEND S 


Case Control Deck 


METHOD 

FREQ 

RANDOM 
LOADS ET 

DLOAD 

SDMAP 


Specifies method by which the eigenvalues and eigenvectors will be 
extracted. 

Selects the set of frequencies to be solved in frequency response 
problems. 

Random Analysis set selection 

Selects a sequence of load sets referenced by dynamic load cards to 
be applied to the structural model. 

Selects the dynamic load to be applied in a frequency response 
problem. 

Selects table, which defines dampings as function of frequency. 
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SET 1*41 


SET 2 = 1,2, 3,4 ,5 
DISPLACEMENT = I 
STRESS(FIBER) = 2 
ECHO . PUNCH 
TITLE = FLAT PLATE DEMO 
SPC - 5« 

METHOD - 1 
FREQ = 20 
RANDOM - 59 
LOADS ET = 100 
DLOAD - 301 
SDAMP = 400 
$ OUTPUT REQUESTS 
OUTPUTPCYPLOT) 

XYPLC1T D1SP PSDF/ 4 1(T3) 
XYPLOT STRESS PSDF/ 3(3) 
BEGIN BULK 


The rms displacement can only be extracted from the data base PSDF in module 
RANDOM. The first card after OUTPUT of either XYPRINT DISP or XYPLOT DISP 
is needed in the Control Deck. If the rms displacement for the first card is zero, the process 
will stop and the fatal error message will be given. The output from the SEMELRR solution 
sequence is long, since it prints the output information for each iteration. 


Bulk Data Deck 


LSEQ 


DLOAD 

RLOAD1 

FREQ1 

RANDPS 

TABRND1 

TABDMPS 


Defines a sequence of load sets referenced by dynamic load cards to 
be applied to the structural model. In this case, it is used to apply a 
unit pressure load to the plate since the Dynamic Load Scale Factor 
(DAREA) card can only handle the point load. 

Defines a dynamic loading condition for frequency response. 

Defines a frequency-dependent dynamic load for use in frequency 
response problems. 

Defines a set of frequencies to be analyzed. 

Defines load set power spectral density factors for use in random 
analysis. 

Defines power spectral density as a tabular function of frequency for 
use in random analysis. Referenced by the RANDPS entry. 

Defines modal damping as a tabular function of frequency. 
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EIGRL 1 0.1 

LSEQ 100 U 

PLOAD2 400 1 . 

$ 

DLOAD <01 1 . 0 

RLOAD1 204 12 

TABLED1 13 

♦ QR 0.0 1.0 

S 

FREOl 20 00.1 

S 

RAND PS 59 1 

TABRND1 63 

♦TR -1-0 0.0 

♦ TR2 3001.0 0.0 

$ 

TABDMP1 400 

♦DPI 0.0 0.01 

$ 

PA RAM NOCOMPS -1 

$ 


<000.0 B 
400 

1 THRU 64 

1.0 204 

13 

30 on .o l.o ENcrr 

1.0 1500 

1 1.0 Ci 

0.0 d. 4215-5 3000. 0 H . 42 1 S -53000 . 0 0.0 

ENDT 


3000.0 0.01 ENDT 


♦ V'k 


♦TR 

♦TR2 

♦ DPI 


Problem Output 


$ /cy/»t /mco/ni a/nn*cna»/»>K</ct.hf /riai,-r67r<' 

1 NAS TRAN SYSTEM PARAMETER ECHO 


MARCH 29. 199.1 MSC/NASTRAN 0/ u/ o PA ( ; 


NASTRAN SY STEM (14 9) * 1 S THIS IS FOR Raw IO 
NASTRAN RHFFS1 2E* 1 22H9 5 CR1 ASKED FOR THIS 
S 

1 

0 


NSC / NA STRAW 
VERSION - G7R2 
ML 29. 1992 
Ciay Rwt:«arrli Iiw. 
MODEL CRAY Y-MP 
UNI cos (. ,i 
t'RAY-YMP 


MARCH 


u 

0 


NAS TRAN EXECUTIVE CONTROL DECK 


E 


H o 


HSC/NACTRAN 7/2^/‘<2 PACE 


ID CHIANG , NASTRAN 
$ 

SOL 111 $ 

COMPILE SHBDKAP = SEDRCVR SOUIN.MSCSOU NOUST NoREE $ 
ALTER 16R . 169 $ 

ENDALTER $ 

TIME 30 S 
CEND $ 

1 TPS RESULTS ANALYSIS 


MARCH 29, 1991 MSC/NASTRAN 7/29/92 PAGE 


CASE CONTROL DECK ECHO 

CARD 

COUNT 


2 SET 1 * 41 

3 SET 2 = 1,2.1,* ,5 

3 DISPLACEMENT. 1 

* STRAIN ( FI RER) = 2 

G ECHOs PUNCH $ 

7 TITLE * TPS RESULTS ANALYSIS 

H SPC = 56 

9 METHOD s 1 

IV DLOAD *301 

11 FREy * 20 

12 RANDOM * 59 

11 SDAMP □ 400 

14 L0ADSET*100 

15 $ 

16 ( turn rr<XY OUT) 

17 XYPLOT DISP PSDF/ 4UT1) 

1 « s 

19 $ ELEMENT RMS STRESSES OUTPUT 

20 S 

21 XYPLOT STRESS PSDF/ 3(1) 


26 


0 

0 


2 2 XYPLOT STRESS PSDF/ J) ( F. > 

23 XYPLOT STRESS PSDF/ 3(7) 

24 $ 

25 BEGIN BULK 

INPUT BULK DATA CARD COUNT * 
TOTAL COUNT* IB 2 
TPS RESULTS ANALYSIS 


202 


0 

0 SEQUENCE PROCESSOR OUTPUT 


MARCH 29. 1993 MSC/NASTRAN 7/29/92 PAGE 3 


0 THERE ARE 


O ELEN ENT TYPE 


81 POINTS DIVIDED INTO 1 GROUP(S). 
CONNECTION DATA 

NUMBER ASSEMBLY TIME (SEC) 


TOTAL MATRIX ASSEMBLY TIME FOR M ELEMENTS IS 

OORIGINAL PERFORMANCE DATA 

OSUPER (GROUP) ID MO GRIDS AV. CONNECTIVITY C-AVERAUE C-RMS C-MAXIMUM P-GROUPS P-AVERAGE DECOMP TIME (SECS) 

(C.O DOF/GRID ) 

11 0 U.00 0.037 


00 B 1 7.72 

ORESEQUENCED PERFORMANCE DATA 

OSUPER(GROUP) ID NO. GRIDS AV. CONNECTIVITY 


.22 SECONDS 


10. 15 


TPS RESULTS ANALYSIS 


C-AVERAGE C-RMS C-HAX1MUM P-GROUPS P-AVERAGE DECOMP TIME (SECS) HETWiD 

K> . 0 DOF/ GRID ) 

11 0 U.0O U.03C ACTIVE 

MARCH 29. 1993 MSC/NASTRAN 7/29/92 PAGE 4 


9.54 


1 0. 0000000£*00 
TPS RESULTS ANALYSIS 


OLOAD 

T2 T3 

. OOUOOOOE*00 1 . 800000UE*U2 


RESULTANT 

R1 

. U8UU000E*U.< 


. 35l>0lK)UE«03 0 , UUUUOUOE*UU 

MARCH 29 . 1 9 9 3 MSC / NAS TRAN 


LANCZOS PARAMETER VERIFICATION 


INITIAL PROBLEM SPECIFICATION 

DEGREES OF FREEDOM * 24‘. 

NUMBER OF MODES * I 

MODE FLAG = t 

PROBLEM TYPE * 1 

SHIFTING SCALE = 3.407«E*U5 


MESSAGE LEVEL 
OUTPUT UNIT 
SIZE OF WORKSPACE 
MAXIMUM BLOCK SIZE 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 
ACCURACY REQUIRED 


. 94UE-UUU1 
. OOOE* 240 H 

0 . UUU0E*0U 

1. t.Uf.uE-lM 


AFTER PROBLEM SPECIFICATION CHECKING 
NUMBER OF MODES = 1 

PROBLEM TYPE = 1 

SHIFTING SCALE = 3.4i>?hE*0! 

WORKSPACE ALLOCATION 

LANCZOS BLOCK SIZE * 2 

MAX. BLOCK STEPS * 100 

MAX. MODES t 245 


LEFT END Pol NT 
RIGHT END POINT 
CENTER FREQUENCY 


MAX. RITZ VALUES 
MAX. S.O. VECTORS 
WORKSPACE USED 


3 . 94 HE- OOU 1 
1 OOOE* 2 403 
0. 0U00E*00 


2 45 
44159 


CP TIHE ALLOWED 


MAX. TRUST REGION:: 


NUMBER OF USER SUPPLIED VECTORS * 0 

NEW SHIFT * 3.947HE-01 NODES STILL NEEDED * 1 

' USER INFORMATION MESSAGE 5UlU. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE * 3 .947M42E-U1, CYCLES = 1 . OUOOOOE-Ol NUMBER OF EIGENVALUES BELOW THIS VALUE r 

ACCEPTED EIGENVALUES 
2 . 549?E*05 

NEW SHIFT * 5. 2415E*05 MODES STILL NEED El' = 0 

USER INFORMATION MESSAGE 415B STATISTICS FOR SYMMETRIC DECOMPOSITION OF DATA BLOCK SCRATCH FOLLOW 

NUMBER OF NEGATIVE TERMS ON FACTOR DIAGONAL * 1 

USER INFORMATION MESSAGE 5010, STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION 

TRIAL EIGENVALUE * 5 . 241507E*05, CYCLES * 1.152254E+02 NUMBER OF EIGENVALUES BELOW THIS VALUE * 

END OF LANCZOS RUN 

WARNING FLAG ■ 0 

HO. OF MODES COMPUTED s 1 


COMPUTED MODES 

2 . 5496 970502503 BE* 05 
TPS RESULTS ANALYSIS 


MARCH 29, 1993 MSC/NASTRAN 7/29/92 PAGE 


EIGENVALUE ANALYSIS SUMMARY (LANCZOS ITERATION) 


BLOCK SIZE USED 2 

NUMBER OF DECOKPOSITl (*JS 2 

NUMBER OF ROOTS FOUND 1 

NUMBER OF SOLVES REQUIRED 6 


TERMINATION MESSAGE REQUIRED NUMBER OF EIGENVALUES FOUND. 

1 TPS RESULTS ANALYSIS MARCH 29, 1993 MSC/NASTRAN 7/29/92 PAGE 

0 


MODE EXTRACTION EIGENVALUE 

NO. ORDER 

1 1 2 . 549697E*05 

1 TPS RESULTS ANALYSIS 


REAL EIGENVALUES 
RADIANS CYCLES 

5 . 049452E+02 BU3t>453E*01 

MARCH 


GENERALIZED 

MASS 

1 . OOOOOOE+QO 
29. 1993 MSC/NASTRAN 


GENERALI ZED 
STIFFNESS 
2 . 549697E+05 
7/29/92 PACE 
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°*** USER INFORMATION MESSAGE 5222 , UNCOUPLED SOLUTION ALGORITHM USED. 


1 TPS RESULTS ANALYSIS 

0 

1 TPS RESULTS ANALYSIS 

0 

0 XY - OUTPUT SUMMARY 

0 PLOT CURVE FRAME RMS NO. POSITIVE 

TYPE TYPE NO CURVE ID. VALUE CROSSINGS 

0 PSDF DISP 1 411 5* fc. G382f.8E-01 8 . 02873UE+ U 1 

0 PSDF EL STR 2 3 ( 3> 4 . 744B67E* 0.1 ft . 02B73 OE+ U1 

1 * * * END OF JOB • * * 


MARCH 28. 1883 MSC/NASTRAN 7/28/82 PAGE 


MARCH 28, 1883 MSC/NASTOAN 7/28/82 PAGE 


(AUTO OR PSDF) 

XMIN FOR XKAX FOR YNIN FOR X FOR YMAX FOR X FOR* 

ALL DATA ALL DATA ALL DATA YNIN ALL DATA YMAX 

1.000E-01 1 . *i0(JE*0.1 2.867E-1U 1 S00E+03 2.513E-01 «.U10E*01 

1.000E-01 l.MH)E+03 1.516E-02 1.5OOE*03 1.283E+07 ft.010E*Ol 


Figure 8: Output from linear random analysis 


ALTER for Strain 

Executive Control Deck 

The spectral density of the strains can be obtained by using an ALTER of the modal 
frequency response solution sequence (SOL 111). This ALTER is placed in the Executive 
Control Deck and the solution sequence is compiled as shown. 

ID PL ATE, DEMO $ 

SOL 1 1 1 $ 

COMPILE SUBDMAP=SEDRCVR SOUlN=MSCSOU NOLIST NOKEF $ 

ALTER 26 $ 

TYPE PARM,.LN.crrAPE2 $ 

FILE PSDF=OVRWRT $ 

ALTER IK6 $ 

RANDOM XYCDBDR,DrT.MPSDL.OUUV2,OKj2,OQG2 > OSTR2 l OEF2,CASEDR/ 

PSDF.AUTO/S.N.NORAND $ 

END ALTER $ 

TIME 10 S 
CEND S 


Case Control Deck 

With the ALTER in the Executive Control Deck to acquire the strains, only one card 
needs to be changed. The change is on the STRAIN (FIBER) card. To obtain the strain, use 
the XYPLOT STRESS card. Use ol the XYPLOT STRAIN card will cause the compiler to 

produce a fatal error. The other cards are used in the same order as for the linear random 
analysis. 


SET 1 - 41 

SET 2 - 1, 2,3,4, 5 

DISPLACEMENT =. I 

STRAIN(F!BEK) = 2 

ECHO « PUNCH 

TITLE - FLAT PLATE DEMO 

SPC - 56 
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METHOD = I 


FREQ = 20 
RANDOM - 59 
LOADS ET - 100 
DLOAD = 301 
SDAMP - 400 
$ OUTPUT REQUESTS 

output; XYPLOTl 

XYPLOT DISP PSDF/ 41(T3) 
XYPLOT STRESS PSDF / 3(3) 
BEGIN BULK 


5.3 Nonlinear Random Analysis 

Strain and Displacement Spectra 

If the XYPRINT is used in the Case Control Deck, the output spectra with frequency 
increment can be found in the *.f()6 file. The module FREQ1 in the Bulk Data Deck controls 
the starting frequency, frequency interval, and the number of frequency increments. 

Problem Description 

The nonlinear random response of the plate subjected to broad-band acoustic excitation is 
next demonstrated. The spectral density functions of the selected displacements and element 
stresses are computed. 

Executive Control Deck 


The Executive Control Deck specifies that the modified DMAP Modal Frequency Re- 
sponses Structured Solution Sequence is to be compiled, linked, and used to analyze the 
plate response under random loads. 


ACQUIRE NDDL $ 

ID PLATEJ3EMO $ 

COMPILE SEMELRR SOUOUT*USROUT OBJOUT-USROBJ NOLIST NOREF S 
INCLUDE SEMELRR, DMAP $ 

COMPILE SEDRCVR SOUOCT»USROUT ODJOUT-USROBJ NOLIST NOREF $ 
INCLUDE SEDRCVR .DMAP S 

COMPILE SUPER 3 SOUOUT^USROUT OBJOUT=USROBJ NOLIST NOREF $ 
INCLUDE SUPER 3. DMAP $ 

COMPILE SUBDMAP=SEMO SOUIN-MSCSOU NOLIST NOREF $ 

ALTER 32 S 

EQUIVX EST/ESTUALWAYS t 
END ALTER S 
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SOL SEMELRR $ 


LINK SEMELRR $ 
TIME 10$ 

CEND $ 


Case Control Deck 

Use the same Case Control Deck as for the linear random analysis. 
Bulk Data Deck 


Use the same Bulk Data Deck as for the linear random analysis. The following 
Parameters are needed in the Bulk Data Deck to proceed with the SEMELRR solution 
sequence. 


PARAMeters 


LGDISP 

RMSTRAIN 


MAXITER 

ABSNORM 

BETA 


If linearized analysis is performed, set LGDISP=1. (default=-l) 

If rms strain is needed and print no output for STRESS=ALL in 
control deck, set RMSTRAIN=1. 

If rms strain is needed and print output for STRESS=ALL in control 
deck, set RMSTRAIN=2. 

Maximum number of iterations. (default=5) 

Absolute norm for convergence test. 

Convergence enhancement factor, ranging from 0.0 to 1.0, but not for 
0.0. (default=0.5) 


PARAM.RM STRAIN, 1 S IF RMS STRAIN IS NEEDED, RMSTRAIN=l 
PAR AM .MAX ITER, 3 $ MAX. NUMBER OF ITERATION 
PARAM.ABSNORM.2.0E 2 $ ABS. NORM FOR CONVERGENCE TEST 

PARAM.BETAO.S $ SCALE FOR BETTER CONVERGENCE. RANGE FROM 0.0 TO 1.0 {BUT NrTTO.O) 
PAR AM .LGDISP, 1 $ FOR LARGE DISPLACEMENT ANALYSIS, LGDISP*! 


Problem Output 


1 TPS RESULTS ANALYSIS 


KAfiCH 29, 1993 NSC / NAS TRAN 7/29/92 PAGE 


CASE CONTROL DECK ECHO 

CARO 

coirwT 

1 s 

2 DISPLACEMENT . ALL 

3 ECHO = NONE S 

* TITLE « TPS RESULTS ANALYSIS 

5 SPC = 56 

€ METHOD 5 1 

7 ('LOAD *301 

« FREQ * 20 

9 RANDOM * 59 
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SDAMP = 400 
LOADS ET= 100 


10 
11 

12 $ 

13 output* xyout) 

14 XYPLOT PISP PSDF/ 41 |T3) 

15 S 

16 BEGIN BULK 

INPUT BULK DATA CARD COUHT = 
TOTAL COUNT* 187 
TPS RESULTS ANALYSIS 


212 


0 

0 SEQUENCE PROCESSOR OUTPUT 


NARCH 2*. 1883 MSC/NASTRAN 7/29/92 


OTHERE ARE 81 POINTS DIVIDED INTO 1 GROUP (S) 

0 CONNECT! ON DATA 

0 ELEN ENT TYPE NUMBER ASSEMBLY TIME(SEC> 


TOTAL MATRIX ASSEMBLY TIME FOR 
OORIGINAL PERFORMANCE DATA 
0 SUPER (GROUP] ID NO. GRIDS 


0 RESEQUENCE!) PERFORMANCE DATA 
OSUPER< CROUP! ID NO. GRIDS 


TPS RESULTS ANALYSIS 


€4 ELEMENTS IS 0. 22 SECONDS. 
AV. CONNECTIVITY C- AVERAGE C-RMS 
7. 72 ■?* . H?* 10.15 


C-MAXIMUN P'GROUPS P -AVERAGE DECOMP TIME(SECS) 
(6.0 DOF/GRID > 


AV. CONNECTIVITY C-AVERAGE C-RMS C -MAXIMUM P- GROUPS P-AVERAGE DECOMP TIMEfSECS) METHOD 

(€.0 DOF/GRID ) 


0 0 . UO 

MARCH 2 V , 1991 MSC/NASTRAN 


0 . 03 6 
7/29/92 


1 0 . OOOOOOOE+OO 0. UOUOOUOE*! 

TPS RESULTS ANALYSIS 


1 . HOUOUUtl£*U2 1 . 0H00000E*03 


Y>0UOOoE*O.l o . UOOUUOUE*UO 

MARCH 29, 1993 MSC/NASTRAN 


LANCZOS PARAMETER VERIFICATION 


INITIAL PROBLEM SPECIFICATION 
DEGREES OF FREEDOM * 
NUMBER OF MODES * 

MODE FLAG * 

PROBLEH TYPE * 

SHIFTING SCALE * 3. 


407KE-I 


MESSAGE LEVEL 
OUTPUT UNIT 
SIZE OF WORKSPACE 
MAXIMUM BLOCK SIZE 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 
ACCURACY REQUIRE!' 


4 . 94HE-U0U1 
1 . OOOE*24>- 1 
0 . OUUOE-UU 
3 . 696UE-1U 


AFTER PROBLEM SPECIFICATION CHECKING 
NUMBER OF MODES * 1 

PROBLEM TYPE = 1 

SHIFTING SCALE * J.407bE*u! 

WORKSPACE ALLOCATION 

LANCZOS BLOCK SIZE - 2 

MAX. BLOCK STEPS * 11)0 

MAX. MODES * 245 

NUMBER OF USER SUPPL1 El' VECTORS » 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 


MAX. RITZ VALUES 
MAX. S.O. VECTORS 
WORKSPACE USED 


3 . !*4 8E-0001 
1 . 000 E* 246.1 
U.(J000E«0U 


CP TIME ALLOW EX 


MAX TRUST REGION. 1 : - 


NEW SHIFT = 3.947BE-01 NODES STILL NEEDEI' * 1 

0*** USER INFORMATION MESSAGE 5010, STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVA1.UE * 3 . 947642E-01 . CYCLES = l . UOOOOUE-Ol NtIMBER OF EIGENVALUES BELOW THIS VALUE = 

ACCEPTED EIGENVALUES 
2.5497E+05 


NEW SHIFT * 5.2415E*05 MODES STILL NEEDED = u 

0**» USER INFORMATION MESSAGE 4158 STATISTICS FOR SYMMETRIC DECOMPOSITION OF DATA BLOCK SCRATCH FOLLOW 

NUMBER OF NEGATIVE TERMS ON FACTOR DIAGONAL = 1 

0*«* USER INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION 

TRIAL EIGENVALUE * 5 . 2415U7E*05, CYCLES s 1.152254E*U2 NUMBER OF EIGENVALUES BELOW THIS VALUE * 

END OF LANCZOS RUN 

WARNING FLAG = 0 

NO. OF MODES COMPUTED * 1 


COMPUTED MODES 

2 . 54 9697 050 2 50 3 BE* 05 
1 TPS RESULTS ANALYSIS 

0 

0 


MARCH 2 9. 1993 MSC/NASTRAN 7/29/92 PAGE 10 


EIGENVALUE ANALYSIS 


SUMMARY (LANCZOS ITERATION) 


1 TPS RESULTS ANALYSIS 
0 


BLOCK SIZE USED 2 

NUMBER OF DECOMPOSITIONS 2 

NUMBER OF ROOTS FOUND 1 

NUMBER OF SOLVES REQUIRED 6 

TERMINATION MESSAGE : REQUIRED NUMBER OF EIGENVALUES FOUND. 

MARCH 29, 1993 MSC/NASTRAN 7/29/92 PAGE 11 


MODE EXTRACTION EIGENVALUE 


REAL EIGENVALUES 
RADIANS CYCLES 


GENERALIZED GENERALIZED 


31 



HO. ORDER 

1 X 

TPS RESULTS ANALYSIS 


2 549497E*o5 


5 . 049452E*02 


4. 034453E»01 


NASS 

1 . 0000001*00 


STIFFNESS 
2 . S49497E»05 


MARCH 29, 1993 MSC/NASTRAN 7/29/92 PACE 


' USER INFORMATION MESSAGE 5222 .UNCOUPLED SOLUTION ALGORITHM USED. 
TPS RESULTS ANALYSIS 


MARCH 29, 199J MSC/NASTRAN 7/29/92 PACE 13 


‘“THIS IS NODE I 
*** (09*22-92) KNT» 

‘"‘NO. OF COLUMNS « 

1 TPS RESULTS ANALYSIS 


MARCH 29. 1993 MSC/NASTRAN 7/29/92 


"OSTR2 (09-14-92) a 

TPS RESULTS ANALYSIS 


MARCH 29, 2993 MSC/NASTRAN 7/29/92 


XY-OUTPUT SUMMARY (AUTO OR PSDF| 

RMS NO. POSITIVE XMIH FOR XMAX FOR YNIN FOR X FOR 

VALUE CROSSINGS ALL DATA ALL DATA ALL DATA YNIN 

4.439254E-01 4.0247J0E*01 2.0001-01 1.500E*03 2.947E-10 1.500E.03 2 523E- 

4.4392S4E-01 


0 PLOT CURVE FRAME 

TYPE TYPE NO CURVE ID 
0 PSX>F DISP 1 41 ( 5) 

RHSDIS1 ■ 4 . 4J9254E-01£5a 

NOOE5 • 410 

‘“HOOENO CONI NOD EDI Sa 41 243 

‘FREQUENCY a 4.024730E*01 

"RMSDI SI a 4.439254E-01 

‘TOE RMS DIS. AT POINT 42 IS 

LOCALNAJCa 4 . 4392541-02 LOCALNINa 

‘“THE MAX R MS DIS. IS 4.639254E-QI 

0* * * USER INFORMATION MESSAGE *110 (OUTPX2! END -OF -DATA SIMULATION ON FORTRAN UNIT 12 


YMAX POR X FOR* 
ALL DATA YMAX 

4.0lUE*ai 


1 . 0000001*00 


4.439254E-01 WITH TOTAL OF 
1 . 640173S-01 


TPS RESULTS ANALYSIS 


I MAXI HUM SHE OF FORTRAN RECORDS WRITTEN . 
(MUMPER OF FORTRAN RECORDS WRITTEN a 
(TOTAL DATA WRITTEN POR EOF MARKER * 


1 WORDS . ) 

1 RECORDS.) 

I WORDS. ) 

MARCH 29, 1993 MSC/NASTRAN 


LANC20S PARAMETER VERIFICATION 


INITIAL PROBLEM SPECIFICATION 
DECREES OF FREEDOM a 
NUMBER OF MOOES 
MODS FLAG a 

PROBLEM TYPE a 


SHIFTING SCALE 


4. 257 J E*05 


MESSAGE LEVEL a 
OUTPirr UNIT ■ 
SHE OF WORKSPACE a 
MAXIMUM BLOCK SHE a 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 
ACCURACY REQUIRED 


J . 948E-Q00 1 
1 . OUOE+2463 
0. 0Q0QE*00 
J.4960E-10 


AFTER PROBLEM SPECIFICATION CHECKING 
NUMBER OF MODES a 1 

PROBLEM TYPE a t 

SHIFTING SCALE . 4.2S73E*0S 


WORKSPACE ALLOCATION 

LANCIOS BLOCK SHE > 
MAX. BLOCK STEPS 
MAX. MODES 


100 

245 


LEFT END POINT 
RIGHT END POINT 
C ENTER FREQUENCY 


MAX. RITt VALUES 
MAX. SO. VECTORS 
WORKSPACE USED 


3 . 94 BE -0001 
2 . OOOE*2443 
O.OOOGE*OQ 


200 

245 

44159 


CP TINE ALLOWED 


MAX. TRUST REGIONS 


1 . 7?»OC*03 


NUMBER OF USER SUPPLIED VECTORS 


NEW SHIFT a 3.947«E-0l MODES STILL NEEDED » 1 

USER INFORMATION MESSAGE 5010. STURM SEQUENCE DATA TOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE a 3 9474421-01, CYCLES * 1.0O0000E-01 NUMBER OF EIGENVALUES BELOW THIS VALUE * 

ACCEPTED EIGENVALUES 
2. 1944 E* 07 


3 . 2574E*07 


NEW SHIFT • 2.73 19E*07 MODES STILL NEEDED a -2 

0**» USER INFORMATION MESSAGE 4 15B- - -STATISTICS FOR SYMMETRIC DECOMPOSITION OF DATA BLOCK SCRATCH FOLLOW 
NUMBER OF NEGATIVE TERMS ON FACTOR DIAGONAL a l 

0*»* USER INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE a 2 73 14991*07 . CYCLES . 4.314443E*02 HUMBER Of El C DEVALUES BELOW THIS VALUE a 

END OF LANCIOS RUM 

WARNING FLAG . 0 

NO. OF MODES COMPUTED a 1 

COMPUTED MODES 

2. 194395402 17577 E* 07 
1 TPS RESULTS ANALYSIS 


MARCH 29, 1993 MSC/MASTRAM 7/29/92 PAGE 17 


EIGENVALUE ANALYSIS SUMMARY 


( LANCIOS ITERATION) 


BLOCK SHE USED 2 

NUMBER OF DECOMPOSITIONS j 

NUMBER OF ROOTS FOUND 1 

NUMBER OF SOLVES REQUIRED « 


1 TPS RESULTS ANALYSIS 
0 


TERMINATION MESSAGE t REQUIRED NUMBER OF EIGENVALUES POUND. 

MARCH 29, 1993 MSC/MASTRAM 7/29/92 PAGE IB 


MOOE EXTRACTION EIGENVALUE 

NO. ORDER 

1 X 2.194396B+07 

1 TPS RESULTS ANALYSIS 


REAL E I C 
RADIANS 


KNVALUES 

CYCLES 


GENERALIZED 


4 . 4B4572E*03 


7.4519111*02 

MARCH 


NAM 

1 . (KKKKK>l»00 

29. 1993 KK/NASTRAM 


mmtsa 

2.1943941*07 
7/29/92 PAGE 


19 


32 


0 

<>••• USE* INFORMATION MESSAGE SHI .UNCOUPLED SOLUTION ALGOfU THN USED. 

1 TPS RESULTS ANALYSIS MARCH 15. 1953 MSC/NASTRAN 7/29/92 PACE 20 


‘THIS IS MODE • 

" (09-22-92) KMT* 

‘NO. OF COLUMNS. 

TPS RESULTS ANALYSIS 


MARCH 29, 1993 MSC/NASTRAN 7/29/92 PACE 


1 ‘OSTR2 (09-14-92)* 

TPS RESULTS ANALYSIS 


29. 1993 MSC/NASTRAN 7/29/92 PACE 


X Y - OUTPUT SUMMARY 
RMS NO. POSITIVE 

VALUE CROSSINGS 

1.723170E-02 7.445943E*02 

1 . 723 170E-02 


(AUTO OR 
XM2N FOR XI 
ALL DATA 
1 . OOOE-Ol 


ALL DATA 
. S00E+03 


) PLOT CURVE FRAME 
TYPE TYPE MO. CURVE ID. 

I PSDF DISP 2 411 3 

“‘RKSD1S1* 1.723170E-021S. 

‘NODES* 410 

“'NOOENO CONI NODEDIS* 41 243 

“-FREQUENCY. 7.445943E.02 

‘RKSDIS1 • 1 . 7S1962E-02 

*“THI RMS DIS. AT POINT 41 16 

“‘LOCALKAX* 1.7519421-02 LOCALM1N. 

"“THE MAX RMS DIS. IS 1.751942E-02 

*** USER INFORMATION MESSAGE 4110 (OUTPX2 1 O2D-0F-DATA SIMULATION ON FORTRAN UNIT 12 


PSDF) 

IX FOR YKIN FOR 
ALL DATA 
2 . 735E-10 


X FOR 
YNIN 

1 . 500E*0 3 


YKAX FOR X FOR* 
ALL DATA YKAX 
2 . S24E-OS 7 . 46 1 E*02 


.•33437E-01 


. 7231 7OE-02 WITH TOTAL OF 
. 410041 E-03 


“ I TDUTI OH NO . ( - 1 * LI N EAR I * 
•‘UGNIMAX* 3.31 942 9E-01 

' MOGMIHAX* 1.731942E-U2 

‘ABSOLUTE NORM . 1 794B06E.01 

TPS RESULTS ANALYSIS 


(MAXIMUM SIZE OF FORTRAN RECORDS WRITTEN 
(NUMBER OF FORTRAN RECORDS WRITTEN 
(TOTAL DATA WRITTEN FOR EOF MARXES 
0 XNORM. 3. 14 73 1 IE-01 


1 WORDS. | 

1 RECORDS. ) 
l WORDS . ) 


MARCH 29. 1993 MSC/NASTSAM 7/29/92 


LANCZOS PARAMETER VERIFICATION 


INITIAL PROBLEM SPECIFICATION 
DECREES OF FREEDOM * 
NUMBER OF NODES 
MODE FLAG 

PROBLEM TYPE • 


SHIFTING SCALE 


4.24791*05 


MESSAGE LEVEL 
OUTPUT UNIT « 

SIZE OF WORKSPACE • 
MAXIMUM BLOCK SIZE * 


LEFT END POINT 
RIGHT DID POINT 
CENTER FREQUENCY 
ACCURACY REQUIRED 


3 . 34 BE- 0001 
1 . QOO£*2463 
0. OOOOE.UO 
3 . 69601- 10 


AFTER PROBLEM SPECIFICATION CHECKING 
NUMBER OP MODES * 1 

PROBLEM TYPE « 1 

SHIFTING SCALE ■ 4.2479K.03 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 


3 . 94 BE-0001 
l.VOOE*2463 
0 . 0000E*00 


CP TIME ALLOWED 


1.7660E.01 


WORKSPACE ALLOCATION 

LAMCZOS BLOCK SIZE • 
MAX. BLOCK STEPS • 
MAX. MOOES • 


100 

245 


MAX RJTZ VALUES 
MAX. SO VECTORS 
WORKSPACE USED 


200 

245 

44159 


MAX. TRUST REGIONS 


OF USER SUPPLIED VICTORS i 


NEW SHIFT • 3. 947*1-01 MODES STILL NEEDED • 1 

USER INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE ■ 3 . 94744 2E-01 , CYCLES * 1.00Q0001-PI NUMBER OF EIGENVALUES BELOW THIS VALUE < 


AC CEP TE D EIGENVALUES 
9 . 3066E+06 

HEW SHIFT • 1. 1172 E* 07 NODES STILL NEEDED * 0 

0*** USER INFORMATION MESSAGE 4 154- - -STATISTICS FOR SYMMETRIC DECOMPOSITION OF DATA BLOCK SCRATCH FOLLOW 
NUMBER OF NEGATIVE TERMS ON FACTOR DIAGONAL > 1 

0*** USER INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE • 1.101*41*07. CYCLES ■ 3.31966*1*02 NUMBER OF EIGENVALUES BELOW THIS VALUE > 1 


END OF LANCZOS RUN 

WARMING FLAG • 0 

NO. OF MODES COMPUTED - 1 

COMPUTED MODES 

9 306344249076901*04 

1 TPS RESULTS ANALYSIS MARCH 29. 1993 MSC/NAfUUN 7/29/93 PACE 24 

0 

0 


EIGENVALUE ANALYSIS SUMMARY (LANCZOS ITERATION I 


BLOCK SIZE USED 2 

NUMBER OF DECOMPOSI TI (4IS 2 

NUMBER OF ROOTS FOUND I 

NUMBER OF SOLVES REQUIRED • 

TERMINATION MESSAGE : REQUIRED NUMBER OF EIGENVALUES POUND 

1 TPS RESULTS ANALYSIS MARCH 2*. 1*93 MSC/NASTRAN 7/29/92 PACE 35 

0 

REAL EIGENVALUES 

MOM EXTRACTION S3 C DIVALUE RADIANS CYCLES 

MO. ORDER 

1 1 9. 3045*61*06 3.0S0670E+03 4. *532921*02 

Z IRS RESULTS ANALYSIS MARCH 


GENERALIZED GENERALI ZED 

MASS STIFFNESS 

1 . 0000001*00 9 . 3065461*06 

29. 1993 MSC/NASTRAN 7/29/92 PACE 24 
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<>•** USER INFORMATION MESSACE 5222 .UNCOUPLED SOLDTIOH ALCOR1 THN USED. 

1 TPS RESULTS ANALYSIS 


KAROt 29, 1992 MSC/NASTRAN 7/29/92 PACE 27 


"THIS IS NOOE t 
“‘(09-22-92) KMT. 

“NO. OP COLUMNS* 

TPS RESULTS ANALYSIS 


KAACH 29, 1993 MSC/NASTRAN 7/2 9/92 PACE 2* 


* OSTR2 (09-14-921* 

TPS RESULTS ANALYSIS 


MAACH 29, 199J MSC/NASTRAN 7/29/92 PACE 29 


OUTPUT SUMMARY 
RMS NO. POSITIVE 

VALUE CROSSINGS 

3.744705E-O2 4 . 8501051*02 

3.744705E-02 


(AUTO OR P S D F ) 

XMIN FOR XMAX FOR YMIN FOR 

ALL DATA ALL DATA 

. S00E*0J 2. 51 BE- 10 


ALL DATA 
1. OOOE-Ol 


> X Y 

) PLOT CURVE FRAME 

TYPE TYPE NO. CURVE ID. 

> PSDP DISP 3 «t( 5) 

‘“RMSDISl- 3,7447051-0225. 

* NODES* 410 

“‘NODENO CONI HOC EDI £• 41 

““‘FREQUENCY. 4.B50105E*Q2 

“‘RMSDISl. J. 7447051-02 

“‘THE RMS DIS. AT POINT 41 IS 

“' LOCALKAXa 1.1439421-02 LOCALMIM* 

‘“THE MAX RMS DIS. IS 3.744705E-02 
>•'• user INFORMATION MESSAGE 4110 ( OOTPX2 ) Off) -OF- DATA SIMULATION ON FORTRAN UNIT 12 


X FOR 
YMIN 

1 . SOOE.OJ 


YKAX FOR X FOR* 
ALL DATA YMAJt 
l . 7B6E-04 4.851E*02 


243 -1 . OOQOUOE*00 


3.7447051-02 WITH TOTAL OF 
3 . 744705E-02 


ITERATION NO. (-1. LINEAR! • 
“‘UCNINAJt* 1.74597JE-01 

* “MUCN1 MAX. J.744705E-02 

“"ABSOLUTE NORM . 3 . 6625 1 1E.00 

TPS RESULTS ANALYSIS 


(MAXIMUM Sill OF FORTRAN RECORDS WRITTEN i 
(NUMBS* OF FORTRAN RECORDS WRITTEN . 
(TOTAL DATA WRITTEN FOR EOF MARKER i 
i XNORM* 1 . 371502E-01 


1 WORDS. I 
1 RECORDS. ) 


MARCH 29, 1993 MSC/NASTRAN 7/29/92 PACE 


LAMCIOS PARAMETER VERIFICATION 

INITIAL PROBLEM SPEC I FI CATION 
DECREES or FREEDOM ■ 
HUMBER OF MODES a 

NODE FLAG . 

PROBLEM TYPE • 


SHIFTING SCALE 


3 . 71B4E*05 


MESSAGE LEVEL 
OUTPtrr UNIT 
SHE OF WORKSPACE 
MAXIMUM BLOCK SHE 


LEFT DID POINT • 1.940E-UOO1 

RIGHT END POINT a 1.0001*2463 

CENTER FREQUENCY a O.OOOOE*UO 

ACCURACY REQUIRE) . J.696QE-1U 


AFTER PROBLEM SPECIFICATION CHECKING 
NUMBER OF MODES . 1 

PROBLEM TYPE . 1 

SHIFTING SCALE a 3.71B4E*05 


WORKSPACE ALLOCATION 

LANCZOS BLOCK SHE » 
MAX. BLOCK STEPS * 
MAX. MODES ■ 


100 

245 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 


MAX. AIT2 VALUES 
MAX. S.O. VECTORS 
WORKSPACE USED 


3 . 948E-0001 
1 OOOE+246J 
0 . 0000£«00 


200 

245 

44159 


CP TIME ALLOWED a 1 . 7550E*03 


X. TRUST REGIONS < 


NUMBER OP USER SUPPLIED VECTORS a g 

NEW SHIFT > 3.947BE-01 MOOES STILL HEEDED . 1 

0*** USSR INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE a 3 . 947842E-01 , CYCLES * 1 0000001-01 HUMBER OF EIGENVALUES BELOW THIS VALUE a 

ACCEPTED EIGENVALUES 
5 . 0OB2E-O6 

NEW SHIFT • 5 96 1 SE* 06 MODES STILL NEEDED • 0 

0*** USER INFORMATION MESSAGE 4 IS*--- STATISTICS FOR SYMMETRIC DECOMPOSITION OF DATA BLOCK SCRATCH FOLLOW 
NUMBER OF NEGATIVE TERMS ON FACTOR DIAGONAL . 1 

0*»* USER INFORMATION MESSAGE 5010, STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION 

TRIAL EICXNVALUt a 5. ><1B34E»04. CYCLES . 3.BB6066C«92 HUMBER OF EIGENVALUES BELOW THIS VALUE « 

END OF LANCZOS RUN 

WARNING FLAG a 0 

MO. OF NODES COMPUTED • 1 


COMPUTED MODES 

5 . 00# 1 7975026906 E* 06 
TPS RESULTS ANALYSIS 


MARCH 29. 1993 MSC/NASTRAN 


EIGENVALUE AMALYSIS SU 


(LANCZOS ITERATION} 


1 TM RESULTS ANALYSIS 
0 


•LOCK SIZE USED 2 

NUMBER OF DECOMPOSITIONS 2 

NUMBER OF ROOTS FOUND 1 

HUMBER OF SOLVES REQUIRED « 

TERMINATION MESSAGE i REQUIRED NUMBER OF EIGENVALUES FOUW>. 

MARCH 19, 1993 MSC/NASTRAN 7/29/92 PAGE 3 2 


NODE EXTRACTION EIGENVALUE 

NO. ORDER 

I 1 S.00«1B0E*06 

l TPS RESULTS AMALYSIS 


REAL EIGENVALUES 
RADIANS CYCLES 

2 . 237#96E*03 J. 5417231*02 

MARCH 


GENERALIZED 


1 . 0000001*00 

29. 1993 MSC/NASTRAN 


GENERALIZED 
STIFFNESS 
5 . U0B1R0E*06 
7/29/92 PAGE 


33 


34 


<)••• DEER IMFORMATIOM MESSAGE 5222 , UNCOUPLED SOLUTIOK ALGORITHM USED. 
1 TPS RESULTS ANALYSIS 


MARCH 29, 1993 MSC/NASTRAM 7/29/92 PACE 24 


“•-THIS IS MODE * 

(03-22-92) KMT. 

NO. OP COLUMNS- 

1 TPS RESULTS ANALYSTS 


MARCH 29. 1993 MSC/NASTRAM 7/29/92 PACE 35 


***OSTR2 (09-14-92) • 

1 TPS RESULTS ANALYSIS 


MARCH 29, 1993 MSC/NASTRAM 7/29/92 PACE 34 


ALL DATA 

l.OOOE-01 


► XY-OUTPBT 

I PLOT CURVE FRAME RMS 

TYPE TYPE NO. CURVE ID. 

I PSDF DISP 4 41 ( 51 

“‘RMSDIS1. 7.0S43J9E-Q21S. 

***NOOES« 410 

-'•NOOEHO CONI MODEDIS* 4 

*• * FREQUENCY* 3.5589831*02 

•'•RMSDIfil* 7 . 05433 BE- 02 

THE RMS DIS. AT POINT 41 IS 

* LOCALMAX* 1.S49021E-02 LOCALMIH* 

'"THE MAX RMS DIS. IS 7.0S4J38E-02 

**• 0*** IMTORMATIOM MESSAGE 4110 ( OUTPX 2 ) END -OF - DATA SIMULATION ON FORTRAN UNIT 12 


SUMMARY 
RMS NO. POSITIVE 

VALUE CROSSINGS 

7.0543381-02 3 . 5583831*02 

7.0543381-02 

243 -1 0000001*00 


7.054338E-02 WITH TOTAL OF 
7 . 054338E-02 


(AUTO OR PSDF) 

XMIN FOR XMAX FOR YMIN FOR X FOR 
ALL DATA ALL DATA YMIN 
1 . 500E.0J 3.175E-10 1 . 50UE+03 


YMAX FOR X FOR* 
ALL DATA YMAX 
8. 8B4I-04 3 . 54 11*02 


ITERATION MO. ( -1 .LINEAR) ■ 

■■‘UCNIMAX. 1 . 0402 12 E- 01 

* * “ MUCNI MAX. 7. 05433 8E-02 

"•ARSOLUTE NORM . 5.029357E-01 

TPS RESULTS ANALYSIS 


(MAXIMUM SUE OF FORTRAN RECORDS WRITTEN * 
INUMBER OF FORTRAN RECORDS WRITTEN . 
(TOTAL DATA WRJTTDi FOR EOF MARKER • 
2 XMORM- 3.5478791-02 


1 WORDS.) 

1 RECORDS.) 
I WORDS.) 


MARCH 29. 1993 MSC/NASTRAM 7/29/92 PACE 


LAWCIOS PARAMETER VERIFICATION 


INITIAL PR08LEM SPECIFICATION 
DECREES OF FREEDOM . 
NUMBER OF MOOES • 

MODE FLAG • 

PROBLEM TYPE . 

SHIFTING SCALE . 3. 


MESSAGE LEVEL 
OUTPUT UNIT 
SUE OF WORKSPACE < 
MAXIMUM BLOCK SUE . 


6230E-.0S 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 
ACCURACY REQUIRED 


. 3.94RE-0001 
> I . OUUE* 24 63 
■ O.OOOOE*UO 
. 3 . 69C0E- 1 0 


AFTER PROBLEM SPECIFICATION CHECKING 
NUMBER OF MODES ■ 1 

PROBLEM TYPE - l 

SHIFTING SCALE * 3.6230E*05 

WORKSPACE ALLOCATION 

LANCZOS BLOCK SUE > 2 

MAX. BLOCK STEPS - 100 

MAX. MODES • 245 

HUMBER OF USER SUPPLIED VECTORS ■ 


LEFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 


MAX RITt VALUES 
MAX. S.O. VECTORS 
WORKSPACE USED 


3 . 9481-0001 
1 . 000E.24O3 
0 . 0000E*QQ 


200 

245 

44159 


CP TIME ALLOWED 


MAX. TRUST REGIONS > 


1.7«30E*0J 


3. 94781-01 


MODES STILL NEEDED * 


0*** USER INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE • 3 . 947S42E-01 . CYCLES • 1 . OOOOOOE-Qt NUMBER OF EIGENVALUES BELOW THIS VALUE ■ 

ACCEPTED EIGENVALUES 
3.7«61E*06 

NEW SHIFT . 4 . 58B6E* 06 MOOES STILL HEEDED » 0 

0*»* USER INFORMATION MESSAGE 41 58- --STATISTICS FOR SYMMETRIC DECOMPOSITION OF DATA BLOCK SCRATCH FOLD 
NUMBER OF NEGATIVE TERNS ON FACTOR DIAGONAL • 1 

©**• USSR INFORMATION MESSAGE 5010, STURM SEQUENCE DATA FOR EIGENVALUE EXTRACT! ON. 

TRIAL EIGENVALUE • 4 . 58843 01*06, CYCLES * 3.4093721*02 NUMBER OF EIGENVALUES BELOW INIS VALUE . 

END OF LANCZOS RUN 

WARNING FLAG • 0 

NO. OF MODES COMPUTED . 1 


3.766129250286801*06 
TPS RESULTS ANALYSIS 


MARCH 29. 1993 MSC/NASTRAM 7/29/92 PACE 38 


eigenvalue analys 


S SUMMARY 


(LANCZOS ITERATION I 


BLOCK SUE USED 2 

NUMBER OF DECOMPOSITIONS 2 

NUMBER OF ROOTS FOUND I 

NUMBER OF SOLVES REQUIRED S 


1 TPS RESULTS ANALYSIS 


TERMINATION MESSAGE i REQUIRED NUMBER OF EIGENVALUES FOUND. 

MARCH 29. 1993 MSC/NASTRAM 7/39/92 PACE 39 


NODE EXTRACTION EIGENVALUE 

NO. ORDER 

1 1 3 . 744129 E* 04 

1 TPS RESULTS ANALYSIS 


REAL EIGENVALUES 
RADIANS CYCLES 


1.940452t*03 J.088643E*U2 

MARCH 


GENERALIZED 

MASS 

1 . 0000001*00 

29. 1993 NSC/MASTRAN 


GENERALIZED 

STIFFNESS 

3.7C«129E*U« 

7/39/92 PACE 40 
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o* * • USER INFORMATION MESSAGE 5222 UNCOUPLED SOLUTION ALGORITHM USED. 
1 TPS RESULTS ANALYSIS 


MARCH 2 9. 1993 MSC/KASTRAN 1/29/92 PACE 41 


'THIS IS MODE • 

<09-22*92) KKT* 

“"NO. OP COLUMNS. 

1 TPS RESULTS ANALYSIS 


MARCH 29. 1993 MSC/HASTRAN 7/29/92 PAGE 42 


"'051312 I 09- 1 4*92 > a 
1 TPS RESULTS ANALYSIS 


MARCH 29, 1993 MSC/KASTRAN 7/29/92 PACE 43 


XY-OUTPUT SUMMARY 
RMS NO. POSITIVE 

VALUE CROSSINGS 

9 . 195 J 19E-02 3.0B4597E.02 

9.1953191-02 


AUTO OR PSDF) 

ENIN FOR XMAJt FOR YMIM FOR 

ALL DATA ALL DATA 
1.500E+03 J.414E-10 


ALL DATA 
1 . Q00E-01 


) PLOT CURVE FRAME 
TYPE TYPE NO. CURVE ID. 

> PSDF DISP 5 4W 5) 

‘"RNSDISl- 9. 19SJ19E-02Z5- 

'"MOOE5- 410 

' "NOODW CONI NOCEDIS- 4 

■"FREQUENCY. 3.0i«591E*02 

'"RNSDIS1- 9 . 195J19E-02 

‘"THE RMS DIS. AT POINT 41 IS 

"'LOCALMAX- 9 . 19531 9E-02 LOCALNIN. 

"‘THE MAX RNS DIS. IS 9.195319E-02 

*** USER INFORMATION MESSAGE 4110 (0UTPX2 | KND-OF-DATA SIMULATION ON FORTRAN UNIT 12 

(MAXIMUM SIZE OF FORTRAN RECORDS WRITTEN - 1 WORDS. ) 

(NUMBER OF FORTRAN RECORDS WRITTEN - 1 RECORDS.) 

(TOTAL DATA WRITTEN FOR EOF MARKER . 1 WORDS . ) 

*" ITERATION NO. (-1- LINEAR)- 3XNORM- J.47O41SE-0J 

'"UGNIMAX- • . 024274E-O2 

'"MUCN1MAX- 9. 195319E-02 

ABSOLUTE NORN - -3 . 99 14 1 3 E- 02 

TPS RESULTS ANALYSIS 


X FOR 
YMIM 

l . 500E+03 


YMAX FOR X FOR* 
ALL DATA YMAX 
1 . 701E-03 3 . 091E*02 


243 l.OOOOOOE-OO 


9. 1953 19E-02 WITH TOTAL OF 
1.972773E-02 


MARCH 29, 1993 MSC/KASTRAN 7/29/92 PACE 44 


LANCZ06 PARAMETER VERIFICATION 


INITIAL PROBLEM SPECIFICATION 
DECREES OF FREEDOM a 
NUMBER OF MOOES 
MOCE FLAG • 

PROBLEM TYPE 


SHIFTING SCALE 


3 . 4252E*05 


MESSAGE LEVEL 
OUTPUT UNIT 
SIZE OF WORKSPACE . 
MAXIMUM BLOCK SIZE < 


LEFT Off) POINT 
RIGHT END POINT 
CENTER FREQUENCY 
ACCURACY REQUIRED 


3 . 94ME-0001 
1 . 000E.2 46 ) 
0. OQOOE.OO 
3 . 6960E-1G 


AFTER PROBLEM SPECIFICATION CHECKING 
NUMBER OF NODES ■ 1 

PROBLEM TYPE . I 

SHIFTING SCALE . 3.42S2E*05 


WORKSPACE ALLOCATION 

LMCZOC BLOCK SIZE 
MAX. BLOCK STEPS 
MAX. MOOES 


100 

245 


LIFT END POINT 
RIGHT END POINT 
CENTER FREQUENCY 


MAX. RITZ VALUES 
MAX. SO. VECTORS 
WORKSPACE USED 


3 . 946E-0001 
1.0Q0C*244J 
0. 0000 E* 00 


200 

245 

44159 


CP TIME ALLOWED • 1.731OE*03 


MAX. TRUST REGIONS i 


NUMBER OF USER SUPPLIED VECTORS 


3. 947BE-Q1 


MOOES STILL NEEDED - 


0*** USER INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACT! ON . 

TRIAL EIGENVALUE • 3 . 947B42E-01 , CYCLES • 1. 0000001-01 NUMBER OF EIGENVALUES BELOW THIS VALUE « 

ACCEPTED EIGENVALUES 

J.HUttOC 

NEW SHIFT • 4 . 5465E-04 NODES STILL NEEDED • 0 

0*** USE* INFORMATION MESSAGE 415B-- -STATISTICS FOR SYMMETRIC DECOMPOSITION OF DATA BLOCK SCRATCH FOLLOW 
MUNBIR OF NEGATIVE TERNS ON FACTOR DIAGONAL • 1 

0*** UM* INFORMATION MESSAGE 5010. STURM SEQUENCE DATA FOR EIGENVALUE EXTRACTION. 

TRIAL EIGENVALUE * 4 . 5444741*04. CYCLES « 3.4010321-02 NUMBER OF EIGENVALUES BELOW THIS VALUE ■ 

END OF LANCZOS RUN 

MARRING FLAG . 0 

NO. OF NODES COMPUTED ■ l 


COMPUTED NODES 

3 . 741?4B539B4742E*04 
TPS RESULTS ANALYSIS 


MARCH 29. 1993 MSC/KASTRAN 7/29/92 PACE 


EIGENVALUE ANALYSIS SUMMARY (LANCZOS ITERATION) 


BLOCK SIZE USED 2 

NUMBSR OF DECOMPOSITIONS 2 

NUMBER OF ROOTS FOUND 1 

NUMBIR OF SOLVES REQUIRED ft 


1 TPS RESULTS ANALYSIS 

0 


TERMINATION MESSAGE t REQUIRED NUMBER OF EIGENVALUES FOUND. 

MARCH 29. 1993 MSC/KASTRAN 7/29/92 PAGE 4« 


NODE EXTRACTION EIGENVALUE 

NO. ORDER 

1 1 3.7*17491.04 

t TPS RESULTS ANALYSIS 


REAL EIGENVALUES 
RADIANS 


1. 93952*1*03 


J . 0444551-02 

MARCH 


NASS 

1. 0000001*00 

29. 1993 NSC/NASTRAM 


STIFFNESS 

3.7417491*04 

7/29/92 PACE 


47 
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MARCH 29, 1993 MSC/NASTRAN 


7/29/92 PACE 48 


0 

<>•** USER INFORMATION MESSAGE 5222 .UNCOUPLED SOUTH ON ALGORITHM USED. 
1 TVS RESULTS ANALYSIS 


““THIS IS MOCK • i 

*“* <09-22-92] KXT» 4LMODES* 

““NO. OF COLUMNS. I 

1 TVS RESULTS ANALYSIS 


0 

““OSTR2 (09-14-92) . -I 

1 TVS RESULTS ANALYSIS 


MARCH 29. 1993 MSC/NASTRAN 7/29/92 PAGE 49 


MARCH 29. 3993 MSC/NASTRAN 7/29/92 PAGE 50 


0 

0 XY-OUTPUT SUMMARY | A U T 0 OR P S D F ) 

0 PLOT CURVE FRAME RMS NO POSITIVE XMIN FOR XMAX FOR YKIN FOR 

TYPE TYPE NO. CURVE ID VALUE CROSSINGS ALL DATA ALL DATA ALL DATA 

0 PSDF DISP t til 5) B.917840E-02 J.084409E.02 1.000E-01 1.500E.0J 3.207E-1Q 

““RMSDrSl. • . 9170COE-Q225. B.917B40E-02 

“‘NODES- 410 

““NODINO CONI NOOEDIS. 41 243 -1 . 000G00E*O0 

‘“FREQUENCY. 3.084809E.02 

‘“RKSD1S1. 8 . 91 7860E-02 

““THE RMS DIS. AT POINT II JS 8.917840E-02 WITH TOTAL OF 1 NODES 

‘“LOCALMAX. 1 . 954 13SE-02 LOCALMIH. 8.91704QE-O2 

““THE MAX RMS DIS. IS B.917I40E-Q2 


0*«. USER INFORMATION MESSAGE 4110 (OUTPX2 ) END-OF-DATA SIMULATION ON FORTRAN UNIT 12 

(MAXIMUM SUE OF FORTRAN RECORDS WRITTEN ■ l WORDS. > 

I MUM ICR OF FORTRAN RECORDS WRITTEN • 1 RECORDS. ] 

{TOTAL DATA WRITTEN FOR EOF MARKER « 1 WORDS . ) 

ITERATION NO. { -1-L1 NEAR) • 4 XMORM. 1.4245441-03 

* “UCNIMAX- 9.01 1798E-02 

““MJGNIKAX. 8.917840E-02 

““ARSOUTTE NORM • 1.053379E-O2 

0 ... USER INFORMATION MESSAGE *109 (OUTPX2) 


THE LABEL IS PDALABEL FOR FORTRAN UNIT 12 


(MAXIMUM SUE OF FORTRAN RECORDS WRITTEN 

(NUMBER OF FORTRAN RECORDS WRITTEN 

(TOTAL DATA WRITTEN FOR TAPE LABEL 
1 USER INFORMATION MESSAGE 4114 (0UTPX2) 

DATA BLOCK OUCV1PAT WRITTEN ON FORTRAN UNIT 12. TRL ■ 

101 0 448 0 

I MAXI HUN POSSIBLE FORTRAN RECORD SUE 
(MAXIMUM SUE OF FORTRAN RECORDS WRITTEN 

(NUMBER Of FORTRAN RECORDS WRITTEN 

(TOTAL DATA WRITTEN FOR DATA BLOCK 

USER INFORMATION MESSAGE *110 (OUTPX2] EMD-OF-DATA SIMULATION ON FORTRAN UNIT 12 

I MAXIMUM SUE OF FORTRAN RECORDS WRITTEN • 1 WORDS.! 

(NUMBER OF FORTRAN RECORDS WRITTEN ■ 1 RECORDS. | 

(TOTAL DATA WRITTEN FOR EOF MARKER . 1 WORDS . | 

* * * END OF JOB • • • 


7 WORDS. | 

I RECORDS. ) 
17 WORDS. ) 


24578 WORDS. ) 
<48 WORDS. ) 

24 RECORDS] 
• 29 WORDS.] 


X FOR YMAX FOR X FOR* 
YMIN ALL DATA YMAX 
1. 5001. 03 1 . 5261-03 3 . 091E*Q2 


0 


0 


Overall Rms Displacements 

When the overall rms displacements are requested, DISP=ALL is required in the Control 
Deck. The overall rms displacement is formed by multiplying the maximum rms displacement 
by the updated mode shape. The overall rms displacement can be found in the *.f06 output 
file. It also can be extracted from an OUTPUT2 file by including a “PARAM.POST,- 1” card 
in either the Case Control or Bulk Data Decks. 

5.4 Static and Nonlinear Random Analysis 

The modified MSC/NASTRAN SEMELRR solution sequence can be combined with 
other solution sequences to handle the effect of static mechanical and thermal loads on the 
nonlinear random response. An example is shown for a hexagonal thermal protection system 
panel subjected to combined thermal and random acoustic loads. The deformed shape due 
to the static thermal load is first obtained using the SOL 101 procedure. A RESTART of the 
SEMELRR solution sequence is then run to obtain the dynamic response due to the combined 
load. The results shown were obtained by post-processing the output with PDA Patran [20], 
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SOL 101 data cards 


3 Ttva diMUion o( tha inotropic plata i« 15 * 12 * 0.04 

$ For tha applied Ova raal load, run tr804-l.dat rirat. 

3 thw MPtMTon tr«0«-l.d*t 

$ IN own* TO err the combined thermal and acoustic loadings 

3 FUa ium t P101.dat 

s BUCK DAT* la In tha PLATE. RDF 

ID NSC, CH2ANC 
TINS 100 J 
SOL 101 3 
CXMD 5 
3 

TITL1 - TPS RESULTS ANALYSIS 
SUPER a ALL 
METHOD a 4 
LOAD a 200 

M UP ARM a 10 

3 

ECHO* NONE 3 
SK>S( i 
BEGIN BULK S 
INCLUDE PLATT- 8DF 

$ 2 3 i * « 7 • » °- 

NLPARM 10 1 AUTO 1 PW YES 

$ 2 3 * 5 < 7 * * 0 

P LOAD 2 200 1. 1 "TWO 44 

E1GRL 4 0.01 0 

S 2 3 « 5 4 T 8 9 0 

PARAM COUPMASS1 
PARAH NMLOOP I 

PARAM X4ROT 100. 

PARAM LMODES 1 

PARAM NOCOMPS - 1 

PARAM AUTOSPC YES 

EMDOATA 


In order to run SEMELRR with the data base stored from SOL 101, the following 
command is needed: 

nastran Plll.dat dbs=P101 


SEMELRR data cards 


RESTART VERSIONa l, KEEP $ 

S 

3 Flla naaa : Plll.dat 

I 

S THIS IS THE M001FIED NASTRAN RANDOM RESPONSE SOLUTION SEQUENCE 
S 

ACQUIRE KDOL 

S 

ID CHI AMO. NASTRAN 

$ 

COMPILE SEMELRR SOOOUTaUSRSOU OBJOUTa USRORJ NOLIST NONET 
INCLUDE SEMELRR. DNAP 

COMPILE SEDRCVA SOUOUTaOSRSOU OBJOUTa USHOBJ NO LI ST NONET 
INCLUDE SEDRCVR.DMAP 

COMPILE SUPER3 SOOOUT'USRSOU OBJOUTa USROBJ NOUST NOREF 
INCLUDE 80PERJ.DMAP 

COMPILE SUBDMApmSENG SOUIN*MSCSOO NOUST NOREF J 
ALTER 33 3 

EQU3VX EST/ESTWALHAYS 3 
IMDALTER $ 

SOL SEMELRR S 
UN* SEMELRR 3 
TIME 100 $ 
cm> 

TITLE a PLATE RESULTS ANALYSIS 

SET 1 a 1 THRU 44 

DXSP • ALL S 

STRAIN (FIBER) a 1 $ 

METHOD a « » 

SPCa»« 

SDAMP a 102 3 
DLOAD a J01 3 

loads rr . 300 3 
FREQ > 3 3 
RANDOM* 34 3 

3 

OUTPUT (XT OUT) 

XYPLOT DIStP PSDF/ 41(T3> 

3 

BEGIN BULK 3 

3 2- 3 4 3 4 T • 9 0- 
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$ FREQUENCY RESPONSE INPUT 


)0Q 

201 


* 2 -2 

FRIQI 3 0.00 

DLOAD 201 1.0 

RLOAD1 204 200 

LSEQ 500 300 

PLOAD2 100 1. 

S 

TAJ LED 1 10 

♦QR 0.0 1.0 

* 

RAND PS 4» 1 

TAJRND1 «] 

♦in -1.0 0.0 

♦TR2 2001.0 0.0 

$ 

TABDMR1 102 

♦DPI 0.0 0.04 

t 


100 

1 THRU $4 


2000.0 1.0 ENDT 

1 1.0 43 

0.0 8. 4215*52000. 0 8.4215-52000.0 0.0 

1HDT 


2000.0 0.04 ENOT 


PARAM.RMSTRAIN, 1 S IF MS STRAIN IS NEEDED RESTRAIN. 1 
RARAM , MAXI TSR , 3 S HAJt NUMBER OF ITERATION 
PARAM.AJSNORN. 2.01-2 S AJS. NORN FOR CONVERGE TEST 

PARAN, BETA. 0.5 $ SCALE FOR BETTER CONVERGENCE, RANGE FROM 0.0 TO 1.0 
PAKAM , LCD I SP . 1 
$ 


♦ QR 


♦ TR 

♦ TR2 


♦ DPI 


$ 


Hexagonal Panel 

A hexagonal thermal protection system (TPS) panel similar to the cutout shown in figure 
1 was subjected to both thermal and acoustic loads. The structure is composed of an eight- 
ply carbon-carbon TPS panel with built-up substructure. The TPS panel is connected to 
the substructure with seven titanium rods (posts). The substructure has an aluminum core 
sandwiched between an aluminum and a graphite/epoxy face sheet. The dimensions of the 
panel are given in Table 1, and the finite element mesh is shown in figure 9. The finite 
element model is comprised of 804 triangular elements and seven bar elements with a total 
of 622 nodes. 

The boundary conditions imposed on the panel were designed to minimize thermal 
stresses, and are summarized for each component The edges of the TPS panel are constrained 
in the perpendicular and tangential directions. The edges of the substructure are constrained 
in all rotations and translations. The post connections to the TPS panel were modeled as 
pinned joints using MPC Bulk Data cards. The three translations at the top of the posts 
were equivalent to the three translations at the adjoining locations of the TPS panel. The 
connections between the posts and the substructure were also modeled using MPCs. The 
center post connection was modeled as a rigid link, i.e. all three translations and the two 
rotations at the lower-end node of the post were equivalent to the translation and rotations of 
the adjoining node of the substructure. The remaining post connections to the substructure 
were also modeled as pin joints. 

Table 1 Dimensions for hexagonal TPS panel example problem 


Radius 

13.0 in. 

Overall height 

2.5 in. 

Radius to posts 

8.0 in. 

Carbon-carbon thickness 

0.091 in. 

Substructure thickness 

0.375 in. 
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Center post radius 
Outer post radii 


0.1875 in. 
0.125 in. 


A 2000 °F temperature load was applied to the TPS panel and a 200 °F load was applied 
to the posts and substructure. The thermal displacements and stresses were predicted using 
SOL 101, and are plotted in figures 10 and 1 1. The TPS panel results are essentially those of 
a stress-free thermal expansion while the substructure shows a moderate compressive thermal 
stress with little thermal displacement. The equivalent linearization solution sequence was 
restarted using the data base from the static thermal solution with the initial stresses and 
displacements. The rms thermal-acoustic displacements and stresses were predicted for a 
broadband acoustic excitation ol 150 dB uniformly distributed over the carbon-carbon panel. 
These rms displacements and stresses are plotted in figures 12 and 13. The solution sequence 
converged in four iterations with the convergence enhancement parameter BETA set to 0.5 
and the default convergence criteria. 

The level ol nonlinearity in the response is typically measured in several ways. The 
two most common are the ratio ot the equivalent linear fundamental frequency to the linear 
fundamental frequency (frequency ratio) and the ratio of the equivalent linear maximum 
rms displacement to the linear maximum displacement (amplitude ratio). For this particular 
problem, these ratios were 1.19 and 0.414, respectively, and are typical of moderate to 
extreme geometric nonlinearity. 



Figure 9: Finite element model 
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Figure 10: Deformed plot of the thermal displacement 
vector. Displacements are given in inches. 



Figure 11: Thermal stresses in the radial direction e r in psi. 
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Figure 12. Deformed plot ot the root-mean-square thermal-acoustic 
displacement vector. Displacements are given in inches. 



in psi. 
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Summary 


An equivalent linearization solution sequence used to predict the nonlinear random 
response of structures has been incorporated into MSC/NASTRAN version 61x2. A new 
main SUBDMAP, SEMELRR, and a significantly modified MSC/NASTRAN SUBDMAP 
SEDRCVR are compiled with the MSC/NASTRAN delivered library of SUBDMAPs to 
create the new solution sequence. 

The equivalent linear rms displacements, strains or stresses, and frequencies are calcu- 
lated by an iterative solution method. The numerical results obtained were in good agreement 
with existing solutions. The output requests and the iterative solution method are controlled 
by several new user defined PARAMeters. The versatility of the implementation will en- 
able the analyst to determine the nonlinear random responses for complex structures under 
combined loads. 
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Appendix A Main SUBDMAP 


summer SMIL HA t PROM SOL 111 - SUFBRXLSNBNT MODAL FREQUENCY 

* 

t PHASE 1 

TTFI DR. KJJ.K7J, BJJ.USST, GOA, GOAT. KAA. DYNAMICS, KOO, DIOXIN. DEQ If®, AOS. 
CFLS, SILJ. SLT. CM. CASES. BQBXXNS . MATFOOL, Of . LOO . XOO, PC DM , 

DIT. BCPOT* . CSTNB , Wt, NPTB . XYCDBS . LLL , KLL. ULL. 

MAP, MAPS . ITT. pvrrs , KPS. JUS, KSP , ItW, EOT. 

GEOK2S. GXOK3S, BCTS, 1PTB. DOTH. KBLM, KOICT, SLIST. 

OPSCT.VSLBH, BQBXXNX . BCTX, BCPDTX, SXLX.CS0M48 . VCPS .KPT.DIT.CSTM. 
f STATIC LOAD GENERATION 

PJ. PTELIM.PO, PA, YS.QR.POS.PSS.PL. 
t DYNAMICS ONLY I SIMA) 

•OS , I4CG . MPP . MLA . KLL , MM . OOAMA , CMPHA . CMPHO . MLAA1 . 

MA, OQAQ, KLAA, KLAA. MAA. MAR, MSA, BAA . MAA , MGC, 

S PHASS 2 

PKA . LAMA , CMDD, CX2 DO . CB2 DO , PHDH , IC4DD , VPHPS , 

ULP . ISTML . UOMI , ISTL , 

PHPA, LAMAP . PMSA, USBTNL, PCHI , 
f PNASI 3 

QPF.PUO.UPP, 

0UPP1 , 0PPF1 . 0QPP1 , OPSP1 .OPSS1 . OPUI . 

0UFF2 . OPPP2 , 0QPP2 , OPIP2 .OPSS2 , OPBS2 , SPOTS , 

I DYNAMIC SENSITIVITY 

ULPDS . UPPDS . QPPDfi, 

S »<2UIV SCRATCH AND OBXSBDR 

BFXH, PLMDDPP , FMPf , F0L1 , POL2 , CKD, 000. EOOMX . CASIDR . XYCD*OR . 

KPHIt. LDNDOPP, UPF, MPHH. MCAA. KKAA. OA. PCDB. ABM, 

PPMP, PH FA l , PNPP. PNLMODPP. PPPI , SMPP . CTMODPF. UHP1 . 

UHFF. IMPS . XYCOR. ABKHX , PHSAP . PM FAS. MDOII . USITl . 

KDOD. KDGJ . KRDO J 

TYPS DA. ZUZA01 , ZUZR02 . ZUXA03 . IUIWJ4 . 102 AOS, 

ZUZR04 , ZUXA07. ZU2A0I . ZUIR09 . ZUZA10 , ZUZA11 S 
TYPS DB. ZUZR12 . ZUZR13 , ZUZA14 % 
f SAVED PARAMETERS 
TYPE FAAM. NODL, I. M. LANCLSS $ 

TYPE PAAM.NDDL, I. N. LUSBTS. HNNLK, FHOOS.KPNPLC I 
TYPE PAAM.NDDL. LOGICAL, N.S0PH2. SKIPS! S 
TYPE PAAM, , LOGICAL. M. RSTRT S 
TYPE PAAM. NDDC. CHARE. M. PAMAKE t 
I QUALIFY A* 

TYPE PAAM.NDDL. I. N. SEID, NTEKP LOAD. TBMPLD.DSPORH.KPC. CPC. PEID, HTtH. 

NLOOP.DSSITIA: DYAD.MFLU1D.MSTHP.TPL. OLOAD.PREQ. SDAMP. r PAM EL $ 
TYPE PAAM. HDOL, CMAAS. M. K2GG. K20G. B2QG. P2C, K2PP.M2PP. E2PP. APMCM t 
TYPE PAAM. M3 OL. LOGICAL. N. P8C0UP t 
t USE A PARAMSTtAS 
TYPE PAAM, . AS, Y, K4ROT«0 . $ 

TYPS PAAM, , CHASE, Y, AUTOS PC ■' YES '.DYNSBH-Vo* % 

TYPE PAAM.NDDL . CHASE . Y. DBXSBDR % 

TYPE PAAM.NDDL. I. Y. AS IMG , DDAMi . LMODBS . RDAMP . COUPMASS . 

MOOACC . DBDICT. DSDRPAJ . OEDRVEA. DBDftOPT. 1XTACV. 

I ACOUSTIC ANALYSIS 

PLUIDMP. STRUCTW, PAMBLMP. «rror J 
TYPS PAAM.NDDL, AS, Y.G.LPABQ, HP REO.WTMASS $ 

S LOCAL PAAAMSTBAS 

TYPE PAAM, , CHARE, N. SUBDMAP* HIDFUQ ' . APP* ' FREQAESP' . APP1« ' PRSQASSP - . 
SOLTYPa MODAL * S 

TYPE PAAM. .LOGICAL, N. ASONLY- FALSE. PS. FALSI. N0PH2 -FALSE. D0PH2M.TRUE J 
TYPS PAAM. .X.N. CKAPHO*211 , PPILE-O. CARDNO-O , LASTSBID, DROPT»0. 

MOUSET, MO A SET. NOBS IT. NOCSBT. NOGSET, NOLSET. MOOSET. NOQSET. 
HORS ST, MOSS IT. MOTBET, NOVSET.MOA. MOSET. NOPC. MOMS IT. 

NOS ASET . MOFAKST-O . NSEIC . NPEIG, NPSBIG . LMOOEX . 

I FREQ, XPRSQX . HOOC. USBTABS , MB(P»* 1 , MBWK««l .MOKOGML. 
EOLCUR, SOLPRB. (MVP ABE f 

TYPE PRAM. MkOL . AS . Y . BPPAT. BPZSRO. MAXRATIO. PAPA. PAPJ. 

BPSMALU. EPSMALC, SPSEIG. EPSAC. PRPKXVZ . RMAX . SIGMA . 
TINY. TABS S 

TYPS PAAM. NDOL. I, Y. SBQOUT, NBWBBQ. SUPER, FACTOR. >»CX. START, 

IRES, GRSPWT. SPCCIM. GPBCT. EST. USETPRT, USVTEEL. 

BAILOUT, ZNALM. PATPCM. WECI. NOSMADE. 

GBOMU.LOADU. POSTU. DBCDIAG. FIXED*, LGOIKP 5 
TYPE PAAM.NDDL, CE.Y.CM1. CM2. CK1, CK3.CE) .CRl.CRZ.CPl.CP? S 
TYPE PAAM.NDDL. QUAE. Y.PEGPST, ORCCOt*/ . DSCWVRT $ 

TYPE PAAM, .I,Y.E>rr»-l.LOOrnM).SLOOPID— l.AKSTRAIN.MAXITBR** f 

TYPE PAAM, .Z.N.MBKIP-0 S 

TYPE PAAM, , AS, N.OPASQ. XNOAMaX .0 f 

TYPE PAAM. . AS. Y. BETA-0. S,ABSN0AN«0. 05, MAXNOAM-1. 0-3 S 

TYPE PAAM. .AS. M.UGMtMAX , MUCBIIHAI, tl , *4 f 

TYPE PAAM. .CS.N.GAMOh* (1.0. 0.0). 0KBGA2.X0WGA ( 

TYPE PAAM, ,CS. N.CGMSTla ( 1 - 0. 0 . 0) ,C0NBT2« (1 .0 . 0.0) ( 

I 

S PERFORM PHASE 0 HMD 1 OPERATIONS 

< 

NLOOPaO S SET QUALIFIER 

CALL SUPSA1 /CASSCC, USBT. GOAT. KAA, DYNAMICS. OPLS. SILS.SLT. DM. 

CASES . SQSXCNS . MkTTOCC , CM . DIT. BGPDT8 . C8TNB . 1ST . MPTB , 

XYCDBS. LLL.KLL.ULL. , PCD*. SNAP. XYCOB, MkPS. SIT, FVTS. 

KPS, ESS . KSF. IDT, EDOKX.GBOM2E . GBOK3S. FOSTCOE ECTS . EFTS . 
INDIA . KELM. EDICT. GPBCT. VBLSH. FORCE. EQflXZNX, BCTX. BOPOTR. 
EILX. ..PJ.FTBLEM. iO.PA.YB.QA.POS.PSS.PL. . 

. KW, KJJ. MO, KOO. LOO, UJ. EGO, K4GC. MJJ. MPF.MLA. 

NLL , MAA , CMLAMt , C)#HA . 04IMO, JCJU1 , MA , OOAQ . KLAA , HLAA , MAA . 
MAA.KEA.BAA.MAA.MOO, .. , , GSCM4S, VCPS. 

PCDAS . KSTL . ESTM. , AOG . ONI . DMXNDX , DTI . CTINDX , DRQATN, DEQIND. 
KIM IT. BQXONO, PLIST2 . PUSTl.COMEEL. DPLDXt . OX DM IT. 

DT0S2 . DT0S3 . CTOS 4 , DSC REM. TABOSQ. DVP, CVTN. DTB, 

DBQCON. DAQPOS , DESGRD. DESOID.COOADM, DCLDXT.COM/ 
APP/DMRMO/C , RSOM. Y / SOL TYP/ FALSE /FALSI /FALSI /TRUE/ 

APACE /* .SEID /S, FEZ D/S ,MTE>#/S, LORD/S.TEWLD/S . DEFORM/ 


S , MFC/S . SPC/S, K20G/E , M20G/S , B2GG/S , P2G tZ , OYRD/S . WTH/ 

8 .KPLUI D/S. LASTS BID /S , PFILE/S , LUSBTE/S . FSCOOP/E* ROT/ 
AOTOSFC/S , N0FM2/S, DOPH2K/S, DESITER $ 

DO WHILE ( KHToNAXITIR ) i 
IF < tXTICV.O ) THEN $ 

s 

f PERFORM PHASE II OPERATIONS 

I 

PUTSY«(0. 109 | $ FORCE EXECUTION ON RESTART 

PUTSYS(0.S2) I DISALLOWS DNAP TO FIELD MOGO FLAGS 

FILE DLTaSAVB. OVRWRT/ VAFSaSA VI. OVAWRT /LAMAS ■SAVB.OVRWRT/ 

MHHaSAVE. OVRMRT/ WH -SAVE. GVWRT/XHH«SAVB. OVWRT/ 

PPF-SAVE, 0VRWAT/PSF«SAV1. OVHNRT/ PDF -SAVE, OVAWAT/ 

FOL -SAVE. OVRWRT/ Pit F -SAVE . OVRWRT/UHF-SAVE, OWRWRT/ 
CASBDR*SAVE/XYCDBDR«SAVI. CK/AWRT/UPFaSAVI . OVRMRT / 

NKAA «£AVB , OVRWRT/CSTW -SAVE . CR/RWAT/ KDQG > SAVE . OVRWRT / 
uset-savb , ovwnrr $ 

FILS K1AA-SAVE.0VWRT/KLAA*SAVE.GVRWAT S 

FILE KEUBTL-SAVE. OVRWRT/KDICTKL-CAVI . CVRNRT S 

FILE MUCNr-SAVE.OVRWRT/KESTNL-SAVl, WRWRT/ES'INL-SAVE . OVRMRT i 
FILE EST .SAVE. OVRMRT / BSTL-SAVS , OVRWRT/UCNI- SAVE. OVRWRT S 

FILE CAS BSX2Ha SAVE $ 

FILE XUZR13-(A/RWAT/IUEA14«OVRNRT/IUZA1$«OVRWRT S 

FILE KAA«SAVB.OVRMAT/MAA«SAVE.OVPWRT/CPECT»SAVB,OVRWRT 5 

s 

% NONLINEAR OUTPUTS 

< 

IF ( NOT GOPM2 ) CALL SRAPK2 //'PHASE I /O i ERROR IN PHASE I 
PAAAML CASSCC/ /’DTI ’ /1/20C //S, N, MANUAL t FLAG FOR MANUAL RESTARTS 
IF < MANUEL* ►- 1 AND LASTEBIDoO ) EXIT S EVERYTHING IS OK 
$ IN PHASE I BUT USER NADI NO PHASE I REQUESTS FOR R.S. 

j AND REQUESTED MANUAL ANALYSIS ONLY 

t 

I PERFORM EIGEN SOLUTION ON RESIDUAL STRUCTURE 
S 

CALL PMLUSIT USBT/ /S . NOASET/S . NOBSET/S . MOCSET/S .MOOSET /S . NOLSET/ 

5 . NOOSET/S , NOQSET/S . MORSET/S . MOSSIT/S . M0TS1T/ 

6. NCWSET/S.NOA/B. NOSET/S, MORC/S, NOHSET/S, NOUKBT S 
IP | MOOSET* * 1 ) THEN $ 

MESSAGE //' DMAP FATAL MESSAGE 901S (MLSEMFPBQ) 

' THE USBT TABLE IS NOT PRESENT . ' $ 

END IP S 

$ 

DBXSSDRaDBXSBDR S SET LOCATION FOR DATABLOCKS NEEDED FOR EXTERNAL SB DR 

t SIT UP TABLES FOR DYNAMIC ANALYSIS - UN 3.3.24 

DPD DYNAMICS . GPLS.SXIS, USBT. SLT, PG/GPLD. SILD. USBTD. 

TFPOOL . DLT . PSDL. PAL . NLPT . TAL . BED . EQ01N / - 1 /S . N . L USBTD/ 
O/O/O/O/O/ 0 /l /123/S. H. NOUS S 
S 

PARAML CASES// DTI 1 /l/ IS// S.N.TFL i 

PAAAML CASES// DTI' /1/13H////S.N.K2PPA $ 

PAAAML CASES// ‘DTI '/l/I 40 ////£. H. K2PPB $ 

K2PP*K2FPA 4 K2PPB 

PARAML CASES// *DTI' /1/141/// /6.N.X2PPA $ 

PARAML CASES// DTI' /1/142////S.M.M2PPB S 

M2PP>M2PPA 4 M2PPB 

PARAML CASES//DTI/1/143////S.N.B2PPA $ 

PARAML CASES// *DTX'/1/144////S.N.B2PP6 $ 

B2PP*B2PPA 4 B2PPB 

PARAML CASES// ‘ DTI ’ /1/14R//S . M.SDAW* 5 

PARAML CASSS//DTI/1/13 //S.M.DLOA D t 

PARAML CASES// 'DTI* /1/14 //S.N.FRSQ $ 

PARAML CASES//’ DTI' /1/22S//S.N.MITHF S METHOD ( FLUID > 

* 

PARAML VCPS/ /' FEES' ////E, M.MOVCFS $ 

IF ( MOVCP*>-l ) THEM | 

BQUIVX VGPE/VAFS / MOA i 

IF ( NOA» - 1 ) UPARTH USET. VGPE/VAPS. , . / ’G' /’ A’ /'S’ /I | 

PARAML VAPS/ / 'TRAILER ' /S/S ,N, NOP ASET S 

PARAML VAPS// TRAILER' /4/S. N.PRBC t 

NOP ASET •MOPAStT/PRSC i NO. OF FLUID DOFS IN A- BET 

BM>XP S 

NOSAfET«NOA£ET>NOFASr.' $ NO. OF STRUCTURE DOFS IN A-SET 

rsm (NOSA£«T>0 AM> N0FASIT>0) t BOTH FLUID AIO STRUCTURE EXIST 
t 

IP ( MOT MCPH2 ) THEN S 

* 

t COMBINE STATIC AND DYNAMIC PORTION OF THE RESIDUAL STRUCTURE 

S 

CALL SBOOA GOAT/GOA /MOOSET /NOQSET /HOTS ET ( TRAMS FORMATZOM 
IF ( NOQSST—1 ) THEN t 

BQUIVX KAA /MKAA/ ALWAYS f 

ELSE IP < N0TSET--1 ) THEN $ 

BQUIVX KLAA/MCAA/ALt*YS $ 

ELSE ( 

ADDS KAA. KLAA. . . /MCAA I 
BM3IF $ 

ADOS MLAA.MUl, . ./W4AA t 
t 

BQUIVX GOA/GOO/NOUE ( 

BQUIVX CM/GMD/N0U1 $ 

BQUIVX B4AA/K4DO/NOUB t 
IF | MOUB>- 1 ) THEN t 

UMRCE1 UsrPD.GM, . ,/CKD/'NE /’M'/’B’/2 $ 

UNBRCB1 USSTD.GOA. . , /OOO/'D' / 'A' / B ’ /2 « 

UWBRCB1 UESTD.K4AA, . ./MDO/'D’/'A'/'I* I 
INDIP $ 

DEVIBW XISTM, « ESTM. (MIERE KLOOF . SLOOP ID) $ 
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DSVIBW YUGNI a MUGNI IMIIM NLGOP . SLOOP ID) $ 

DEVIBW YPGW a POl (Mil RE NLOOP a SLOOPID) $ 

DBVI1W Y1STNL - 1STNL (WERE KLOOF a LOOPID) % 

4 

$ INITIALIZE SOL PRS AND SOCCUR 
SOL PUS* LOOPID i 
SOLCURaSOLPRS 4 

USTUTa ( LOOPID>tJ J | LOOPIDaO FOP COLD START AND >0 FOP RESTART 
PAPAML YUOMI// 'PUSS' / / / /S . N. NOYUQU 5 

IP { NOYUGNX>- X ) THIN S 
NLOOP-O | 

SQUIVX YUCNI/VUC/NOA 5 

squivx ypgni/plx/noa s 

IP ( NOA»- 1 ) T»ON 4 

UPARTN osst, Y0C»I/LHJC, , , / 'O' / 'A' /-S-/1 s 
6S02 USST. «. , KFS.COA.DM. YPCNI/, . .PAX.PLX $ 

SNDIF S NQAxl 

SQUIVX ULX /(JUITI / ALWAYS $ 

SQUIVX PLX/IFO/AU*YS S 

COPY IPD/IPDH | 

COPY ULNTI /ULNTH $ 

BM>XP $ NOYUGNI>- 1 

DSL STS /HESTNL J 

COPY YBS1ML / MSTWL 1 

DSVISW CNO-CN OMBRE NLOOP-O) f 

DKVISW RDOOaMKAA (WHERE KLOOPaO) | 

4 

t IP THERE IS NO SXTRA POINT. A- KIT a D-SST 

I 

DBviBw usrrOaussT {wnx nloop-0) $ 

t 

IP ( KNT>- 1 ) THEN $ 

$ 

4 TOP OP NON-LINEAR LOOP 

< 

SOLCURaSOLPRX+1 $ 

NLOOPaSOLCtm-1 $ 

$ 

$ FIND FOG POP NONLINEAR BLSMSNTSj ALSO KDGC IP GBOM NONLINEAR 

5 REDUCE THESE TO A-SIZE 

4 

BQUIVX KDDO/KAA/ ALWAYS $ 

IP I LGDISP*- 1 | THEN $ FOP CBOKBTPIC NONLINEAR 
* GENERATE KDJJ AND REDUCE TO KDOD 

5 

IP ( UHVFABS*-! ) THEN S JUMP ONLY FOP COLD START W/O IC 
$ EXTRACT AND FORK SJLNI FROM ULNTI 

PAPAML ULNTH// TRAILER' /1/S.N.DCOL f DCOL S IN ULNTI 
IP ) DCOL*2 ) THEN ( 

SQUIVX ULNTH /ULWTX /ALWAYS $ 

ELSE 4 

IP | NBWP>-1 ) DCOL* DCOL - 1 4 SUBTRACT TO BACK UP FOP ACCEL 
DCOLaDCOL-1 S NOW DCOL CORRESPONDS TO DISPL 
S EXTRACT DCOL-TH COUREi VICTOR 

HATMOD ULNTH /ULNTX , / I/DCOL J 

*»IP 4 DCOL< 2 

BQUIVX ULNTX/ULNX /ALWAYS $ 

PORI USET, .ULNX, . .GQA.GHO, , . . /ULNIX. . /1/ FR10' $ 

I 

$ END OP FORKING DISPL VECTOR POP X MATRIX 

( 

CASE CASBCC, /CASSXX/' FREQ' /S,N. REPEAT /0 $ 

V EC PLOT NUCM1 , ZUZR11 , BQBXZNS, CSTMS. CAS EX X. /UCVBAS /l) /Q/ J S 
BQUIVX EGPCrrs/ZUI RIO /ALWAYS $ 

DELETE / EGPDTS . . . , | 

METHOD UCVBAS, ZUZRIO /BGPDTS. /II $ 

DELETE /ESTHL, BST, . , $ 

TA1 KPTS.ECTS, EFTS, BCFDTS. SILS. BTT. CSTMS. DIT/ 

SST , ISTNL, GEI , GPBCT. SSTL/ 

LUSm/ 1 /S , N . NOS IMP/2 /E , N. NOGEM./SBID/ 

I /NLA YE RS $ 

EHS SSTJ*,. CSTMS, MPTE, DIT, -ULHIX. . /ICDEUPi, KDOICTH, . , , / 

1/0/0// / / / ///////// /LAMGLBS $ 

^ ENA GPECT, KDOICTH. KDBLMH, BGfVTB. SILS, CSTNE/KDJJ. /-I $ 

$ REDUCE DIFFERENTIAL STIFFNESS AND COMPUTE 
* TANGENTIAL STIPPWPS POP GEOMETRIC NONLXMAP 
I 

EQUXVX KDJJ / KDMt/NOMSET 4 

IP | NGMSlTk-1 ) MCI 2 USBTO ,040, KDJJ . , , /KIMN, , , ( 

EQUXVX KDNN/XDPP/NOSSET $ 

IP < HQSSST>- 1 ) SCSI USBTO, KDNN. .. /KDPP J 

EQUXVX RDPF/KDU./MOOPET $ 

IP ( MOOSBT>-1 ) THEN f 
f ROOM, MST EE NULL 

UPARTTN USBTO , KDFF /KDCL . , . KDOO/ • F ' / ' A ' / ' O' | 

PAPAML XDOO/ /' TRAILER/ S/S. N.ONOTMULL//S.N, HP 4 
f ERROR 4423 a> NONLINEAR ELEMENT ATTACHED TO OMITTED DOF 

IP t ON01MA.L»0 > CALL ERRPH2 / /8UBCMAP/44 23 J RPSRRM 
ENDIF f NOOS*T>-l 
BQUIVX XDLL/KDDO/ ALWAYS t 
BQUIVX KDCD/KROO/ALkMYS | 

1*0 IF t UHVFABS*-1 
ENDIF I LGDISP»-1 
■MJtr « XNT » *1 11/20/R1 

$ 

IP < NOQSETa- l | THEM | 

ADOS RAA , ROLL. , , /MKAA// ( 3 . 0, 0 . 0) $ 

BLEB IP ( NOTStTa-i ) THEN ( 

ADDS KLAA.KDLL, .. /MKAA// <3. 0.0.0) f 
BLEB S 

ADOS KAA.KLAA.KCCL. , /MCAA/// <3.0,0 .0) S 

EM>IP | 

< 

IP { DOPH2M ) <Wl«c« /pJu, Iwtr.phM $ 

d*l«C« /1«M». ... t 
CALL MODBPSRS. 


MR. US BT. DM. CAS IS. DYNAMICS, KMAA .MCAA, CPLS, SILS. BED 
BQBXINS.VAFS/ 


$ 

IP 


PHA. LAMA. PHPA, LANAF. PHSA, LANAS/ 
NORSST/NOLSBT/AS D*G /NQASET/ FS / NOFASET/NOG ASET/ 
METH/KBTHF $ 

(KETH>0 AND NO6ASBT>0) THEN $ 

IP { FE ) THEN S 

MERGE PHSA VAPS/FHSAP/l 5 


BQUIVX PHSA/PHSAF/-1 $ 

ENDIF $ 

VDR CASES, EQEXINS.USBT.PKSAP, LANAS,. /OfHSA. /'RBIG'/ 

* DIRECT ’ /0/S. N, NOSVBCT/ 1/1 $ 

IP ( NOSVBCT * -1 ) OFF OPHSA//S.N.CARENO 4 

I SV1CTOR OUTPUT FOR RESIDUAL STRUCTURE ONLY 

ENDIF S NSTH>0 AND NtWASBT»0 

$ 

IP (HETHF>0 AND N OF AS ST >0 > THEN 4 
IF ( FS | THEN $ 


ELSE f 

BQUIVX PHFA/PHFAS/-1 S 
ENDIF $ 

VDR CASES . BQBXINS. USBT. PHFAS, LANAF . , /OfHFA, / ’ REIG ' / 

DIRECT' /0/S. N. NOSVBCT/ 1/1 $ 

IF l NOSVBCT > -1 > OFP OPHFA / / S . N , CARENO S 

S S VECTOR OUTPUT FOR RESIDUAL STRUCTURE ONLY 

ENDIF S XSTHFoO AND N0FASHT»0 

S 


DELETE /ESTNL S 

DELETE /MCK, BHH. KHH. PHIDH, PHDH/ $ 

CALL OKA CASES , KATPOOL, EQDYN, TFPOOL 

MKAA, GOD. GMD.USBTD. MKAA, BAA, JUDO, PHA, LANA, DIT, VAFS, 
PHSA, PH FA. LAKES, AGG/ 

. .CMDD.CK2DD.CB2DD.MtH.BKH. KHH, PHDH, VPHFS/ 
lusbtd/soltyp/app/noub/noa/g/o . /U . /TRUE/ 

0 /LNO DBS /L FREQ /H FREQ/ S , FMODB/ KDAMP /FALSE/ FALSE /FS S 


DELETE /PPF. PSP. PDP.POL.PHF/ $ 

DELETE /(JtVF.UHF, . , S AVOID OUTPUT TWICE 
CALL MFREQRS CASES . USBTD. DLT, FRL . GKD, GOD, DIT, fHIW , KHH. BHH MtH 


PPF.PSF.PDF, POL, fHP.UHP/ 

SOL TYP/APP/NOA /FALSI/0. /0./0. /0 S 

$ 

SNDIP S NOT (NO PH 2 1 

$ 

PARAKL UHP// 'PRES* ////S.N, NOUNP $ 

IP (NOCHPa.l) EXIT $ IP NO SOLUTION TO PROCESS. 

f 

PARAKL PHPA// 'TRAILER' /1/S.N, NFBZC//S.N.NOPHF A 5 

PARAKL PHSA/ / 'TRAILER' / 1/S, N.N6BIG//E, N.NOPHSA S 

IF \ NOPHPAa-1 ) THEN J 

BQUIVX UHP /UHPS/ ALWAYS 5 

ELSE IP { MOPHSAa-1 ) THEN $ 

BQUIVX UHP /UHF F / ALWAYS J 

ELSE S 

5 EXTRACT FLUID AND STRUCTURE SOLUTIONS 

PARTN UHP, . VPHFS/UHFS, UHFF. , /I $ 

ENDIF $ 

$ 

IF ( NOPJff A>- 1 ) CALL VDR1. 

CASES, BQOYN, USETD.UHPS, FOL, XYCDfcS , PSDL , DIT, // 

^ RFP/SOLTYP/S , CARDNO/S , PFILE / ' HSXT ' /PMDDt/ FALSE { 

IF ( NOPHPA*- 1 ) CALL VIM1, 

CASES, EQDYN. USBTD. UHFF, FOL, XYCDES, , PSDL. DIT, // 
RFF/SOLTYP/S, CARDNO/S, PFILE /'HSET ' /FMDM/ FALSE $ 

IF | FS 1 DOMMta-1 S 

IF ( DDRMtxQ ) APP1 a 'MOIEXG ' $ 

I 

DELETE /POL1 . PPF1 . ULF. PSP1 , UHP1 4 AVOID OUTPUT TWICE 
CALL MOOACC CASES . POL. UHP, PPP. PDF, PSP, PHDH. 

USETD.CK2DD.CB2DD.OCiD, LLL, CM/ 

FOL1, PPP1, ULP, PSP1. UHP1/ 

^ APP/APP1/NOUI/J4DDACC/SOLTYP/FALS1/ NO ' J 

t It MMBBR OF FREQUENCIES FOR FLUID. PUIIM*, IS -1 THEN 
t COMPUTE FLUID MODAL PARTICIPATION FOR ALL FREQUENCIES 
PARAKL F0L1//' TRAILER' /1/S.N. NOFRSQ $ 

IF ( FLUIDWUO OR FLUIDMP>NOPRjQ > PLUX0MP-KOPRBQ $ 

IF ( FLUIDHPX3 ) MESSAGE // 

' DMAP INPORMATZON MESSAGE 4013 (NLSBMFREQ) - FLUID MODAL'/ 

' PARTICIPATION IS REQUESTED POR' /PUTIDNP/ ' FREQUENCIES . ' f 
J IF NUtEBR OF FREQUENCIES POR STRUCTURE. STRUCTMP. IS -1 THEN 
4 COMPUTE STRUCTURAL MODAL PARTICIPATION POR ALL FREQUENCIES 
IF ( STRUCTMP* 0 OR STRlXTHP>NOFRSQ ) STRUCTMP* NOP REQ 4 
IF ( STRUCTMP>0 ) MESSAGE // 

DMAP INFORMATION MESSAGE 905 3 INLSBKPRBQ) - STRUCTURAL AND LOAD’/ 
' MODAL PARTICIPATION IS REQUESTED FOR • /STRUCTMP/ ' FREQUENCIES. ' | 

IF i STRlCTKP*0 OR FLUIDMP>0 OR PANELMPvQ ) THEN 4 
4 COMPUTE MODAL PARTICIPATION IP REQUESTED 

PARAKL CASES//' DTI '/-!/ 1C l//S,N.YISOG 4 

NOOGaLTOI (NDTL (YB80G) ) 4 

IP ( NOOO-1 | HATMOD BQBXINS, USBT.SILS. CASES , / 

OG. /17//1/S.N.NOOC S 

EQUIVX VAPS/QA/NOOO 4 
SQUIVX PHPA/FHPA1/MOOO 4 
IP ( NOOC>-I | THIN f 

UPAR3W USBT. OG/OA. . . / *G' / ' A' / 'S’ f\ $ 

PARTN OA. .VAPS/.OAP. . /I 4 
PARTN PHPA. .0AF/.PHPA1. . /I 4 
DO IF 4 


A2 


4 

DELETE /UHPE.UHFF. , . $ 

IP < MOPHPA--1 > THIN $ 

BQUIVX UHF1/UHPS/ALVAYS $ 

BUB IT ( NOPKSA--1 > THIN $ 

BQUIVX UHF1/UHFP/ ALWAYS $ 

BUB t 

f EXTRACT FLUID AND STRUCTURE SOLUTIONS 

PANTO UKPI, , VPHPS /UM PS , UMFP , , /I $ 

ENDIP % 

4 

IP ( £TRUCTKP>0 J HOOACC CASES . POL , UHF, . PKF. / POLX. UHFX PHfl , /APP $ 
IP ( N0PKSA.-1 AND STRUCTMP>0 ) TH8N $ 

4 RESIDUAL MODES API ALL FLUID MODES 

BQUIVX MHH/KFHH/ ALWAYS 4 

BQUIVX KHH/KPHH/ ALWAYS J 

BQUIVX BHH/BFHH/ ALWAYS 4 

BQUIVX PHP1/PPMP/ ALWAYS $ 

BUB IP ( NOPHPA*-l ) THEN S 

( RESIDUAL MODES ARB BOTH PLUID AND STRUCTURE MODES 

IF { STRUCTMP>0 ) THEN $ DO THIS ONCE ONLY IP REQUESTED 

PARTO MHH , VPHFS , / , AM . . KFHH S EXTRACT COUPLING AND PLUID MASS 

PARTO XHH , VPHPS . I .. . KFHH $ EXTRACT PLUID STIFFNESS 

PARTO BHH. VPHPS. /,, - BPHH $ EXTRACT PLUID DAMPING 

PARTO PHP1. .VPHPS/, PPHP, . /I % EXTRACT PLUID LOADS 
END IP 4 STRUCTHP>0 

IP < PANSLMP>- 1 AND KFNPLC>0 I THEN $ 

GP4 BCTS, BGPDTS, BQBXINS, BDT, SXLS/ 

HFNSLT. KBQAST. MNRTAB/ 

MFNPLG/fi. N.KUMPAN/S, N, MATCH $ 

I PANEL- I $ 

DO WHILE l IPAN1L< *NUKPAN ) % 

A CMC MPNSLT. BGPOTS. CSTMS . SILS. ECTS . KBQAST, MNRTAB/ 

ABB/LUSBTS/MFNFLG/NUMPAN/S, N. PANAMB/ I PANEL /MATCH $ 

$ THE ABB MATRIX IS NEGATED BBLCW, BY 2ND PARAMETER OF SMPYAD 

KATRBDU ABE.USBT, GM. OOA/ABBA/S . N.NOABE S 
PARTO ABBA, VAFS. /, .ABBSF. S 

SMPYAD PHFA1.A8ESP, FHSA, , , /ABBKX/3/-1///1/1 $ 

IP ( LM3DES.0 AND LMODB£< NC BIG > THEN $ 

LM0D8X -NS EIC- INODES S 

KATGEN , /LMODBS/6 /NSBIG/LMODES/LMODEX $ 

PARTO ABSHX . LMO DBS , / AB EH , . . / 1 $ 

ELSE $ 

BQUIVX ABSHX/ABKH/-1 $ 

END IF $ 

I PANEL -I PANEL* 1 $ 

END DO $ I PANE L< -NUMPAN 
END IF $ PANELMP*-1 AND MTOPLG»0 
ENDIP $ 

4 

IPRBQ-1 4 

DO WHILE ( IPRBQoFLUXDMP ) $ 

MATMDD UHPP /UHPPI, /1/IPREQ $ EXTRACT A COLUMN ( FREQUENCY l 

PARAML UKPPI// 'NULL ' / ///S . N, ZUHPI $ 

IPREQX- IFREQ* 2 $ 

PARAML POLl/ /' CTI ’ /0 / IFREQX/S, N. OPREQ $ 

IP ( 1UHPI ►- 1 AND 0PRBQ>0. ) THEN $ 

MATMDD UHPPI /UKPFD./2* $ DIAGONALIZE COLUMN 

KPYAD PKPA1.UHPFD. /FMPF $ FLUID NODAL PARTICIPATION 
MESSAGE //* DMAP INFORMATION MESSAGE 9054 (NLSBMFREQ) - / 

' PLUID MODAL PARTICIPATION FACTORS FOR FREQUENCY. /OPREQ/ 

' SHOW BELOW IN MATRIX PLMOOPP J 
4 IP NSBDBD. EXPAND TO A-SIZt SO THAT HA TCP R UNDERSTANDS 

BQUIVX FMPF/ FLMODPP/ NOG ASET $ 

IP { NOS AS IT*- 1 ) MERGE ..PMPP OA/FLMOOPP/X % 

KATGPR GPU, UEIT. SXLS . PLMODPP/ / 'H ' / * A ' $ 

ENDIP 4 ZUHPI >- 1 OR OPREQ>0. 

IFREQ. IFREQ* 1 $ 

DODO $ IPRBQoFLUIDMP 
4 

IPRBQ-1 ( 

DO MilLI ( I PRBQ< ■ STRUCTOP ) $ 

MATMDD IMPS /UHPSI, /!/ XPRBQ $ EXTRACT A COLIMN OP DISPLACEMENT 

MATMDD PPHP /PFHFI. /!/ IFREQ $ EXTRACT A COLUMN OP LOAD 

PARAML UKFSX// 'NULL* ////8,M. ZUHPI 4 

IPREQX -IFREQ* 2 t SET UP POINTER TO FREQUENCY LIST 

PARAML FOCI/ / ' DTI ' /0/IFREQX/S, N, OPREQ $ EXTRACT INPUT FREQUENCY 

IP ( ZUHPI*- 1 AND OPRBQ»0. > THEN % 

MATMDD UHPSI /UKPSD. /2fl $ DIAGONAL! ZI COLUMN 

MPYAD AN.UMPSD. /A1 $ MODAL COUPLING X DIAGONAL COLUMN 
4 CONVERT FREQ TO COMPLEX -OMEGA* *2 

0MBGA2-CMPLX {- ( (2 . *PX (1 ) * OPREQ > ■ *2| J 4 
ICMEGA.O<PLX(0 .. (2. •PIUI • OPREQ ) > 4 
4 -OMEGA* *2 (MASS) « I -OMEGA (DAMP INC) * STIFFNESS 

ADOS MPMM, BFHH, JCFHH .. /22X/GMBGA2 /IOMEGA $ 

DECOMP Z2X/Z2L , Z2U//////S.H, SING $ DECOMPOSE Z2X 
IP ( SING*- 1 ) THEN 4 IP Z2X IS NOT SINGULAR THEN PROCEED 
t 

I COMPUTE STRUCTURAL PARTICIPATION 

t 

MESSAGE //' DMAP INFORMATION MESSAGE »0S4 (MLSEMPR1Q) - / 

' STRUCTURAL MODAL PARTICIPATION FACTORS FOR FREQUENCY. ' / 
OPREQ/' SHOW BELOW IN MATRIX STOOOPf • J 
PBS Z2L.Z2U.A1/XSKPP 4 SOLVE FOR PARTIAL PARTICIPATION 

OMEGAS - -OMBGA2 { CONVERT OMEGA “2 ABOVE TO POSITIVE 
ADDS PKFA1 .... /W3 PH FA /OMEGA 2 $ OMEGA* *2 X FLUID MODES 

MPYAD W2PHPA. XSKPP, /SMPP 4 STRUCTURAL MODAL PARTICIPATION 
$ IP NEEDED , EXPAND TO A -Sill SO THAT MATCPR UNDERSTANDS 

EQUXVX 8NPP/STMDDPP /NOP ASET $ 

IP ( NOPAS ET>- 1 ) MERGE ..SMTP QA/STOOOPP/ 1 J 

4 PRINT STMOOPF IN EXTERNAL ORDER 

MATCPR GPLS. US *T. SILS . STMODPF/ / ' H‘ / ' A ' $ 

< 

4 COMPUTE PANEL PARTICIPATION 

4 

IP ( KTOFLC-0 AND PANBLMF>-1 > THEN f 
I PANEL- 1 | 


DO WHILE ( I PANBL< -NUMPAN ) $ 

MPYAD ABEH.UHPSD, /A1P J 

PBS Z2L. Z2U. A1P/XPMPP S 

MPYAD W2PKPA , XPMPP . / FKPP { PANEL MODAL PARTICIPATION 
t IT NEEDED. EXPAND TO A -SIZE SO THAT MATCPR UW>*RSTANDS 

BQUIVX PMPP/ TO LMO DPP /NOPAS BT S 

IP < NOP ASST*- 1 > MERGE .. PMPP, ,,, OA /PNLMODPP /I $ 

4 PRINT PNLMODPP IN EXTERNAL ORDER 

MESSAGE //' DMAP INFORMATION MESSAGE 9054 (NLSEMPRIQI -'/ 
' PANEL MODAL PARTICIPATION PACTORS FOR PANEL • '/ 
PANAMB/' AT FREQUENCY -'/OPREQ/' SHOW BELOW IN'/ 

' MATRIX * f PNLMODPP' $ 

KATGPR GPLS. USBT, SILS, PNLMODPP // 'H'/ 'A' % 

I PANEL. I PANEL* 1 $ 

BNDDO 4 I PANEL-. NUKP AN 
ENDIP $ MTOPLG>0 AND PANELMP>-1 
4 

$ COMPUTE LOAD PARTICIPATION 

4 

MESSAGE //' DMAP INFORMATION MESSAGE 9054 (NLSBMPREQ) -*/ 

' LOAD MODAL PARTICIPATION FACTORS FOR FREQUENCY- ' /OPREQ/ 

’ SHOW BELOW IN MATRIX LDMOOPP' $ 

PBS Z2L, Z2U , PFHPI/XLKPF $ GET THE INTERMEDIATE VICTOR 

KPYAD PHPAl.XLMPP, /LMPP % 

$ IP NEEDED, EXPAND TO A-SIZB SO THAT KATGPR UNDERSTANDS 

BQUIVX LKPP/LDMODPP /NOPAS BT $ 

IP ( NOP ASET*— 1 ) MERGE , , LMPP OA/LDMOOPF/1 S 

J PRINT LMO DPP IN EXTERNAL ORDER 

MATCPR GPLS. USBT. SILS, LDHODPP// 'H' / A' $ 

ENDIP $ SING* - 1 
ENDIP $ ZUHPI >-l OR OPREQ >0 . 

IFREQ. IFREQ* 1 $ 

END DO $ IFRBQoSTRUCTOP 
ENDIP $ STRUCTMP>0 OR FLUIDMP>0 
SNDIP $ SXTKV.O 
$ 

4 PERFORM PHASE III OPERATIONS 

4 

BQUIVX POL 1 / POL2 / DDRMM S 

IP ( DDRW>-0 ) BQUIVX LAMA/ POL2 /ALWAYS S 
$ 

DELETE /CASEDR , UPP, QPF. XYCDB DR, PUG 5 

CALL SUPER! CASBCC. . ULP , POL2 . PPP1 . PSP l , . FOLI . UHP 1 . DLT . 

4 CRX . CASBBK1 , QRG 

PCD8. XYCDB. POGTCDB. /KUC.NI . 

UPF , QPP, CASE DR , XYCDEDR, PUG, . . 

OUPF1 . OPPPl, 0QPP1. OFBF1 . OFBSt . OPBE1, , 

0UPP2 . OPPF2 . OQPP2. 0FEF2 , OFES2 , OFEB2, , 


4 /CYCLIC /NLHBAT/AERO 

APP / A PP 1 / RSONL Y / FALSE /FALSE /PALSI/O/PS/ 

S , PPILB/S , GARDNO/O / PALS B/RMSTRA IN /BETA /LMOOES $ 

IP ( MAX ITER > -It THEN 4 

ADD MUGNI , / ZUZR1 3/11.0.0.0) $ 

BUB IP ( KVT • -1 J THEN 4 
C0NST1 «OU> LX I BETA .0.0) 5 
ADO MUGNI , /2U2R1 3 /COMST1 $ 

IUE S 

A 00 ZU2R14 , MUGNI /DIP/ ( - 1 . Q, 0 . 0t i 
NORN DIP/DI PN///S.N, XNORM 5 

MESSAGE //' ITERATION NO. ( - 1 . 1 ifi—r I - ' /KNT/ ' XNORM- '/XNORM 4 
non ZUZR14/ugnin// fm.n , ugniMX $ 
non aucni /Kuonin ///•,!). Kugnimitx $ 

IMMN / / ’ U0Tli.JM»X» ' /U9D1.M % 

■*«»Ag« //'Bugnisax- ’ /augnitux 4 
24. (UCTItMAX-MUCWDIAXt/KUGNIMAX $ 

IP ( ABS< 24 ) < ABSNORM ) ICNT-MAXITBR* 1 $ 

MESSAGE //‘ABSOLUTE NORM • /Z4 $ 

C0NST1.CKPLX (BETA. 0.0) $ 

C0NST2-CMPLX( 1.0 -BETA, 0.0) $ 

ADO MUGNI , ZUZR14 /2U2R1 J /CONST 1/ CONST 2 5 
ENDIP $ 

DELETE /MUGNI $ 

BQUIVX 2UZRU/ MUGNI /ALWAYS 4 
BQUIVX MUCNI/2UZR 14 /ALWAYS $ 

BQUIVX BCPOTS/ZUZRll/RMT 4 
ENT-RNT* 1 4 

dfcviww tHtnliMthl (wKar* I9di1p.ll 4 
EM? DO $ 

NLOOP.O 4 SET QUALIFIER FOR PST AND DATA RECOVERY 

OUTPUT2 //- 1/1 2/LA BEL- PDA LA BEL 4 

8 DR 2 CASEDR, CSTMS. . .IQEXZNS. , , POL2 , BGPDTS , , .WON I, ./ 

. , OUGV1PAT. . ,/APPl $ 

OUTPUT2 OUGV1PAT//0/12 $ 

ovmm //-9/12 4 

IP ( DYNSIN. YES' ) THEN 4 

EQUXVX UPP/UPPDC/DORNM 4 

BQUIVX QPP/QPPDS/DORMM $ 

FURGEX /ULPDS /DDRW C 

IP < DDRW>-1 ) THIN $ 

DELETE /POLl , PPP1 . PSP1, UHF1, 4 AVOID OUTPUT TWICE 

CALL MOOACC CASES. POL. LMP.PPP.PDF. PSP, PHDH , 

US1TD / 

POLl. PPP1, ULPDS. PSP1, UHP1/ 
APP/APP/MOU1/-1/SOLTYP/PAUE/ YES ' | 

CALL SUPERJ CASBCC. .ULPDS. POLl. PPP1.PST1, . 

POLl, UKPI. DLT PCDB, XYCDB. POSTCDB, /MUGNI. 

UPPDS, QPPDS. CASEDR. XYCDBDR. PUG. . . 


APP /APP /RS ONLY /FALSE /PAUI /FALSE/ 2 /PS/ 
B,PPIU/S.CARDN0/0/PAUB/RMSTRAIN/BrEA/LM0Dt8 4 

ENDIP 4 DDRW*- 1 

PUTSYS ( 1 , 109) S FORCE EXECUTION OH RESTART 

CALL DTOSEN UPPDS. QPPDS. ECTS, EFTS, EST.ETT. DIT. MPTB.GPBCT, CASES, 

CSTMS, BQKXINS. POLl, . BCPDTS, PPP1. VBLSM. IDOWt . CBOtOS . S ILS . 
EDICT . KKLM . USBTD , GMD . GOD . TOOK . KMH , EMM , MM . . 
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usbt.qoat, ot.xrs.KSfi, , .wqato.dbsind.cbomss.slt, 

XINIT, ICIDWO. PLISTJ , PLIST3 , C0HS6L. DPLDXI , DXDXIT, 
DTOS3 , DT083, DTOC4 , DPCR1H, TABCBQ, DVP. DVTM. OTB, 
D8QC0N, DBQPOS. DBSCAIX DBSGID. COOAON, DCLOXT, CCW// 
SBID/G/NOUl/WP/SOLTYP/0 . /0 . /0/WTMAJ»S/ 

covpk«s/k6act s 

BNDIP S DYNSIM* ’ YB* ' 

t 

IF ( DBCICT>-2 ) DBOIK t /DBORPftJ t DeDPVBB/DBOBOPT J 


Appendix B Additional User Defined Parameters 


PARAMETER 

NAME 

DEFAULT 

DESCRIPTION 

ABSNORM 

0.05 

Absolute norm 

BETA 

0.5 

Specifies control factor for converge enhancement. 
Ranges from 0.0 to 1.0 

KNT 

-1 

Sets iteration counter 

LGDISP 

-1 

Selects large displacement effects 

LMODES 

1 

Requests the number of modes (uses with EIGRL 
card selection) 

MAXITER 

5 

Requests the maximum number of iterations 

MAXNORM 

1.0E-3 

Defines converged rms displacement norm 

RM STRAIN 

-1 

Requests the rms strains 

XNORM 

I.0E-3 

Sets norm of rms displacements 


B1 




























Appendix C SUBDMAP SEDRCVR 


imauf sedkvr ucvs. qos, bcfets, iqixink.cstmb. cxp«d«.kpts. dit, 
itt.olbi . pji , est, xycdbdr. cioms. cecmxs , postcdb, 
ects , cpu , efts , s iut , north . xilm . edict . cpect. velem . 

FORCE , XYCDB, PC . PCDBDA . USBT. SLT . UNVF , OLBH . PCO. OCT , 

PPL . SP8EL, DYNAMICS , CRX. QOB, OBFMXX. QSSMLXX/ 
MUGHl.PUCV. 

QUGV1 , QPC1 , QQC1 , QSP1X . QBE IX . QETR1 . QOS 1 , 

QUGV2 , QPG2 , QQG2 - QEP2 . QBS2. QCTR2,QGS2. 

QBS1M, QBS1C, QCTR1X, QSTR1Q. ICPBT* . BGPSF , GMRGY1 . QCPf B1 / 
GRDPWT/APP/APP1 /NOCCMPS /CURVPLOT/PFI LI /MUWOUT1 / 
MUM0UT2 /> XCER1 / BICER2 /NUWOUT/ » IGEA/SRTOPT /LSTRN / 
simLTYP/cuiw/<xrTOPT/og mzNTrrs/siM/sic /s iam/ 
SUO/DOPT/TI NY /NOSLOP /NOSLOP/TABS /SBID/S1 /CARDNO/ 
PDRHBC/SCRSPBC/RSPBCTRA/R*PRINT/TABID/ INREL/CPPDR/ 
NLHIAT/ AERO/ ZCYCLIC / PS/OBCMU/ LOADU / POSTU / 
DBCDIAO/DBCPAOO/POST/CP/DBCOVWRT/ 

OUMIT2/OTAPB2 /OQC/OUQ/OBF/OBS /OEB/OCMP/OOPS / 
OESR/OUNU/OGFP /OUGCORD/ DSEITER/RMETRAIN /BETA/ 

LHDOES 5 

S 

TYPE DB BGPDP, HOSP 1, OBPI.OBFIX, OB P3 , OBP2Y, 0SS1 , 

OBS1G.OES1N. OBSIX, OES1X1, OBS3 . OBS2Y. 0CPPB1 , 

OGB 1 , 0GS2 - 0NRGY1 , 0PC1 , 0PC2 , 0QC1 , 0QG2 , 

0QG2Y ,OSTRl . 0STR1C, 0CTR1M. OSTR2 . 0STR2Y . OOCV1 , 

0UGV2 . OU5V2Y, PCI ,QEF1, SIP ft SCRATCH 
TYPB DB GB0K1B. PTBLKM. KJJ ft POP POST-O • DBC 
TYPE DB PJC S CYCLIC STATICS - POP CPFDR 
TYPE PARM. . I.N, CP, ZCYCLIC, CEOIRJ. LOAOO, POSTU ft 
TYPB PARM. , X . Y. D6CD1AG, RMSTAAZN $ 

TYPE PARM, .CHAR8, H, DBCPROG ft 

TYPE PARM. , CKAAS, Y. DBCOVWRT, WCF. ACOUT* PBAX' ft 
TYPE PARM. .RS.Y.PRBPDB-l. $ 

TYPE PARM. NDCL, X. N. ACOUSTIC ft ACOUSTIC ELEMENT FLAG - SET IN SBPABP2 
TYPB PARM. MCL. I . N. DBCPATH S DUMMY PARAMETER TO PASS QUALIFIERS 
TYPE DB.MLM, KDICT, LAMA ft FOP POST--3 
S L OCA L PARAMETERS 

TYPB PARM, . I. N. PFZLI.SI ZD. CARDMO, RICORD . MH. MOXOUT- -1 . DBS I TER $ 

TYPB PARM. .CHARE. N, APP, AFP 1.APP3 i 

TYPB PARM. , LOGICAL i N. SOAT2 .STATICS . GPFDR, fR-HBAT , STATMH. AERO. F8 ft 
| USER PARAMETERS 

TYPB PARM. , I. Y, POST, GRDPOT.HOBLOP.MOBLOP. CURVPLOT. SIAM. SIAG.LSTRM f 
TYPE PARM. , I , Y . CURV.NUMOUT1 . NUN0UT2 . NUMOUT. SRTOPT. SRTICTYP t 
TYPE PARM, , I.Y.MOMSCSTR. SI , SIM, SIC, OUTOPT. OG, NXNTPTE . DOPT. PDRMSG 5 
TYPB PARM. , I. Y. REPRINT, TABID i 

TYPE PARM. ,RS. Y.BZGIRl. BIGBR2. BIGER. TINY, TABS S 

$ USER PARAMETERS A PPL I ID TO ALL EE AMD RESOLVED IM SlEERl 

TYPB PARM. .I.N.ENRSL $ 

ft USER PARAMETERS RESOLVED IN SUPER) 

TYPE PARM, . I.N, SCRSPBC, RSPECTRA , HOCOMPS S 
TYPE PARM. .CB.M, ALPHA- (10. 0 0) ft 
TYPE FARM. .CS.N.ALPKA1* (1.0. 0.0) ft 
TYPB PARM. .RS.N.RWOISl. PREQ1 $ 

TYPE PARM. .RJ.N.ll. 12. IS $ 

TYPB PARM. ,RS,Y, BETA S 

TYPB PABM, . I , Y . CONI , LMODBS . I $ 

f **+,»**«»*»*,*•••*•«.«**»*•• *•••**»»*•«» 

$ owrwn CONTROLS 

* 

s THE Otmvn DMkP INSTRUCTIONS IN THIS SUBONAP HAVE BEEN SUPPLIED 
| AS A COURTESY TO POA/PATRAN AM SDRC/I-DBAS USERS AND QUEETIOMS 
$ AS TO HOW IT IKTERFACIS WITH THE PATRAM OR I- DBAS PROGRAM SHOULD 
ft BE DIRECTED TOi 
ft 

ft PATRAM USERS CONTACT. X-DEAS USERS CONTACT, 

ft 

( PDA ENGINES RING SDRC 

S 2*7ft RICH ILL AVI. 2000 EASTMAN DR. 

| COSTA MESA. CA 9 2ft 2* MILFORD. OH 4S1S0 

ft 

I PARAM. POST. -1 , OUTPUTS THE APPROPXATB PILES FOR 
% THE PM/PATRAM MASPAT PROGRAM VERSXCM 2.0. 

< 

ft PARAM. POST, -2. OUTPUTS THE APPROPIATE FILES FOR 
s THE EDRC/I-DSAS DATA LOADER PROGRAM VERSION ).Q. 

S 

TYPE PARM, .CHARE, Y, OQC, OUG.OEP .OBS . OSE. OCNP, OOPS . OESt $ 

TYPE PARM, , CHAR* , Y , CUMU , OCPP , OUGCORD I 
TYPE PARM. , I , Y , 0UNZT2 , OMAXR ft 
TYPE PARM. , I.N. CTAPE2 t 
PILE PEDP-SAVE.OWBBfRT $ 

*•♦♦♦♦♦•♦♦♦*♦•*♦******•*•♦*♦♦♦♦+*•♦ 
STATICS* (APP* ' STATICS ' OR AFP-MST' OR 

(APP* 'CYC ‘ AMD APP1.’ STATICS 1 M t 

* 

PARAML SILS// 'TRAILER' /l/S.M.NSILS $ 

PARAKL S ILS//' TRAILER ' /2/S.H, LUSETS $ 

IP ( NOT (STATICS OR KLKEAT) AMD (NBILS'ft) <»LUSETS ) THEN t 
PLTTRAN BCPCTE, SILS/BGPDP. SIP/LUSBTS/8,M,LUSEP I 

ELSE t 

BQUXVX BGPDT8/ BCPDP /ALWAYS $ 

BQMZVE SILS /SIP/ ALWAYS S 

L US BF -LUSETS % 

EM) IP I 

I 

IP | APP- 'STATICS ' AM) NOT(NUflAT) OR 

(APP- ‘ASIC ' AM) APPlo'BELI ' > ) THEN t 
VBCPLOT QGE.BOPDP.aOBXINS.CrrNB.CASSDR./gGSUM/ORDPNr/O/l/ 
'EPCPOMCE' S 

VBCPLOT QCS.BGPDP, SQEXXMS. CSTW.CASEDR, /QGB//0/S/ ' MAXDSJH' / 

’ SPCPORCE' / ’ S' I 

VBCPLOT U0V1, BCPDP. BQCXIMB . CSTMS. CASBDR. AXR/B/ /0/5/ 'NAXIMM' / 

’ DISPLACE ' / ' MBMTS ' I 


IP < APP* STATICS ' AND NOT (NLMtAT) ) VBCPLOT. 

PJI .BCPDP. BQEX IMS. CSTMS. CASED*, /PJB//0/5/ MAXIMUM'/ 
APPLIED' /'LOADS' $ 

END IF f 

t 

IF ( SCRSFBC*- 1 ) RETURN $ 

( 

EQUIVX UCVS/ZU2RI9 /ALWAYS $ 

EQUIVX UGVS/ZUIA20/ ALWAYS $ 

MATMDD UOVS /ZUZA18./1/1 $ 

DELETE /ZUIR1S, ZUIR20 , , . ft 
ZS-0.0 ft 
1*1 ft 

ft EXTRACT THE N-TH NODE ( CONTROLLED BY VARIABLE I| 
ft SO WE CAM HAVE THE S BPS RATE RMS VALUES FOR EACH MODE 
ft TOTAL RMS DISPLACEMENT IS IN ZS 
DO WHILE ( X< -LMODBS ) ft 

MESSAGE //'THIS IS MODE • '/I ft 
PARAML UGVS// 'TRAILER'/ 1/S. N.NOCOLS 5 NOCOLS'S IN UGVS 
//'LMODBS* ' /LMODBS/ 'no. of col»n*. ' /NOCOLS ft 
KATGSN ./C/1 /NOCOLS ft 

MATNOD C /Cl, /I /I ft 

NATCSN , /A/7/1 /NOCOLS $ 

MATMOD a /B, /12/S.N, NULLS/V, Y. NIM*1 ft 

ADD B.C1/D/U. 0,0. 01/1-1.0,0.0) ft 

MATNOD UCVS /UCVSXYI, /l/I ft 

MERCS UCVSXYI D./UCVSX/l ft 

d*l*t* /ugvs, , , . ft 
•<Rjivx ugv«x /ugvi /always ft 

ft 

SDR2 CASBDR. CSTMS. NPTS . DIT. 8QBXIN6 , . BTT , OLB1 . BCPDP. 

PJl.QCS.UGVS.BST , XYCDBDR/ 

OPG1. OQG1 , 0UGV1, OES1, OBF1 , PUGV/APP1/S. N, NOSORT2 / 

NOCOMPS// ///ACOUT/PREPDB ft SORT! O LOADS. S PC FORCES, 
ft NOCOMPS //ACOUT/PREFDB ft SORT1 OLOADfi. S PC FORCES 

ft DISPLACEMENTS. ELEMENT STRESSES AND FORCES, PLOT VECTORS 

SDR2 CASBDR. CSTMS, MFTS. DIT. BQSXINS . , BTT , OLB 1 , BCPDP , 

PJI , QGS , UGVS , 1ST . 

/ . , , OSTR1 . . / 

APP1/S, N.NOSTR2/3//// /ACOUT/PREFDB ft SORT1 ELEMENT STRAINS 
ft APP1/S, N. NOSTR2/1//ACOUT/PRSFDB ft SORT1 ELEMENT STRAINS 

IF | APP* FREQRBSP- AND NOT (AERO) AND ACOUSTICS ) SDR2. 

CASBDR, CSTMS , MFTS , . BQEXZNS, ,. OLB1 , BCPDP, , UGVS. 1ST, / 

, .OUG1F. , , / 

APP// /ACOUSTIC// //ACOUT/PREFDB ft SORT1 ACOUSTIC PRESSURE OUTPUT 
ft APP// /ACOUSTIC /ACOUT/PREFDB ft SORT! ACOUSTIC PRESSURE OUTPUT 
EQUIVX OPCl/gpGl /ALWAYS ft 
EQUIVX OQG1/QQC1/ ALWAYS ft 
BQUXVX OUCV1/QUCV1/ ALWAYS ft 
EQUIVX OSTR1/QCTR1/ ALWAYS ft 
f 

IP ( MpCF- YES' AND STATICS ) THEN ft 
ADOS PJ1.QGS. , . /PQ ft 

MPT AD KJJ.UOVS, PQ/QOM///-1 ft 

SDR2 CASBDR. CSTH8. MFTS. DIT, EQBXINS . . BTT.OLB1 .BGPDP. 

, QGM, , BST. / 

, OQOf /APP1/S.N.NP /NOCOMPS ft 

MESSAGE / / ' ' / 

'MULTIPOINT PORCIS OP '/ 
■CONSTRAINT' ft 
OPP OQGMV / ft 

SM>IP ft 

s 

ft OLD MAP -- IP FIXEDB < 0 THEN STORE PUGV?? 

ft EQUIV PUGV WILL PULL PATH TO PUCVS WITH PATH SBC BOM 

ft 

IP | POST-O AND OQG* ‘ YES ' ) THEN ft 
ft •• OUTPUT2 OQC1//OTAPI2/OUMXT3 //OMAXR ft 
OTXPE2-0 ft 
DO IP ft 
ft 

IP ( POST* -2 AMD STATICS Al® NOT (AERO) AM) 

GETSY6(Mf. Sf |«>0 ) SDR2. 

CASBDR, CSTMS , MFTS , DIT , EQBXINS . . ETT.OLBl . BCPDP. 

PJI, QOS. UGVS. SET .XYCDBDR/ 

. .TOUCVt. . ./APP1/S.N.NOSORT2 /NOCOMPS ft 

ft 

IP ( POST-O AND OUG*‘ YES* ) THEN ft 

IP ( PQST--2 AMI STATICS AM NOT (AERO) AM 
GSTSYS ( Ml . 5S ) - >0 ) THEN ft 
ft** OUTPUTS TOOCV l / /OTAPE2 /OUMIT2 / /OMAXR ft 

ELBE IP < OUGCORD* BASIC' ) THEN ft 

VBCPLOT UGVS, BGPDP. SQEXZME, CSTMS, CASBDR, /UWSB//0/1 ft 
IP ( POST* - 1 ) THEM ft 

8DR2 CABSDR, CSTMS .EQBXINS. , . OLB1 . BOPDTS , . . UOVSB . ./ 

. , OUGV1PAT. . , ZAP FI ft 

ft** OUTPUTS OUCV1PAT/ /OTAPE2 /OUMIT2 / /OMAXR ft 

BUB IP ( POST* *2 1 THEN ft 
IP ( APP1**REZG- > THEM ft 

SDR2 CASBDR, CETM, . .BQVXZMfl. . , OLB1 . BOPDTS , . .UOVSB. ,/ 

. , BOPHIG. . ,/APPl ft 

ft** outrun BOPHIG/ /OTAPE2/OUMIT3 //OMAXR ft 

BUB ft 

8DR2 CASBDR . CSTMS . . .BQCXZMS, . ,OLBl. BGPDTS. , . USVSB, , / 
, BOUCV1. . ,/APPl ft 

ft** OUTFUT2 BOU3V1/ /OTAFE2/OUMIT2// OMAXR ft 

EM IP ft 
SMZP ft 
BUB ft 

ft** OUTPUT2 OUGV1//GTAPB2/OUNIT2//OMUCR ft 

EM IP ft 
OTAPB3*0 ft 


Cl 



ENDIF $ 

$ 

IF { PO£T<0 AND OBS ■ ' YES ' ) THEN $ 

OVTWT2 OBS1 / /OTAPB2 /OLNIT2 //OMAXR S 
OTXPE2-0 $ 

ENDIF t 
$ 

IF ( < (STATICS AND NOT CAIRO) ) OR NLHEAT) AND 
(CBTSYS (NH. S6> <»0) AND SEID-0 | THEN $ 

sown* sils.uset . ocvs.owi.slt.sst . dit. Q ca. olt. /howi / 

TAES/O $ 

DSL STS /OEF1 .... | 

BQUIVX NOEFl/OBFl /ALWAYS ( 

BQUIVX H0EF1/QW l /ALWAYS $ 

SM)ZP $ 

* 

ZP { OW* ' YSS ' ) THEN $ 

IP ( POST--2 AND STATICS AND NOT I A KAO) AND 

crrsYStm. S6j<x> j tkxn $ 

OUTPUTS HOB FI / /CTAPB2/OUNIT2/ /OKAXR $ 

OTA Pt 2.0 $ 

BLSB IP < POST«0 ) THEN S 
$*• OUTPUTS OSPl//OTAPS2/OOKm//CHAXR $ 

OTAPB2-0 f 


KNDZP $ 

S 

IP | PC6T<0 AND OBS- ' YSS ' ) THIN $ 

5** OUT PUTS OSTR1//OTAPB2/OUNIT2//OMAXR $ 

OTAPX2.0 i 
SNDIP $ 

S 

IP ( APP-NLST AND SiID.0 I THEN $ 

KBRGIOFP OW1, OSPNLXX/OSriX J 
HBRGBOPP OSS1. OBSNLXX/OBS1X $ 

BLSI S 

IT < STATICS OR APP-'RSIC ’ | THEN S 

SWO( CASKDR, 0«P1 , OBEl . GBCH2S. CE0K3S, 8ST. CSTHS. KPTS . DIT/ 
0EP1X.0IS1X/S. N. NOXOUT $ 

ELSE IP < APPlo mMIC’ > THEN $ 

SDRXD CASEDR. OEF 1 . OES1 . C10M2S . GE0IO6 . *ST. CSTHS. 

KPTS, DIT, UGVS, DLT. OLE1/ 

0EF1X.08S1X/S.N. NOXOUT/ APP/COUFKASS ( 

SNDIP $ 

BQUIVX OEP 1/OE FIX /NOXOUT $ 

BQUIVX OBS 1/OSS1X /NOXOUT $ 

BNDIF $ APP-'NLST' AND SEID-Q 
$ 

BQUIVX OBFIX/QEFIX/ ALWAYS $ 

BQUIVX OBSIX/QBS1X/ ALWAYS $ 

$ 

STATNH* (STATICS AND NOT(NUtBAT)) $ 

IP ( APP1< > ' IMRIIG' AND (STAT»f OR 

(APP-'REIG* AND APPlo'WLl I OR APP* TRANAESP' | ) THEN $ 
PARAKL POETCDB/ / ' PRESENCE' ////S.N.NOPOST i 
IP ( NOPOST » -1 | THEN $ $ CAID POINT STRESS FACTORS 

PLTSBT POS7C » , XQEX INS . ECTS / PLTP PLTPARP . CPSBTP . BLfi ETP / 

S , N.NSXLS $ 

CPSTRi POSTC DB.BGPOTS.ICTS, CSTHS. BLSITP UTS. BQBXINS/ECPSF 

S BGPEF - ILEHBNT TO CRID POINT SURFACE FACTORS 

SNDIP $ CRID POINT STRESS FACTORS 
PARAKL SGPSF// ’ PRESENCE ' / / / /S, N. NOBGPSF $ 

PARAML OBS1X/ / ‘ PRXSBNCE '////£, N. NOOSS1 $ 

IP ( NOBGPSP>- 1 AND NOOESl>-l | THEN $ 

$ SORT1 CRID POINT STRESSES 

GPSTR2 CASEDR. EC PSP . CPUS, OBSIX/OCSl, BGPSTR///APP $ 

BQUIVX OGS1/QCS1 /ALWAYS $ 

$ 

S KSSH ERROR ANALYSIS - LINEAR STATIC ANALYSIS ONLY 

PARAKL ECPSTR// ’ PRESENCE' /// /S.H, NOECPSTR $ 

IP ( N01GPSTR>- 1 AND STATICS AND NOT (APP* NLST ) ) STDCON, 
CASEDR. BCPSP . SQCXINS. OBS IX. ECPSTR. BCTS/ 

OBDSl . OCDS1 . ELDCT. CPDCT/ 

S . N. NOSD61/S . N. N0GD61/6 . N, NOE DTI /S.N, NOCDTl /APP $ 

$ 


IP ( POST- 1 AND OCRS- YES' ) THEN $ 

$*■ OUTPUT! 0061/ /OTAPB2/OUNIT2 / /OMAXR $ 

OTAPB2-0 t 
ENDIP i 

* 

SNDIP $ NOBCPSP>- 1 AND N00E6I>-1 


s 


« 

SORT! -NOT ( ANH> ( NOSORT3 , M0STR2 ) ) $ 

I IP FLUID /STRUCTURE MODEL AND FREQUENCY RESPONSE. THEN 
S DO SORTS EVEN THOUGH DDAJM-- 1 . 

S BUT DO NOT PBRPOJM MATRIX METHOD DATA RECOVERY (OORMfJ 
IP ( APP1* MMREIC ' OR 

(CURVPLCT*-! AND ISORT2 OR I PS AND APP* PR IQ RES P | | ) ) THEN C S0RT2 
SDR 3 OUCVI.OPC1,OQGI.O«F1X.O«S1X.OSTR1/ 

OUCV2.OPC2.OQC2.OBF 2. 0BS2.0STR2 5 SORTS OUTPUT 

SDR3 0UC1P /OUC2P f 

ZP ( STATNH OR APP*'RSIC * OR 

(APP- ' TRANRISP * AND NOT (NLMBAT) ) | SDR3. 

0081. 00061, OBDSl. , . /0CS2.0CDS2 ,0*062 , . . t 
IT ( APP 1* NHRI IQ ’ > THEN ( 

IP (SEZDaO) THEN $ 

DELETE /0PC1.0PC2. . . / t FROM SCR 2 ABOVE 

SCR 2 CAEBDR. , , , 8QIXINS. , ,OLBM. . PCD XYCDEDR/ 

0FG1 /APP $ 

snu opgi /ores / t 

BNDIF S 

SCRIM CASEDR. 1JHVP. OUM, OUGV2, OQC2 ,OB82 .OBF2. XYCOEOR/ 

OOGV2Y, OQC2Y.OBS2Y, 0SP2Y, / S 
BQUIVX OUSV2Y/OUCV3 /ALWAYS S 
BQUIVX 0QC2Y /0QG2 /ALWAYS < 

BQUIVX OI62Y/OB6 2 /ALWAYS $ 


BQUIVX 0EP2Y/0BP2 /ALWAYS $ 

DDRNH CASEDR, UHVF. OLBH. . , OSTR2 . .XYCDEDR/, .0STR2Y. .//3 $ 

BQUIVX OCTR2 Y/06TR2 /ALWAYS $ 

KNDIF S APP1* 'HMRBIC 
S 

BQUIVX OPG2/QPG2/ ALWAYS S 

BQUIVX OQC2/QQG2 /ALWAYS $ 

BQUIVX 0UGV2 /QUGV2/ALWAYS $ 

BQUIVX OBS2/QBS2/ALWAYS $ 

BQUIVX 0BP2/QBP 2 /ALWAYS S 

BQUIVX 0STR2 /QSTR2 /ALWAYS $ 

BQUIVX OGS2/QCS2/ALWAYS $ 

$ 

IP ( RMSTRAXN-2 } THIN S 

OPP 0UCV2 , 0PC2 , 0QG2 . 0SP2 . 0BS2 , 0STR2 / /S , N, CARDNO $ 

OPP 0GDS2 , 0BDS2 . 0UC2F/ /S , N. CARDNO $ 

ENDIF $ 

$ 

IP < STATNH OR APP-'REIC ' OR 

(APP- 'TRANRBSP' AND NCT INLH8AT) ) J OPP OC62/ /S . N. CARDNO $ 
XYTAAN XYCDeDR.OPG2.OQC2. 0UGV2 . 0BS2 , 0BP2 /XYPLTT/APP/ ' PSBT' / 

S.N, PPILB/S. N, CARDNO /S , N. NOXYP $ 

IP ( NOXYP »0 ) XYPLOT XYPLTT// $ 

J 

IP ( RSPSCTRA>*0 ) TKBN $ 

$ SKIP IP NONLINEAR ANALYSIS 
RECORD* 0 $ INITIALIZE 
DO WHILE ( RECORD* »- 1 | $ 

RSPBC PRL . 0UGV2 , SPSBL/OXRBSP/S. N, RECORD $ 

IP ( ABCORD»-0 J THEN S 

IP ( RSPRINT»-0 ) OPP OXRBSP/ /S.N, CARDNO S 

XYTRAN XYCDEDR , OXRBSP /XYPLTSS/ ' RSPBC* /' PSBT' / 

S.N, PPILB/S. N. CARDNO/S , N, NOXYPLT /TABID S 
IP ( NOXYPLT*-(J ) XYPLOT XYPLTSS// 5 
BNDIF $ 

END DO S RECORD* »-l 
S 

BLSI IP { APP * -FRBQRiSP' ) TH8N S RSPECTRAv-0 
DPD DYNAMICS. GPLS, SILS . USBT, . / 

XGPL. XSIL, XUSET. , . HPSDL XBQDYN/ 

-1/DUM2/0/0/S.N.NDPSDL/0/0/ 0 / 1 /123/DUKia S 

S DYNSTAT/NBWDYN 

IP ( NOPSDL > *1 | THEN 
$ DELBTB /PS DP .AUTO. , . 5 

IF I RMS TRAIN * 1) THEN 
$ i«r.pm hPSOL.O IP) // J 

$ Htprii XYCDBOR. HPSDL. 0UGV2 . 0IP2 , CASBDR // S 

RANDOM XYCDEDR. DIT. HPSDL, OUGV2 , 0PG2 , 0QG2 . 06TR2 ,OBP2. CASEDR/ 

PSDP. AUTO/ S.N, NORAND t 
PARAML 0STR2// ' PRESENCE 1 / / / /6, N, NOOSTA2 S 
KBSSACB //’OSTR2 (09-14-921- ' /NOOSTR2 S 
S 0UTFUT2 , . , , //- 1/12 / LABEL* PDALABBL J 
i 0UTPUT2 OSTR1//0/12 $ 

S 0UTPUT2 , . , .//-E/12 J 
ELSE t 

RANDOM XYCDBDR. DIT. HPSDL. 0UGV2 , 0PG2, 0QC2, 0ES2 . OEP2 . CASBDR/ 
PSDP. AUTO /S.N. NORAND $ 

BNDIF $ 

IF ( NORAND > -I > THEN C 

XYTRAN XYCDEDR. PSDP, AUTO, .. / 

XYPLTRS/ 

RAND/ 'PSBT'/S.N, PPILB/S . N. CARDNO/S , N . NOXYPLT J 
IP < NOXYPLT ► -1 ) XYPLOT XYPLTRS// $ 

SNDIP $ 

ENDIF % 

BNDIF S ELSE IF APP-PREQRESP 
t 

IF ( P06T-0 | THEN $ 

DEC 0PG2 . CUGV2 , 0BP2 . 0SS2 , 0QG2 , GPLS , HGPSTR, BGPSP.GPDCT, 

HLDCT // 

‘ OPG ' / ‘ OUG ' / ' OBF ' / (MS ' / OQC’ / ‘ GPL ’ / 'GPS * / ' SVP* / CPDCT'/ 
' ELDCT* /////////// 

- 1/DECPATH/S .H.CP/APPl/ICYCLIC/GKJICI/LOADU/POSTU/ 
DBCDIAC/DECPROC/DECOVWRT/DSsrTER S 

DEC 0STR2 // 

’OB S'//////////////////// 

- 1/DRCPATH/S.H. CP/APP1/ICYCLIC/CBOMJ/LOADU/POSTU/ 

D ECDIAG/DBCP ROC/ DEC OVWRT/ DBS ITER $ 

DEC OUC2P // 

•OUC' /////// ///////////// 

- 1 / DECPATH /S . N . CP/ APPI / ICYCLIC /CEOMU/LOADU/POSTU/ 

DBCDIAC ' DECP ROC/ DEC OVWRT/ DBS ITER S 

ENDIP $ POST-O 
I 

i SCALED SPECTRA RESPONSE NOT INCLUDED HERE 

S 

ELSE $ APPl-'MMRBIG ' OR ( (CUAVPLOT*-l OR PS) AND S0RT2 ) 

OPP 0UGV1. OPGI. OQCl.OBPlX.OSTRl/ /S.N, CARDNO ( 

OPP 0GDS1 , OBDSl, 0UG1P//S, N, CARDNO $ 

IP ( NOT (STATICS OR APP.’RBIG • ) ) OPP OBS IX / / 6 , N , CARINO $ 

IP < STATICS OR APP-REIG ■ OR 

APP* ' TRAN RES P' ) OPP OCS 1 //£. N , CARDNO $ 

ZP ( STATICS OR APP-'REIG * > THEN $ 

IP < NOCOMPS 0 ) THEN $ NOCOMPS »« 0 COMPOSITE PLY STRESS 
SDR 2 CASEDR. CSTHS. KPTS. DIT. IQIXINS . . BTT, 0LB1 , BGPDP, 

PJ1.QGS. UCVS , 1ST. XYCDEDR/ 

. . .OES1A. . /'STATICS '/S.N, NOS/2 t 
SDR COMP CASBDR. MPTS, BPTS. ETT, 1ST, 0CS1A, OSP1, DIT/ 

ES1C.EPIT, /LSTRN $ 

STRSORT ES1C, /OSSIC/NUMOUTI /EIGER 1 | 

ETRSOItr EPIT, /0EPIT/HUM0UT2/EICER2 I 
OPP OEE1C.OBPIT. , . //S. N, CARDNO t 

IP ( POST* -2 AK> OOCP- YES ' > THEN $ 

HATHOD 0ES1C /0STR1C. /13 $ 

OUTPUTS 0SS1C, OBPIT/ /OTAFE2/OCMIT2/ /OMAXR $ 

O1APE2-0 S 
ENDIP ( 


C2 


f 


t 


I 


I 

I 

I 


< 

< 

( 


$ 


m>ir I NOCCNFS ** 0 COMPOSITE PLY STRESS 
PQUIVX OH1X/OU11U /SI | 

xp ( ix >. o j strport ossix, iwta/opsixi/nuhout/bigsr/ 

SKTOFT/EKT1LTYP 5 ELEMENT STRESS SORTING 
OFF 0*61X1/ /S, H. CARDNO | t PRINT KLEMENT S1USSU 

PARAML FORCE// PRESENCE '////*, H.NOPORC* t 

IP ( NCMSGSTR >■ 0 AM> MOFORCS >• 0 > KSGSTAES FORCE . 0BS1X/ 

/S.N.PFILf/NOKSCSTR { 

IP { CURV >■ 0 ) THEN | 8TR1SS/STRAIN TRANS FORMATION TO GRID POINTS 
OR KATL COORD SYSTEM 

ID RF24D40 IN I* ). S FOR PARAMETER EXPLANATION 
CURV OSS IX , HPTS , CSTNB , 1ST .SILS.CPLS/OESIM.OSSIG/OUTOPT/ 

OC/HINTPTS S 

BQUIVX OBJ 1M/Q1J IK/ ALVA Y£ $ 

BQUXVX 0BS1C/QBS1G/ ALWAYS t 

XP { SIM » 0 } THEN $ 

STMORT 0CS1M. IMYTA/0SS1K2 /NUKX7T /BXGRR/SRTOFT/SRTELTYP $ 

OFF OBS1H1/ /£, N, CARDNO t PRINT STRESS** IN NATL COORD SYS 

DO IF t 

xp i sio » o i the n $ 

STRSORT onto. INDTA/OEClGl/NUNOUT/BlGtF/SRTOPT/SRTELTYP J 
OFF OBBLGl//S,M. CARDNO t PRINT STRESSES AT GRID POINTS 

MR) IF | 

DIAOON<23) S 

CURV OSTA1 .KPTS.CSTNJ, 1ST , SIL1 , GPLS/OSTR1M. OfiTRlO/OOTOPT/ 
OO/NXNTFTt f STRAIN! 

D1AOOPF (23 ) I 

PQUIVX OCTIUM/qfTRlM/ADAYS S 

BQUIVX 0CTR1G/Q6TR1C/ ALWAYS $ 

IP ( 11AM >• 0 ) OFF OBTAIN/ /£ , N. CARDNO $ PRINT STRAINS IN 

KATL COORD SYS 

IP ( S1AC >• 0 ) OPT 0STR1C//S.N, CARDNO 5 PRINT STRAINS AT 

GRID POINTS 

BM>ZP $ STRESS / STRA IN TRANSFORMATION TO GRID POINTS OR NATL COORD SY 

PARAML XTCDB// ’ PRESENCE //// 8, N, N0XYCD6 X 
IP ( MOXYCDR >> 0 } THIN ( 

CURVPLOT SQBXIN6 , DC POTS . 0L11 . XYCDBDR , OPG1 , 0QC1 , OUGV1 , OES 1G. / 
OPC2X.OQC2X.OUC3X.OBS2X. /DOPT $ 

XYTRAN XTCDRDR, OPGiX, OQG2X . OOC2X. OES2X. /XYPITS/'SETI' /PSIT/ 
S,N, PFIL1/S. N. CARDNO/ S. M.NOXYP $ 

XP ( WOXYP >« 0 ) XYPLOT XYPLTS/ / $ 

END IF I CUR VP LOT 

IP < OFT DR 1 THIN I GRID POINT FORCE 

XP ( APT* 'CYC' AND APP1. 1 STATICS ' | THEN | 

PQUIVX PJC /PCI /ALWAYS $ 

ELSE ( 

BQUXVX PJl/PGl/ALbAY* f 


drcdiac/dpcproc/dpcowat/ dbsiter i 

DBC OES 1C 

*ORSC '//////////// //////// 

- 1/DPCFATH/l . N . CP/ APP1 / ICYCLIC /OPOHU/LOADU/ POSTU / 
DPCDIAC /DBCPROO/ DRCCVWRT/ DBS ITER $ 

DPC OUOIP // 

*CUG' ///////////// /////// 

- 1/DPCPATH/S , N.CP/APPl/ XCYCLIC/GPOMU/LOADU/FOCTU/ 
DBCDIAG/DBCPAOG/ DPCOVWRT/ DCS ITER S 

END IF $ 

$ 

BNDXF S APP1* 'MMR1XG ' OR (CURVPLOT- -1 AND £0RT2) 

$ 

PV*1 padf //'dtl'/l/fl/in, ra.di.l 5 
ZS.tS.RMBDISl i 

•••Mfl* //’ cwadial- ‘ /raadial/'ZS- ■ /*S $ 
p*df//'dti'/l/5//»,n,nod*no 5 
•••MR //'ikkWS- ' /nodano $ 

NODSNO-NODBNO/ 10 $ 

CCN1-NO0BNO* 6 - > $ 

non* UGVSXYI/UCVSNORM// /a. n, ugviMX S 
P«TMil ugvanora/ / '<A»i ' /1/CONl/a. n. nodwtlia $ 

•«“N //’nwkw co<i 1 nodadis- * /nodano /coni /nodadia { 
XP ( WOOED IS-O. 000 ) THEN $ 

MESSAGE 1 !' FATAL MESSAGE FROM THE MODIFIED MAP' / 

' THE GRID POINT ' /NOOENO/ ' HAS ZERO RMS DIS . . ' / 

' PLEASE CHOOSE ANOTHER POINT IN THE ‘XYPLOT* 1 / 

' OR 'XY PRINT’ . * $ 

EXIT $ 

END IF $ 

I raadi a l.raadial /nodadia $ 

raadial.rwadial/aba (nodadia) $ 
paraal pad*// 'dti' /1/9/a, n, Iraqi $ 

//■ fraquancy. '/Iraqi % 

ALPHAl-CKPLX (-FRBQ1.0 .0) 5 
aaaaaga //ra.dial* ' /raadial X 
ALPHA -CMFLXIRMSD IS 1,0.0) S 

DELETE /P $ 

ADO UCVSNORM, ZUZR1B/P /ALPHA S 

BQUIVX P/ZUZR1R/ ALWAYS S 

NORM UGVSX/UGVSXNORM/ / / * ,n , ugvaxaax % 

DELETE /Z $ 

ADO OCVSXNORM. ZUZR20 / Z /ALPHA S 

DELETE /ZUZR20 $ 

BQUIVX Z/ZUZR20/ ALWAYS S 
I-I.l $ 

BQUIVX ZUZR1 9 /UCVS/ ALWAYS 5 
ENDDO t loLMOOBS $ 
t RMSDIS1 *Z& f 

MESSAGE //'THE RMS DIS. AT POINT 1 /NODBNO/ ' IS ' /ZS/ 

' WITH TOTAL OP ' /LMODBS/ ' MODES' $ 

DELETE /MUGNX $ 

•<RiivK zucrlt /ougni /«lw«ya $ 


APF2«APF1 { 

XP ( APP. NLST' ) APF2* ' STATICS * f 

GPFDR CASED*. UCVS. KXLM. EDICT, PCTS . 1QEXINS. GFECT. PCI , QGS, 
PGFDTS.BILE.CnNS, VBLM.PTILEN/ 

ONRCYl. OGPFB1/APP2 /TINT I 
BQUXVX 0NRGY1/QNRGY1 /ALWAYS $ 

BQUIVX OOPPB1/QGPFB1/ ALWAYS { 

OFF 0MUZY1.0GPPP1//S.N. CARDNO t 

$ 

XP ( POrr<0 ) THEN $ 

IP < FOST--1 AMD OCPP. YES’ ) THEN $ 

OUTFUT2 OGPPE1/ /0TAPE2 /OUNIT2 / /OMAXR | 

OTA P* 2.0 { 

END IF X 

IP ( OESt-YES ) THEN f 

OOmm MRCYl/ /OTAPE2 /0<JNIT2 f /OMAXR S 
OTAPE2-0 $ 
l»tP f 

XP 4 POTT- -2 AMD OUMU.'YBS’ AND APP.’UIG' ) THEN t 

OFF DR CA1SDR. U0V1.MBLM.MDICT, PCTS , BQBXINS, GFBCT, FG1 , QGS, 
BCFDTS, EILS.CSTNE. VELIM. / 

0JPWY2.0GPFP2/AFF1/TXHY S 
DPVI1W XLAMAbLAKA (MIUIE WILDCARD* TRUE) $ 

ovmrra xlaja . omwy 2 / /otafe 2 /ouwtt 2 / /omax* t 
OTRFP2-0 | 

M>XF | 

PDIF f 

XP ( NOPLOP >• tOI NO* LOP >• 0 ) THEN S 

ELF DR OGF FBI. GFECT . CETM*. SILS. CFLS . PGFDTS/OELOP1 . OELOF1/ 

NOPLOF/MOBLOF f 

OFF OBLOP1, O* LOT 1//S.N, CARDNO | 

«*IF S ELEMENT ALI04ED GRID POINT FORCE 
PWXP S GRID POINT FORCE 
END IF S STATICS OR APP. * RSIC 

I 

IP 4 rorr.0 ) THEN t 

DPC 0PC1 . 0UGV1 . 0BP1X , OBS1X . OQG1 . GFLS . ICFSTR , BGPSF . GPOCT, ELCCT. 
Oip xm ■ OGFFP1 // 

*0PC' / ■ OUO * / 'QfP ' / 'OPS 1 / ' OQC • / 'GPL ' / 'CPS ■ / • svr ' / 

'GFDCT' / * SLDCT' / 'OBSS ' / 'OCPP '///////// 

- 1 / DPCFATH /t , N , CP/ AFF 1 / ICYCL IC /OPOMU /LOADU/ POSTU / 

DPCDIAC /DPCFROa/ DPCOVWRT/ DPSITVR X 

DPC OPTR1 // 

OPE'//////////////////// 

•1/DPCFATH/l , N, CF/AFP1/ ICYCUC/CPOMU /LOADU /POSTU/ 


S Thia ia Cor i«p«rat« »od»a 
t 

S o^uivv zuxr20/Rugni/nlwitya $ 

I 

t FIND LARGB1T P08ITXVB AND SMALLEST ( NEGATIVE) TERMS IN MATRICES 
$ 

I INPUT MATRIX IS ZUZR03 

% 

IQUIVX MUGNX /ZUZR03 /ALWAYS $ 

DIAGONAL ZUZR03 / APS / * MfOtLB* / 1 . 5 

ADO APS, ZUZR03/FO1XTIVS/ ( .S . 0 . ) / ( . S , 0 . ) ( AEMCV1 NEGATIVES 

MATMOD POPITIVI /MAX, ft $ MAXIMUM VALUES 

ADO APS. ZUZR03 /NEGATIVE/ I . S , 0 .)/(•. S , Q . ) S REMOVE POSITIVES 

MATMOD NSCATXVI /MIN, (1 $ MINIMUM 

ADO MIN, /KIWI/ (-1. .0.) t 

norm MAX/out»*x///a.n. IocxUax S 

norm MINN/ouuiin/ //*, n. localain $ 

// ' local* ax. ‘/IocmIrmx / 1 local *in. /loon Lain $ 

IF 4 LOCALMIN » LOCALMAXI THEM $ 

ADO ZUZR03. /ZUZR04/4-1. 0,0.0) S 
ELSE S 

ADO ZUZR03. /ZUZR04/ (*1.0, (1.0) $ 

•Norp t 

t 

PQUIVX ZU1R04/MMNX /ALWAYS t 

DELETE /ZUZR04 f 

non. ZUlRli/DUMOf ///a.n. Madia t 
MESSAGE //‘THE KM RMS DIS. XS ’ /MAX DIS $ 

( 

PAPAML PCDRDR// ' PRSSSNCB' /// /S.N. JPLOT ( IS THERE AM OVTFUT(FLOT) 

< FOR THIS SUPBRELWENT f 

IF 4 JPLOT >-0 ) THEN $ PHASE III DSPORM0 AM) CCMTOUR PLOTTING 
PLTSST PCDRDR. PQBXIME , ICTS /PLTXX , PLTFAJlX , GFS1TSX , ELSPTSX/ 
*.M,N*ILS/8.N, JFLOT J 
PRTMSG PLTXX //PORM6G $ 

IF ( JFLOT »«0 ) THEN f CREATE PLOT FILE 

PLOT FLTFARX.CFSrrSX. ELSETSX, CASBDR. PGFDTS, IQKXIMS.SIF . 

FUBV FUCV . GFECT . OES1X/FLOTX2 /MSXLS /LUSPF/ JFLOT/ • 1/ 
S.N.FFXLB S 

PRTMSG FLOTX2 / /FDRMSC S 
END IF S CREATE PLOT FILE 

BWIP $ PHASE XXX DEFORMED AMD CONTOUR P LOTT INC 

RETURN < 

END t 8EDRCVR 
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Appendix D SUBDMAP SUPER3 


SUBEMAP SUP BA .3 CAS BAS, CASBUP. UL, 0L2 , PP 1 , PS 1 , PRL . 0L1 , UH . DLT, 
CAX.CASSBX1, QRC. PJA.OEFNLXX.OBSNLXX, 

PCDB. XYCOB. POSTCDB, POACB/MUQ9I . 

UC, QG.CASBDA.XYCDBPA.FUG BGPSF . OLBRS , 
OUGVl.OPCl.OQCl.OEFlX.OEElX.OSTRl.OCSl. 

OUGV2 . OPG2 , OQG2 , OBP2 , OBS2. 0CTR2.OCS2. 

OES IM , OU 1 0 , QSTR1M, 06 TRIG , BGPSTA, OHACY1 , OGPPB1/ 
APP/APP1/ RSOMLY /CYCLIC /NUHIAT/ ABAC/ DAOPT/ PS / 

ppxlb /caadno / dbs iter/ no nlnr/rkst rain /beta /undoes $ 

t 

TYPE DB 8PT, KPT, CHCWl, CBCM2 .GBOK3 .G80M4 . KMAP, HAPS. BQBX INS, 


DYNAMICS, BGPDTS , CSTMS , KFTS , DIT, BST, GEOK2S , CBOK3S , BCTS . 

■PTS. INDTA , KBLM, KDICT. CPBCT, VBLIX, FOR* , SI LX , 1QEXINX, SCTX. 

BGPOTX. QGI.GBOI44S, PG. BSTL, SPSEL, USBTNL $ 

TYPB DB CASBBX , CAS IX , BSTN, CBOM3BK, GBQM3X , 0LB2. PJ1 . 

ULS.ULS1 S SCRATCH 
I QUALZFSRS 

TYPB PAAM, ,AS. Y.AKSDIS1, BETA $ 

TYPB PAAM, NDDL, I, Y.HIQIQUAL. LMODBS $ 

TYPB PAAM, NDDL, I , H, SBID, XTEMP, LOAD, TBMPLD. DEFORM. HPC . SPC $ 

TYPB PAAM. NDDL. I , M, PI ID. MATH, DYAD, MPLUID, NL99. MLOOP $ 

TYPB PAAM, NDDL , QiAAB , N, K2GG. N2GG , B2CC , P2G . APRCH $ 

TYPB PAAM, NDDL, LOGICAL, N, PSCOUP $ 

$ LOCAL PARAMETERS 

TYPB PAAM, ,CHAA8,N, APP.APPl. SUBDMAP. 'SUPB W .CNTL $ 

TYPB PAAM, .LOGICAL. N, AS ONLY, CYCLIC, CP FDR. NLH EAT. ABftO.PS, NONLNR S 
TYPB PAAM, , I,N,NOS*,IMDDA.NOUPL. INDPLOT.O.PPILB. CAADNO.NOQC.TSBT S 
TYPB PAAM, , I ,N, NOUSBT. NOASBT. NOBSBT. N0CS1T. M0GS8T. NOLSET, NOOEET, NOQGBT $ 
TYPB PAAM, . I ,M, NOASBT , NOSSBT, NOTEBT. MCVSBT, MCA, NOS BT, NOAC , NCKSBT J 
TYPB PAAM, . I ,N , ASPBCTAX . SCRSPECX, WOCOKPX* * I , DBSITBA, BXTRCV1 $ 

$ USBA PARAMETERS 

TYPB PAAM, NDDL, I, Y.PL'IMSC. PORKSG. CADPNT . NOBLOF . NOBLOP . S 1*0 S 
TYPB PAAM, NDDL . I, ¥, CURVPLOT. CURV, MOCOMPS . NUMCUTl . NUHOUT2 . NUMOUT, LSTAN $ 
TYPB PAAM, NDDL, I, Y, SATOPT. SATBLTYP , NOMSGETR, SI, SIM. SIC, OUTOPT, OC , S 1AM $ 
TYPB PAAM.NDOL, I, Y, NINTPTS , DOPT, IRIS. ASPBCTRA . RSPA INT. TABID 5 
TYPB PAAM, NDDL , AS , Y , TABS . BIGBAI . BIGHR2 - BICBA, TINY $ 

TYPB PAAM.NDOL. I. Y. SBSEP . SCRSPBC . IXTRCV $ 
f USSR PAAAM8TBAS APPLIED TO ALL SB 
( AND SAVED ON DATA BASE IN PHASBO 
TYPB PAAM, NDDL, CHAR*. N, ALTABD, HBATSTAT S 
TYPB PAAM, NDDL, I.M. FIXBDB, INABL J 
f SAVED PAAAMXTEAS 
TYPB PAAM, NDDL, LOGICAL, N, SKIPS! $ 

TYPB PAAM. NDDL. I , N, NOUP J FOR POGTRBIG 
I DSC AM) POST ABLATED PARAMETERS 
TYPB PAAM, , I.M, ICYCLIC S 

TYPB PAAM, NDDL. I.M, CP, GIOMU LOADO, POSTU.DBC PATH $ SCRATCH 
TYPB PAAM. NDDL. I, Y. DBCD I AC, POST * 

TYPB PAAM. NDDL, CHAAt, Y, DBCCONV. DBCOVWRT $ 

TYPB PAAM, . CHAAB , N, DBCCONVX $ 

t 

S SBMBLRA FLAGS 

$ 

TYPB PAAM, , I. Y, RWSTRAIN $ 


0UTPUT2 CONTROLS 

THS OUTPUT2 DHAP INSTRUCTIONS IN THIS ALTER HAVE BEXN ENTERED INTO 
TKB MSC/NASTRAN RP ALTER LIBRAAY AS A COURTESY TO PDA/PATRAH AND 
SDRC/I- DBAS USERS AND QUESTIONS AS TO HON IT INTERFACES WITH THE 
PA TRAN OR I -DBAS PROGRAM SHOULD BB DIRECTED TO: 


% 

% 

S 
I 
* 
t 
t 
s 

I 
♦ 

$ 
t 
I 
I 
( 
s 

I 
t 
I 

TYPB FARM. .CKARB. Y.OQC-' YES' .OUC*' YES' ,OBF* 1 YES .OSS* YES' ,OBB*' YES' i 


PATRAN USERS CONTACT ■ 

PDA BNGIN8BRING 
2975 RE DRILL AVI . 
C06TA KBSA, CA 92424 


I -DBAS USERS CONTACT: 
SDftC 

2000 EASTMAN DA. 
MILFORD, OH 4&1S0 


P ARAM. POST. -1, OUTPUTS TKB APPRO FI ATE FILES FOR 

TKB PDA /PATRAN NAS PAT PROGRAM VERSION 2.0. 


P ARAM, POST. - 2, OUTPUTS TKB APPRO PI ATE FILES FOR 

TKB SDAC/I'DBAS DATA LOADER PROGRAM VERSION 3.0. 


TYPB PAAM. .CHARE. Y.Oa**‘YBS' .0GPS-’YB8.0BS«*'Y88'.0CMU*’YBS' C 
TYPB PAAM, . CHARI , T . OCFP» ' YES ' , OUOCORD* '.OCBCM*' * t 
TYPB PAAM, , I , Y , OUNXT1 *11 . 0UNIT2* 12 . CMAXA $ 

TYPB PAAM, , I.M. CTAPE2-0 , OUNITX— 1 $ 



f DESIGN OPTIMIZATION AM) DYNAMIC SENSITIVITY 
TYPB PAAM, . I.M.DROPT t 

s 

( PERFORM PHASE III OPERATIONS 

f 

BQUIVX USET/USBTX/-1 $ 

IP ( APP-NLST- ) BQUIVX USBTNL/UStTX/- l i 
IP ( BXTRCV*0 ) PUTSYSU. 109) $ 

f FOR Cl EXECUTION SO SDA1 IS NOT SKIPPED ON RESTART IN THE EVENT 
$ THAT A NEW R.S. UG IS BEING LOCATED AND NO OTHER CMAMIES AM MADE. 
PAAAML USBTX// 'USBT' ////S.M. NOUSBT// A‘ /S , M. NOASBT/ ' Q' /S , N. NOQSBT $ 
ZP ( NOT ( A PP*' STATICS' AM) N0QSIT*N0A£8T AND N0USET>-1) AMD 
PIXEDB>-0 ) THEN $ CHECK POA R.S. SOLUTION VECTOR 
PAAAML UL// 'TRAILER 1 /1/S. N.NCUL//8.H.NOUL f 
| USB 1X3 FOR EXTERNAL TIP SUFIABLEMENTS WHICH DO NOT NEED UL 
ZP ( MOUL--1 AM) EXTRCVoO ) PAAAML UG// 

'TRAILER' /I/S. N.NCUL//E.M.NOUL S 


ZP ( NOUL<Q > THIN | 


MESSAGE //' DHAP FATAL MESSAGE 90S* ( SUPS AJ ) -'/ 

* THE SOLUTION POE THE RESIDUAL STRUCTUM DOBS NOT EXIST.' { 


EXIT $ 

BNDIP $ 

END IF S CHECK POR A.S. SOLUTION VECTOR 

$ 

PAAAML CASBUP/ / 'TRAILER ' / 1 /S .N . I CASE 1 //S, N. NOCASEUP S 
BQUIVX CASBAS/CASBX/NOCASEUP $ 

IP ( N0CASBUP>-1 ) THEN $ 

APPEND CASBUP, CASERS/CASEX/1 S 
PAAAML CA£SRS//’TRAZLBR'/1/S.N.XCA£B $ 

PAAAML CASSAS/ / 'TRAILER' / 3 /S . N. NAXL $ 

PAAAML CASSAS/ / 'TRAILER' /4 /S . N, IPLT t 

PAAAML CAS BUP//' TRAILER' /3/S, M.KAXLl $ 

PAAAML CASBUP/ /'TRA ILEA '/4/S.H, I PLT1 5 

I CASS ■ I CAS B* IC AS El $ 

MAXL-KXX(KAXL.MAXLl) $ 

IPLT*KAX ( IPLT, IPLTl ) $ 

MODTRL CASBX//1CASB/ICASB /MAXL/IPLT $ 

END IF S NOCASEUP*- 1 

$ 

DBVIBW UGPaUG (WERE SBID* * AM) WILDCAAINTRU1) t 
DBVIBW PUG F* PUG (WERE SBID** AND WILDCARD* TRUE I $ 

DBVIKW QGPaQG I WHERE SBID** AM) WILDCAAD-TRUB) S 

IP I DROPT*S I CMTL* ' ALL ' S 

5EP4 CASBX . PCDB. BMAP. XYCDB, UGF. PUGF, QGF/ 

DRLIST/ 

APP/ ' SBID' /S , N. NOOAALL/S , N. SBID/S, N, NOUPL /OFTL $ 

S NODRALL-- 1 MO DATA RECOVERY REQUESTS 

$ SBID *-l TO INDICATE TO SBDR FIRST TIME INTO LOOP 

S 

S IF SOLS 101 -1S9 THEN DROPT-O { DEFAULT) 

S BUT IN SOLS I0S, 111. OR 112 IF DYNSEN* ' YES' THEN DKJPT.2 

S AND IN SOLS 101. 103, OR 105 IF SENSITY-N THEN DROPT-4 

S IF SOL 200 AND ■ 

S 1. OPT IN* ’NO' THEN DROPT-1 

I 2. OPTIM-rBS' THEN DROPT-2 

$ 3. (NASPRT*- 1 AMD l<OPTBXrT AND OFT EX IT <6 ) OR 

i MOOtDBSrrBRl. 1. NASPRT) *0 THIN DROPT.J 

$ 

IP ( NODRALLaO AM) NOUPL<0 AMD DR0PT<*1 ) RETURN $ 

$ 

S SET FOR ALL SUPBAELEMENTS 

SCASPBCX-LTOKNOTtAPP*' ABIC' AND SCRSPBC**!)) ) $ 

$ DBC CONTROL PARAMETERS 
DBCCONVX* DBCCONV % 

ICYCLIC-LTO I (CYCLIC) $ 

IP < PO6T*0 AM) DROPT<*2 ) DBC. 

CAS EX .//CASBCC' ///////////// /////// 

-1/DBCPATH/S, M.CP/APP1/ ICYCLIC/GBOMU/LOADU/POSTU/ 

D BCDI AG/ DBCCONV/ DBC WWRT/ DBS ITER S 
IP ( POST* -2 AMD (OCEOKa ‘ YES ' OR OGSOH* ’) ) CPI. 

CB0N1 , CSOM2 . /GPt.EQEXrN.GPDT.CSTM, BCPDT.SIL. /0/0/0 S 

$ 

NOGE*LTOI (RSONLY) J 


DO WHILE (ENDOR *« 0) $ SUP BREL SMENT DATA RECOVERY LOOP 
i 

IP ( NOORALL*- 1 AMD DR0PT>1 | THEN C 
S IP DROPT*l IS REQUESTED THEN 

f DO DATA RECOVERY FOB R-S. . EVEN IF MO SUCH REQUESTS EXIST 
SBIDaO 5 
PS1D*0 S 

END DR* -1 $ DATA RECOVERY FOR SEID*0 ONLY 
BQUIVX CASEX /CAS BOB/ ALWAYS % 

BQUIVX UL/ULS/ALWAYS 5 
ILGS S NOORALL*- 1 AND DROPT>l 
S ED ADR DRLIST// 

S.N, ENDOR/S, M. SEID/S.N. FEID/S.N. SEONH/S , N. MOOR/MOCI S 
ZP ( SlIDoO ) THEN $ SET QUALIFIERS 
PSCOUP* FALSE $ 

MLOOP- -1 S 
END IF S 

ZP ( MOOR* - 1 ) THEN $ 

IP ( DROFT»0 ) DELETE /CASlDB, XYCDBDR. , . / t 

PVT PVTB./PVTX/ t UPDATE XPVT -- PVTE POBCXS EXEC ON ASST* 

S INPUTS • UCD - FAMILY OF DOWNSTREAM DISPLACEMENT VECTORS - G-SIZ8 
$ - UL - ABSXMEtL STRUCTURE DISPLACEMENT VECTOR - A -SIZE 

t OUTPUT - UL8 - UPSTREAM BOUNDARY DISPLACEMENT VECTOR OP 
f CURRENT SUPEAELEMENT SBID • A-6XZE 


DBVIBW MAPSP* MAPS ( WHERE SBID** AM) WILDCARD* TRUE ) f 

DBVIBW UGO* UG ( WHERE SBXD*CBDWN AMD WILDCARD- TRUE ) { 

DBVIBW BQSXXNED* BQBXINS ( WHERE PEID-SEDM4 ) { 

IXTRCV1 *IXTRCV | 

IP ( £BZM«<*0 ) EXTRCV1 .0 ( 

DBVIBW EMAPX-EMAP ( WHERE BXTACV. EXTACV1 ) $ 

SEDA EMAPX. CAS »X. PCDB. DRLIST. XTCDB, SLT, ITT, 

MAPSP. UGD. BQBXINSD/ 

ULC.CASBDB, PCDBDR. XYCDBDR/ 

APP/SEID/S.M , NOUP/S , N, NOS ATI/ /S, N. NOOUT/ 

E. M. NOPLOT/S. M.NOKYPLOT/ SE ID ' /NCUL f 
IP ( MOT (APP* 'STATICS' AM) NOQSET.NOASET) AM) 

MOOO*-1 AM) PIXED»**0 AMD SBIDoO ) TUN $ 

MESSAGE //' DMAP IWORMATION MESSAGE 9012 (SUPBA3) -V 

' THE DISPLACEMENTS IN DO WNS TREAM SUPBRBLMBNT ' / 
SBDWN/' DO NOT EXIST. f 

EXIT $ 

BNDIP $ NOGOa-1 AM) PIXBDB**0 

ZP ( SBID.0 ) BQUIVX UL/ULS/ALWAYS $ 

EM) IF S N0DR>-1 

8M>XP S ELSE NOORALL*- 1 AM) DROPT>l 


D1 



IP ( NOOR.-1 AND 

I APPo'MIO ■ OA WOTISBSBF-l AND NOUP.-H J ) THEN 5 

IP (HEATSTAT- 'YES' ) SQUrvX CAS SS /CASE DA /ALWAYS $ 

CAI.D SffTg CASES/ /SI I D/Pl IO/S , MTBHP/S , K2GG/S . N2CJC/S , 82 GC /& , MPC / 
S . SPC/S , LOAD/S , DEFORM /£ . TKKPLD/S , P2C/S , DYRO/S HETH/ 

S , KPLUID $ 

IP ( DBOPTvO ) THEN S 


GP3 

TA1 


DELETE 

DSLBTB 

DELETE 

D1LBTS 


DELETE 

INDIE $ 


/UG.QC.0LB2. ./ S 

/Poe, Bcpsp.oocvi.opci.ogcx/ $ 
/0BF1X.OBS1X.OSTA1.0CS1.0OCV2/ $ 
/0PG2 , OQC2 . 0BF2 . 01*2. 0STR2 / $ 
/0GS2 , 0BS1M, 0181 C. OSTIUM. 06TP1G/ 
/BCFSTR.OMRCYl.OCFPBl. . / $ 


vquivx usrr/us*TX/-i t 

IP i APP- NLST' AM) SEID.0 | sourvx USBTNL/USBTX/- 1 J 
CALL PMLIJSBT usrrx / / s , nor sit/ s , nobs«t/ s , nocsit/ s . nogsht / * , nolsbt / 
£ , M00S8T/S, NOQSBT/ 6 , MORS IT /£ . NQSSBT/S . NOTSIT/ 

S . NOVSET/ 8 . NCA /S . NOSBT/S . NO AC/S . N0KS1T/ 

s.nousbt $ 


USET HAS NOT 


IP { MOUSBT.-l ) THIN J 
THIS CAN HAPPEN FOR PARAK, PXXHDB, - 1 AND A.S. 

BBBN CBN SPA TIO IN SUCK OPERATION 

PRTPAAM //4410/ 'DKAP' //SUBDKAP $ 

BLSB S NOUSXT* - 1 

IP ( SlIDoO ) ML99.0 $ MS ST QUALIFIER IP NOT R.fi. 

FOR PROPER NPTS IN NONLINEAR TRANS I8NT 

PARAKL CAS SDR// ’ DTI ' /-1/3S//S. N, SPCFOR 5 
IP S PC FORCES ARB REQUESTED THEN SPCFOR.- 1. OTHERWISE 0 
PARAKL CASEDR// 'DTI'/ -1/1*7/ /B.N.CPFOR < 

IP CPPORCI S AM RBQU1STBD THBN GPPOR—1. OTHERWISE 0 

NOQG.LTOI ( NOT (OALt SPCFOR, GPPOR) OR (DROPT.1 AND DROFT<»4) | ) $ 
PARAKL CASEDR/ / 'DTI '/-1/17Q/ /S.N. BSE $ 

IP BSB AM REQUESTED THIN ESB--1. OTHERWISE 
CPPDR. (ORL {GPPOR, BSE) > $ 

IP | APP. ‘ CYC ’ AND APEX. -STATICS * ) THEN $ 

CYCLIC! CASEDR, GB0K3S. GB0M4S. DIT. / 

DBI . 010*0 BE .DB3.DB4.DRi. CAS HR K . BACK/ 

PK1/PM2 / ' STATICS ' /PH4 /PNS/PHt S 
EQUIVX CAS ESK /CASEDR/ ALWAYS $ 

BQUIVX CBOM3BM/CECM3X /ALWAYS 5 
KPYAD FORE. BACK. /PB $ 

BLSB IP { CYCLIC AND (APP.'REIG' OR APP. FREQRESP’ ) ) THEN J 
BQUIVX CAS 1BK1 /CASEDR /ALWAYS $ 

BLSB $ APP. ' CYC - AND APPl* 'STATICS ‘ 

BQUIVX CBO*OS/CBO»OX/ ALWAYS J 

IWUP $ APP* ' CYC' AND APP1. STATICS ' 

IP ( APP.'NLST' AND £EID*0 ) THIN * 

CORMCTION POR OLOAD OUTPUT - CRX IS ACTUALLY COKB IN NLSTATIC 


TSBT.GHTSYS (TSBT. 10 ) $ 

MATHOD CHOK3S.CPLS, UC, . , /CB0K3SN. /1B/TBBT $ 

CBOK3SN , BQ8XINS . GBOM2S/SLTN , BTTN/123 /123/ 123 S 
. .BCTS. EPTS. BGPDTS. SILS.ETTN.CSTKS, / 

SSTN. .GBIN.CPSCTN. / 

LUSBT/123/ 12 3/ 1/123/123 $ 

ELSE IP < NONLNR AND SBID.0 ) THBN S 
BQUIVX HSTL/BSTN/ ALWAYS 5 

BLSB $ APP*' STATICS ’ AND NLHBAT 
BQUIVX BST/BSTN /ALWAYS 5 

BNDIP $ BLSB APP. ’ STATICS ' AND NLHBAT 

IP < APP. STATICS' OR APP.'REIG' OR APP.'NLST' OR 

(APP.'CYC' AND APP1. 'STATICS' ) ) HOC OKPX. NOCOMPS J 

RSPBCTRX. LTOI ( NOT ( APP* ' TRAN MSP ' AND 

RSPBCTRA..0 AND NOT (NONLNR ) ) | $ 

WRITE END- OF- DATA IP CURRENT 0UNIT2 IS DIFFERENT FROM 
PREVIOUS owm (ounitx) and this is not thi 
FIRST TIKE THROUGH THE DATA RECOVERY LOOP. 

IP { PQST«0 AND OUNITX<>OUNIT2 AND 

OLWITX< >- 1 ) OUTPUT2 / / - 9 /OUNITX/ /OKAXR 5 
PUT LABEL ON CURRENT UNIT ONLY IP POST. -I 

IP ( POST. - 1 AND OUN ITX< VOUNIT2 ) OTA PE 2. -1 5 

IP ( POST. -2 AND (OGBOK* 'YES’ OR OCEQH-' ') ) THEN S 
IP GEOMETRY DATA BLOCKS ARE REQUESTED THBN OCWITl AND 
OLNm AM THE SAME UNIT. 

0UNIT1.0UNIT2 $ 

0LTPUT2 CSTM.GPL.GPDT, BPT, KPT//O/OUNITI//OKAXR $ 

0UTPUT2 GBOK2 . G BOM3 . GB0M4 , , / / 0 / 0UNIT1 / /OMAXR $ 

BNDIP S POST* -2 AND (OGEOK* ' YBS ‘ OR OGBOK. ' '} 


I OOGCORD. ‘ GLOBAL ' S 


PC.PJ/DPJ//-1. S 
DPJ.CRX, /DPJX S 
PJR , DPJX/PJl // -1 , $ 
OL2/OLB2/ ALWAYS | 


ADD 
KPYAD 
ADO 
BQUIVX 
BLSB $ 

CALL SBDISP COAT. GOAQ. LOO. KOO. LAO. POS. UOX . GPLS , USRTX , SILS, 
PJ. CASEDR, ULS. BDT. YS. 04. PSS . KPS. KSS QR , CMPHO. OL2. 
CMLAMA.CWPHA.PP1 .PS1.HAR.MBA. XYCD6DR. CRX. 

BACK. FB.QRG/ 

UC.QC.0LB2.ULS1.PJ1/ 

P IX EDB/NOQSET/ HOTS BT/ IRES /APP /APP1 /NOOSBT/SBID/ 
NOOUT /CARDNO/PPILB/ XNML/ALTRBD/NORSET/NOQG/ 
NLHBAT / DROPT/NONLNR t 

END IF J 

IP < SBID.0 ) BQUIVX OLB2/OL MS /ALWAYS $ 

IP \ POBT.Q AMD (APP. ' MIC’ OR APP.'CXCGBN' I ANV 
DROPToJ ) DBC. 


.// 


'LAMA'//////////////////// 

-1/DBCPATK/* . N. CP/APP1/ICYCLIC/GBOMU/LOADU/POSTU/ 
OBCDIAC/npCCOMV/DBCOVWirr/DBSITBR * 

IP ( (NODRALL. - 1 Ate DAO FT -4 ) OR DROPT.2 ) THBN $ 

* IP S1NSITY *N I* RBOUUTBD IH SOLS 101. 103. A*K> 10S 

S AMD NO DR REQUESTED THEM RETURN 

I IP PAAAM.DYXSm.YE8 IM SOLS 10B. 1«. AND 112 THEN RETURN 

RETURN S 

ELSE IP < DR0PT.5 ) THEM * 

< SUPERELBMBNT SmsiTIVITY - IGNORE DATA RECOVERY ASQU1STS 

NOOUT.- 1 $ 

END IP S 

S 

IP I RSOM.Y > BQUIVX PCDB/ PC DBDR/ ALWAYS S 

$ THIS IS N BIDED IP R.S.-ONLY NOOIL AMD NO SIP LOT OR 
S SBUPPLOT CARDS AM PRESENT 
{ 

IP ( APP. STATICS ' AND NLHBAT I THEN $ 


IP ( POST* - 1 AND OUGCORD. 

IP < N00UT>-1 ) THIN $ 

CALL SBDRCVR, 

UG, QG . BGPDTS . EQBXINS , CSTMS , CASEDR. KPTS . DIT 
BTT. 0LB2 , P.T1 . BSTN. XYCDBDR , GBOK2S , GBO»OX . POSTCDB . 
ECTS, GPLS. BPTS.SILS , INDTA, KBLN, EDICT. GPICT, VBL3M. 
FORC8.XYCDB. K5. PCDBOR, USITX.SLT. UH.OL1 . PPl DLT 
PRL. SPSEL . DYNAMICS . CRX . QGE , OB PNLXX . OBSNLXX/ 

MUCNI . PUC , OUGV1 . OPG 1 . OQG1 . OIF 1 X . OBS IX . OSTR1 . 0GS1 
OUGV2.OPG2.OQG2.OBP2, OES 2. 0STR2.0GS2, 

0BS1M.0BS1G, OS TRIM. OSTRIG, BCPSTR . EGPSP 
0NRCY1. 0CPPB1/ 

GRDPNT/APP/APP1 /NOCOKPX /CURVPLOT/S , PPILE/NUWOUri / 
NIW0UT2 /BIGBR1/BICBR2 /NUHOUT/BIGSR/SRTOPT/LSTRN/ 

S RTBLTYP / C URV / OUTOPT/ (X3 /N INTPTS / S 1 K/ S 10 / S 1AK/ 

S LAG/ DOPT/ TINY /NOB LOP /NOB LOP/ TABS /SB ID/S l/S.CARDNO/ 
PDRMSG / SCRSPBCX / RS P8CTRX / A£ PRINT /TABID / INRIL / CPPDR / 
NLHBAT/ABRO/ ICYCUC /PS /CS0MU/L0ADU/P06TU/ 
OBCDIAG/DBCCONVX/POST/S , CP/ DBCOVWRT/ 

S . 0UNXT2/S. OTAPB2/OQG/CKJC/OSP/OBS/OBS/ 

OCHP/ OGPS / OBSB /OUMU /OGPP / OUOCORD / D ES IT1R / RMSTRA IN / 

B8TA/LHOD8S $ 

$ 

01BIITX.0UNIT2 $ SAVB UNIT NLMB1R POR COMPARISON ABOVE 
BNDIP S NOOUT >- 1 
BNDIP £ BLSB NOUS IT. - 1 

5 N0DR>-1 AND ( APP. » ' RBIG ' OR NOTISBSBP.l AND NOUP.-1J 
END DO $ SUPIRBLBMBNT DATA RECOVERY LOOP 
$ WRITE END-OP- DATA ON FORTRAN UNIT OUNXT2 
IP l POST.O | OUTPUT2 //- !*/0UNIT2/ /OKAXR J 

SBID.-1 s INITIALIZE 

DO WHILE ( NOUPL.-l AND BNDPLOTo-1 | S 

$ PHASE IV - UPSTREAM DEFORMED PLOTS 

t 

PLVCDR BKAP. DRLIST/ /S ,N, BNDPLOT/S.N, SBID/S.N, PBID S 
IP ( BNDPLCTo-1 | THBN 5 

DBVItW SILSP . SILS (W8M PBID* * ) 5 

DBVIB W CPLSP a GPLS IWKBM PBID. * ) 5 

$ COLLECT UPSTREAM DISPLACEMENT PLOT VICTORS 

PLTVBC EMAP, DRLIST. PCDB. 1QEXINX, SI LX , SILSP, GPLSP, PUGP/ 

POCX. PCDBX/SBID/'SIID' $ 

PARAKL PCDBX// 'PRESENCE ' / ///8, N. JPLOT $ 

IP ( JPLOT>. 0 ) THEN $ 

PLTSBT PCDBX. SQIXINX. BCTX/PLTXY. PL TPARY . GPS ITS Y, KLTSETS Y/ 

S.M.NSILS/S.N. JPLOT $ 

PRTWSO PLTXY//PDRMBG | 

IP ( JPLOT. «0 ) THEM $ 

PLOT PLTPARY , GPSBTSY , ELTSBTSY, CASEX , BGPOTX . BQEXINX 

SILX. PUGX.POCX. . / 

PLOTY2/NSILS/0/JPLOT/-1/S.N.PFILB $ 

PRTNSG PLOTY2/ /PDRMSG $ 
tt®IF S JPLOT.. 0 
BNDIP $ JPLOT*. 0 
BNDIP S BNDPLCTo- 1 

BNDDO $ PHASE IV - UPSTREAM DEFORMED PLOTS 
i 

RETURN $ 

EW> $ SUPER 3 
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